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Zusammenfassung 


Die Arbeitsweise yon kreisringférmigen Ultraschallpfeifen wird beschrieben; bei niede- 
rem Druck: Selbsterregung des Luftstrahls, Angabe der Bedingungen fiir guten Wirkungs- 
grad (30%) und Erzeugung reiner Tone; bei hohem Druck: Schwingungserregung des 
Luftstrahls durch Reflexion an der Basis, einem sehr komplexen Tongemisch iiberlagerter 
Schneidenton, Ultraschallkomponenten sehr hoher Frequenz, aber verschlechterter Wir- 
kungsgrad. Berechnung der Pfeife. Uber Druckkammer angekoppelter Trichter nicht zu 
groBer Offnung. 

Strahlergruppen aus Pfeifen. Notwendigkeit der Synchronisation: Direkte Kopplung oder 
solche iiber Resonanzrohr. Synchronisation mit zusatzlichem Hohlraum (,,akustische 
Triode“), angewandt bei Pfeifen nach Lrvavassrur. Synchronisation mit Pfeifen, Sirenen 
und Lautsprechern. Vor- und Nachteile gegentiber Ultraschallsirenen. 


Summary 

Toroidal whistles derived from the police whistle are described. Their operation at low 
pressure —oscillation of the air jet produced by the emitted sound, conditions to obtain 
high efficiency (30%) and pure sinusoidal tone without harmonics is discussed. Their 
operation at high pressure — oscillations produced by the air jet returning to strike at its 
base, edge sound superposed on the other components of the emitted complex sound. The 
equation of whistles and calculation of their theoretical efficiency are given. Contradiction 
between the theory of horns and the experimental results: no cut off frequency. Practical 
calculation of annular exponential horns for emission of plane waves. Advantages and 
disadvantages of whistles and of ultrasonic sirens. Applications are described. 


Introduction Une autre application trés importante des sons 


Les sons aigus, a partir d’environ 2kHz, et les &t des ultra-sons dans l’air et dans les gaz, est 


ultra-sons émis dans l’air, provoquent une coagula- 
tion rapide et la précipitation de brouillards, fumées, 


de certaines réactions 
1 


l’accélération considérable 


chimiques, telles que la combustion! et la catalyse 


poussiéres, — lorsque leur intensité dépasse 150 dB 
(soit 0,1 W/cm?). 


* Communication au 2™° Congrés International 
d’Acoustique 4 Cambridge, Mass., Etats-Unis, 17 au 
23 Juin 1956. 


[1], [2], [3], [4]. 


1 Essais faits au «Laboratoire National d’Essais des 
Arts et Métiers», par M. Cxarpennet, sous la direc- 
tion du Prof. Veron, pour la société des Chaudiéres 
Babkock et Pillard. 
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De tels sons sont produits par des sirénes tournant 
a grande vitesse ou, beaucoup plus facilement et 
plus économiquement, par des sifflets toroidaux. 


Dans les sifflets industriels de grande puissance, on — 


utilise le méme principe que dans un sifflet de police. 


Fig. 1. Sifflet de police. 


Rappelons qu’un sifflet de police (Fig.1) est 
constitué par une cavité fonctionnant comme un 
résonateur de Helmholtz. Un jet d’air oscillant, émis 
par une buse, est dirigé vers le bord bien aiguisé de 
la cavité resonnante [5]. L’expérience montre qu’il 
existe trois modes de fonctionnement possibles d’un 
tel sifflet, ces modes de fonctionnement pouvant 
d’ailleurs coexister: 


1. Le fonctionnement «normal», a basse 
pression [6], Fig. 2 


Un sifflet de police ordinaire, alimenté sous sa 
pression normale de 40 gr/cm? émet un son trés pur, 
dépourvu d’harmoniques, et son rendement est tres 
bon: de l’ordre de 30%. On croyait autrefois que la 
forme du couteau était sans importance; en réalité, 


Fig. 2. Fonctionnement a basse 
pression. 


elle n’influe pas sur la fréquence, mais influe beau- 
coup sur la puissance et le rendement: le couteau 
doit étre aussi mince et aussi bien aiguisé que pos- 
sible. La puissance acoustique d’un tel sifflet ordi- 
naire est relativement faible (1,15 W environ pour 
une puissance pneumatique de 4 W environ fournie 
par l’agent de police). 


2. Le fonctionnement «en siréne statique», a haute 
pression, Fig. 3 (cas d’un compresseur fournis- 
sant une pression de plusieurs kg/cm?) 


Le jet d’air a vitesse supersonique, contourne 
intérieurement la cavité resonnante et revient frap- 
per a sa base, ce qui provoque périodiquement une 
trés brusque déviation du jet. Le jet d’air est ainsi 
périodiquement interrompu au bord de la cavité 
résonnante. Ceci rappelle le fonctionnement d’une 
siréne ordinaire dans laquelle le jet d’air est périodi- 
quement interrompu par un disque ou un tambour 
perforé tournant. Mais ici, c’est le jet qui se déplace 
par rapport a l’obstacle. 
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Le son émis n’est pas pur, mais riche en harmont- 
ques (ce qui est favorable pour les applications 
industrielles) . 


N 
\ 


Fig. 3. Fonctionnement 4 haute 
pression. 


3. «Son de couteau» 


Ce son est produit par une oscillation trés rapide 
du jet d’air sous l’effet des tourbillons qui l’accom- 
pagnent. Ce phénomeéne n’existe pas dans le sifflets 
a trés basse pression; il apparait cependant, bien 
avant que la vitesse du jet devienne supersonique. 

Le son de couteau donne naissance aux compo- 
santes ultra-sonores intenses de fréquences tres éle- 
vées. Les fréquences de ces composantes sont prati- 
quement indépendantes de la fréquence «fondamen- 
tale» du sifflet (ce ne sont pas des harmoniques) . 


1. Adaptation d’un pavillon exponentiel 
a un sifflet 


Sur un sifflet ordinaire il est impossible d’adapter 
un pavillon: sa présence arréte instantanément le 
fonctionnement du sifflet?. Or, un pavillon est néces- 
saire pour obtenir des ondes planes (directivité, et 
meilleure résonance dans le cas des traitements in- 
dustriels). Et, contrairement aux simples réflecteurs 
(miroirs) permettant tout au plus d’obtenir une 
directivité accrue, les pavillons améliorent con- 
sidérablement le rendement. 

Dans notre Laboratoire, au Centre de Recherches 
Scientifiques Industrielles et Maritimes de Marseille, 
M. R. Levavasseur® a réussi a tourner cette diffi- 
culté. I] utilise une cavité resonnante auxiliaire, 
sensiblement de mémes dimensions, mais un peu 
plus ouverte que la cavité «sifflante», Fig. 4. Sur 
cette figure S est la cavité sifflante et A la caviteé 
auxiliaire. I] devient ainsi possible de munir les 
sifflets d’un pavillon, par exemple, d’un pavillon 
exponentiel. 

Récemment, R. LevavassEur® a apporté un autre 
perfectionnement trés important aux sifflets: il a 
supprimé le jet d’air qui accompagne le son émis 
par le sifflet. Ce jet d’air était trés nuisible car il 
diluait la fumée a coaguler et il la refroidissait, ce 
qui diminuait le tirage des cheminées. 


3 


> Un sifflet de Galton ordinaire, un sifflet de Loco- 
motive, un sifflet de Hartmann peuvent étre munis d’un 
réflecteur (miroir) parabolique, mais pas d’un pavillon 
(tel qu’un pavillon exponentiel) . 

3 Brevets du Centre National de la Recherche Scien- 
tifique. 
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Dans un modele récent du sifflet Levavasseur, 
le jet d’air dés qu'il a dépassé le couteau, est capté, 
comme le montre la Fig. 5, par une fente dans la 
paroi du pavillon. Ainsi l’air est séparé du son et 
renvoyé vers l’arriére par une canalisation spéciale. 


Fig. 4. Sifflet a cavité 


auxiliaire. tion du jet d’air. 


2. Sifflets toroidaux 


Les sifflets industriels de grande puissance sont 
des sifflets toroidaux Fig. 6. Leur cavité sifflante et 
leur cavité auxiliaire sont toutes les deux toroidales, 
tandis que le jet d’air est annulaire. Leur pavillon 
exponentiel est annulaire, semblable a celui des 
sirenes industrielles. Sur la Fig. 7 on voit le plus 
grand et le plus petit sifflets réalisés dans notre 
Laboratoire. Le plus grand («sifflet No. 5») est un 
sifflet a basse pression, débitant environ 200 litres 
d’air par seconde, — mais cet air est séparé du son 
et renvoyé en arriére. I] émet un son de 2,5kHz 
(c’est la fréquence d’un sifflet de police ordinaire*; 
accompagné d'un ultra-son beaucoup moins intense 
de 28 kHz environ. Sa puissance acoustique est de 
1kW environ. Ce sifflet coagule bien la poussiére 
humide des usines de ciment, mais reste sans effet 
lorsque cette poussiére est séche. Le diamétre de 
Youverture du pavillon est de 60 cm. 


Fig. 6. Principe du sifflet toroidal. 


Le petit sifflet qu’on voit sur la Fig.7 est un 
sifflet ultra-sonore de laboratoire il émet un mélange 


4 Mais beaucoup plus efficace ... pour disperser une 
foule par exemple. 


V. GAVREAU: GENERATEURS PNEUMATIQUES D’ULTRA-SONS 


Figs5. Sifflet avec sépara- 


123 


Fig. 7. Grand sifflet industriel et petit sifflet de labo- 
ratoire. 


de deux sons tous les deux trés intenses: un son de 
22 kHz et un son de 460 kHz. Sa puissance est de 
quelques watts acoustiques seulement, ce qui suffit 
Wailleurs pour les essais de laboratoire et pour la 
signalisation °. 


Fig. 8. Différents types de sifflets toroidaux. 


Sur la Fig. 8 on voit toute une gamme de sifflets 
toroidaux de différentes dimensions et de différents 
types, réalisés dans notre Laboratoire. Le plus grand, 
a droite est le sifflet industriel No. 4 a basse pression, 
avec séparation complete du son et du jet d’air. Il 
émet un son riche en composantes ultra sonores, 
convenant bien pour le dépoussiérage et la pré- 
cipitation de fumées (fumée de tabac précipitée en 
2 secondes) les fréquences des composantes les plus 
intenses sont: 2,72kHz et 38kHz. Sa puissance 
acoustique est de 262 W environ. Il consomme 
108 litre par seconde; la pression d’alimentation est 
de 400 gr/cm?. Ce sifflet est utilisé dans la premiére 


5 Les sifflets de commerce, dits «ultrasonores» 
(«pour le dressage de chiens») émettent un son par- 
faitement audible dont la fréquence ne dépasse jamais 


14kHz. 
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installation industrielle, dans une usine de blancs de 
zinc. 
Sur la méme Fig. 8, a l’extréme gauche, on voit 


le «sifflet en bronze» émettant un son de 9,4 kHz’ 


+150 kHz avec une puissance acoustique de 60 W 
environ. C’est un sifflet 4 haute pression: 2,3 kg/cm?; 
il ne consomme que 7 litres par seconde environ 
(tous les débits sont indiqués, ramenés a la pression 
atmosphérique). Des sifflets de ce type, un peu 
modifiés, a trou central, sont utilisés dans des bru- 
leurs industriels de charbon pulérisé et de mazout. 


Fig.9. Rides de sable produites par le petit sifflet de 
laboratoire (la graduation a droite est en cen- 
timétres) . 


Fig. 10. Visualisation strioscopique du son émis par un 
grand sifflet. 


Le plus petit des sifflets de cette figure est le méme 
que sur la Fig. 7. Sur la Fig. 9 on voit les rides de 
sable produites par ce petit sifflet. Tandis que sur 
la Fig. 10 on voit l’image strioscopique (donnée par 
un miroir parabolique de 30cm de diamétre) des 
ondes sonores planes émises par le grand sifflet a 
haute pression (Fig. 8 au centre), dans le voisinage 
immédiat de son pavillon également de 30cm de 
diamétre. 
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3. Caleul des sifflets toroidaux 


Le fonctionnement des sifflets toroidaux est tout 
a fait analogue 4 celui d’un sifflet de police. La fré- 
quence d’un sifflet de police est donnée par la méme 
formule que celle d’un résonateur de Helmholtz 


oo Cee 
I- V4 


Ou c est la vitesse de propagation du son dans I’air, 
V je volume de la cavité resonnante et y un coef- 
ficient qui dépent de la forme et des dimensions de 
Vouverture par laquelle la cavité résonnante com- 
munique avec l’extérieur. 

Or, lexpérience montre que, dans le cas des sif- 
flets cylindriques (sifflets de police), la fréquence ne 
dépend pas de la longueur L de la cavité résonnante 
cylindrique. Son volume V étant proportionnel a 
cette longueur L, le coefficient y doit également étre 
proportionnel a L 


(1) 


= lle 


Cc EE , /L 
eS ee ee fp ve 
f Ve V 


Et, dans le cas d’un sifflet cylindrique (sifflet de 
police), on trouve ainsi: 


&L ; LL eee 
ee hee 207 Leap «Ae 
rece Vck jn@d a 


(2) 


Donc 


=o) 


Ou d est le diamétre de la base, L la longueur de 
la cavité cylindrique et K une constante, facile a 
déterminer expérimentalement. 

Dans le cas d’un sifflet toroidal a@ cavité sifflante 
«intérieure» (tel que celui représenté sur la Fig. 6), 
en désignant par D le diamétre moyen du tore 
(diamétre du cercle décrit par le centre de la section 
qui engendre le tore) et par d le diamétre de la 
section, on a 


x D @ 


(ay et L=a(D + d) = longueur de l’ouverture, 


pak 14+ _2K1/D+d_K 
x D d? |x D d? d 


Dans le cas d’un sifflet toroidal a cavité sifflante 
«extérieure» (sifflet analogue a celui représenté sur 
la Fig. 6, mais dans lequel la cavité sifflante se 
trouve a l’extérieur et la cavité auxiliaire a linté- 
rieur), ona 


eta € et L = x(D — d) = longueur de l’ ouverture, 


(3b) 


D-4_K1/D~d 
D@ d Dems 


donc 
aay MULE) 2k 
xz D d? 


|x 


——_ eo 
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Un grand nombre de mesures faites sur des sifflets 
de police et sur des sifflets toroidaux de différentes 
dimensions, ont donné des résultats extrémement 
concordants: 


Pour les sifflets sans cavité auxiliaire 
K=51 000 environ. 


Pour les sifflets avec cavité auxiliaire 
K =70 350 environ. 


Toutes les dimensions, c’est a dire D et d, étant 
exprimées en millimétres. 


4. Caleul des pavillons annulaires émettant 
des ondes planes 


Pour qu’un sifflet toroidal émette des ondes par- 
faitement planes on calcule son pavillon exponentiel 
annulaire de la fagon suivante: soient (Fig.11) R, 
et Ry les rayons des deux cercles limitant |’entrée 
annulaire du pavillon (c’est a dire les rayons des 
bords de la cavité sifflante et de la cavité auxiliaire) . 

Soit R le rayon du cercle qui partage cette aire 
annulaire en deux parties égales 


10° 2 10° 


Fig. 12. Directivité obtenue avec un pavillon exponen- 
tiel correctement calculé. 
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(Rp? — R?) = a(R? — R,?) 


rea 
R= ee (4) 


d’ot 


On considére le cylindre passant par ce cercle de 
rayon R et ayant ses génératrices paralléles a l’axe 
du sifflet et on calcule les surfaces limitant le pavil- 
lon annulaire exponentiel de telle fagon que les 
sections comprises entre chacune de ces surfaces et 
ce cylindre de rayon R soient toujours égales. 

Désignons par y le rayon de Ja surface exponen- 
tielle intérieure et par z le rayon de la surface 
exponentielle extérieure dans une section se trouvant 
a une distance x de l’entrée du pavillon, on a 


S eee) 
SAT pie Oe epee | Re an ae 
y y/ Saas Zz V a oe (5) 


i359 alee (6) 


avec 


So étant la surface totale de la section annulaire a 
Ventrée du pavillon. 

Théoriquement [7], un tel pavillon devrait avoir 
une «fréquence de coupure» 


fy=2710 m (7) 


(en Hz) et ne devrait pas transmettre les sons de 
fréquence inférieure a fy. En fait, l’expérience mon- 
tre qu’il n’en est rien: la «fréquence de coupure» 
n’existe pas. Et les constructeurs de sirénes ultra- 
sonores utilisant des pavillons trop largement ou- 
verts, conformes 4a cette théorie erronnée, obtiennent 
de trés mauvais rendements: en fait, c’est juste le 
contraire qu’on doit faire, — utiliser des pavillons 
trés allongés dans le cas des ultra-sons, pour éviter 
le décollement du faisceau ultra-sonore des parois 
du pavillon. 

D’ailleurs, le rdle du pavillon exponentiel n’est 
pas du tout le méme dans le cas des haut parleurs et 
dans le cas des sirénes et sifflets: dans le cas d’un 
haut-parleur, le pavillon réalise une «adaptation des 
impédances», permettant a une petite membrane 
métallique de mettre en mouvement un grand vo- 
lume de gaz. Dans le cas des sifflets le son étant 
engendré dans le gaz lui méme, on n’a pas besoin 
d’une telle adaptation. Le réle du pavillon est alors 
complétement différent: il transforme en son ou en 
ultrason [8] un jet d’air simplement interrompu dont 
les troncons auraient tendance a se recoller, freinés 
par l’air ambiant. 

On calcule le pavillon exponentiel en se donnant 
la hauteur de la basse centrale; connaissant la 
diamétre de la base on en déduit la valeur de m. 
Comme exemple de valeur de m, utilisée pour le 
caleul d’un sifflet de 10 kHz, on peut indiquer 
m =0,0491 dans l’un des sifflets réalisés. 
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5. Rendement et puissance des sifflets 


Dans un article publié précédemment [6] nous 


avons établi les formules théoriques donnant le ren-: 


dement et la puissance des sifflets, et montré qu’il 
était possible d’obtenir pratiquement des rendements 
aussi bons que ceux des sirénes. Nous nous con- 
tenterons ici de rappeler trés briévement les résultats 
de cette étude: 

Le rendement théorique d’un sifflet est donné par 
la formule: 


2 qo\" 
[pees ete) Vr 
Pp a? 6°06 | Wo 
Ua pee 3 2a," ) (8) 
Q qo (V+ ia) 4 V7 
; Wo 


ot Py = 1 013 000 dyne/cm? est la pression atmosphé- 
rique, 
e@ =0,001205 gr/cm® est la masse spécifique de 
Lair a 20° C, 
c¢ =34000 cm/s est la vitesse de propagation du 
son dans |’air, 
qo est le débit de l’air en cm?/s (qq = 2000 cm?/s 
pour un sifflet de police ordinaire), 
a est l’épaisseur du jet (a=0,15cm pour un 
sifflet de police ordinaire), 
b est la largeur de la fente (hauteur du volume 
résonnant cylindrique) (b=1,56 cm), 
6 est la hauteur de l’ouverture (6 =0,73), 
f est la fréquence (f = 2900 Hz donc m = 18 230 
dans le cas du sifflet de police étudié), 
Vest le volume de la cavité résonnante (V = 4,15 
cm? dans le cas du sifflet de police étudié) . 


Dans les cas du sifflet de police étudié, on trouve 
comme valeur théorique maximum du rendement 


n=21%, 


et, expérimentalement, on a trouvé 7=13,25%. 
Dans le cas d’un petit sifflet de type «simplifié» 
(émettant un son de 3900 Hz), on a trouve, expéri- 
mentalement, un rendement 7 = 28,4%. 

La puissance acoustique d'un sifflet, par centi- 
métre de longueur de couteau est d’autant plus 
grande que la pression d’alimentation est grande. 
Pour un bon sifflet de police, p=40 gr/cm? et on 
trouve: 


_11BW 


= = 0,725 W/cm. 
1,56 cm 


Wo 


Pour le sifflet Levavasseur No. 4, p= 400 gr/cm? et 
on a trouvé: 

262 W 
ux 13cm 


Wo = 6,42 W/cm . 


Lintensité du son émis dépasse toujours 150 dB 
dans le voisinage du sifflet par exemple pour le 
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sifflet No. 4 de diamétre = 30 cm elle est 


T= ——*— = 0,37 W/cm? soit 156 dB. 


Tw 2 
eg 


Et dans une «colonne de traitement», dans les condi- 
tions de résonance, elle atteint facilement 160 dB °. 


6. Batteries de sifflets («panneaux sifflants») 


Dans le cas des sifflets, tout comme dans le cas 
des sirénes, on est obligé d’utiliser des fréquences 
moins élevées lorsqu’on veut émettre des puissances 
plus grandes. Le son des grands sifflets industriels 
est parfaitement audible et on ne peut pas les utiliser 
a lair libre, par exemple pour dissiper le brouillard: 
Vintensité du son est absolument insupportable. 

C’est pourquoi, pour pouvoir émettre de grandes 
puissances acoustiques aux fréquences élevées, nous 
avons cherché a grouper en batteries’ des sifflets 
Levavasseur de dimensions moyennes. 

Pour fonctionner avec un bon rendement, sans 
court-circuits acoustiques, et émettre des ondes 
planes, directives, tous les élements d’un panneau 
sifflant doivent étre parfaitement synchronisés en 
fréquence et en phase. 

Les sifflets sont des oscillateurs acoustiques et, 
comme tous les oscillateurs, ils peuvent étre syn- 
chronisés par couplage, comme le montre la Fig. 13. 
Ce couplage peut étre réalisé de deux fagons: soit en 
prévoyant une petite ouverture faisant communiquer 
les cavités résonnantes des deux sifflets, soit en 
établissant une communication entre ces cavités 
résonnantes au moyen de tubes de longueur con- 
venable (accordés a la fréquence des sifflets). 


rs 


t 4 f \ 


er 


| 
Fig. 13. Couplage direct et indirect de sifflets. 


Un tel couplage permet la réalisation de sifflets 
Levavasseur a deux cavités sifflantes concentriques, 
permettant d’obtenir des intensités et des puissances 


6 Les puissances acoustiques et les rendements in- 
yraisemblables indiqués par certains constructeurs de 
sirénes, s’expliquent par des mesures erronnées. 

7 Brevet C.N.R.S., Invention de MM. VY. Gavrzau et 
A. Caraora. 
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plus grandes que les sifflets simples. Toutefois, nous 
avons abandonné la construction de tels sifflets 
multiples, car nous avons trouvé un procédé de syn- 
chronisation bien meilleur, beaucoup plus efficace, 
qu'un simple couplage: c’est la synchronisation par 
«buse auxiliaire». 


Fig. 15. Sifflet ordinaire a buse auxiliaire. 


Un sifflet 4 buse auxiliaire, Fig. 14 et 15, est un 
sifflet ordinaire dans lequel un son synchronisant 
est appliqué transversalement a la base du jet d’air; 
ce son est guidé par une canalisation aplatie; la 
«buse auxiliaire» se trouve trés prés de la buse 
principale par laquelle s’échappe le jet d’air. 

Voici le principe de la synchronisation par buse 
auxiliaire: nous avons vu [6], Fig. 2, que dans un 
sifflet ordinaire, l’oscillation périodique du jet d’air 
est due a l’addition géométrique des vitesses du jet 
d’air et des «vitesses des particules» dans les ondes 
sonores émises par la cavité résonnant. En_inté- 
grant les vitesses transversales («vitesses des parti- 
cules dans l’onde sonore») le long du parcours du 
jet d’air entre la buse et le couteau, on obtient le 
déplacement transversal résultant au niveau du 
couteau. 

Comme il s’agit d’une fonction sinusoidale, son 
intégrale le long du parcours du jet, est également 
une fonction sinusoidale dont |’amplitude dépend 
de la distance entre la buse et le couteau et de la 
vitesse du jet d’air (cette derniére déterminant le 
temps de parcours, donc la durée de |’intégration 
des vitesses pour chaque particule d’air). Cette 
amplitude d’oscillation du jet au niveau du couteau 
doit s’annuler pour certaines valeurs de la vitesse 
du jet; ce fait a été confirmé expérimentalement. 

Si, au lieu d’appliquer une onde sonore tout le 
long du parcours du jet d’air entre la buse et le 
couteau, on applique une onde sonore seulement a 
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Vorigine du jet, prés de sa sortie de la buse, le 
phénomeéne est différent: tout le jet (qui est sinu- 
soidal sous l’effet du son émis par le sifflet, ou 
droit, en absence de ce son) s’incline et amplitude 
des déplacements aw niveau du couteau est beaucoup 
plus grande. Ainsi un son trés faible suffit pour 
synchroniser les oscillations d’un jet puissant, a con- 
dition d’étre appliqué trés prés de la buse émettant 
ce jet. Ce résultat est obtenu au moyen d’une 
deuxiéme buse (la buse auxiliaire) amenant le son 
de synchronisation dans le voisinage de l’origine du 
jet, perpendiculairement 4 ce jet. 

Le role de la buse auxiliaire est donc absolument 
analogue a celui de la grille dans une lampe triode: 
de méme que la grille placée dans le voisinage im- 
médiat de la cathode émettant un puissant flux 
d’électrons permet de commander ce flux par des 
tensions électriques trés faibles, — la buse auxi- 
liaire, placée dans le voisinage immédiat de la buse 
principale émettant un puissant jet d’air, permet 
de commander ce jet d’air par des ondes sonores 
trés faibles. Un sifflet 4 buse auxiliaire est donc une 
veritable «triode acoustique». 

L’expérience montre que la buse auxiliaire permet 
une synchronisation parfaite non seulement des sif- 
flets a basse pression mais également des sif- 
flets a haute pression et que le son émis par un sif- 
flet Levavasseur permet la synchronisation d’une 
vingtaine de sifflets Levavasseur de mémes dimen- 
sions munis de buses auxiliaires. I] est donc possible 
de réaliser des panneaux synchronisés de sifflets a 
basse ou a haute pression. Un élément de tel pan- 
neau est représenté sur la Fig. 16. La buse auxiliaire 
est annulaire, tout comme la buse principale. Le 
son de synchronisation est amené par un tube (un 
tube de 18mm de diamétre intérieur dans le cas 
de sifflets de 40 W fonctionnant a 10 kHz). 

Pour que tous les sifflets d’un panneau fonction- 
nent en synchronisme parfait, il faut évidemment 
que les tubes de synchronisation soient accordées a 
la fréquence du sifflet pilote et que leurs longueurs 
différent d’un nombre entier de longueurs d’ondes 
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Fig. 16. Sifflet Levavasseur a buse auxiliaire (élément 
d’un panneau sifflant). 
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de ce son. On voit le schéma d’un tel panneau de 
19 sifflets synchronisés par un 20™° sur la Fig. 17 
et l’aspect extérieur de ce panneau sur la Fig. 18. 


=— Air comprimé 


i es de synchronisation 
; = 


<— Air comprimé 


Panneau de 
sifflets synchronisés 
Coupe Face 
Fig.17. Panneau de 19  sifflets synchronisés par 


un vingtieme. 


Fig. 18. Vue d’un panneau de 19 sifflets synchronisés 
avec quatre éléments en place. Devant le pan- 
neau deux autres éléments synchronisés. 


La synchronisation a été contrélée expérimentale- 
ment en captant le son émis par un microphone con- 
necté a un oscilloscope cathodique: sur la Fig. 19 
on voit les battements (d’ailleurs parfaitement 
audibles) en absence de la synchronisation; sur la 
Fig. 20, les sifflets sont synchronisés, il n’y a plus 


de battements. un son unique est émis. 


Fig. 19. Battement entre le son émis émis par deux 
sifflets non synchronisés (vu sur |’écran d’un 
oscilloscope cathodique). 
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Fig. 20. Son unique émis par deux sifflets synchronisés. 


7. Mesure de la puissance des sifflets 


La puissance acoustique de la composante audible 
est mesurée au moyen d’un microphone étalonné 
placé a une assez grande distance du sifflet, dans 
une chambre sans échos. On tourne le sifflet ce qui 
permet d’obtenir un diagramme de rayonnement et, 
en intégrant les intensités dans différentes directions, 
on détermine la puissance acoustique totale aux 
fréquences audibles. 

La difficulté de cette mesure vient de la trés 
grande intensité du son a mesurer: la membrane 
d’un microphone placé trop prés du sifflet est trés 
rapidement brisée, et méme a une distance de plu- 
sieurs metres le son dépasse 145 dB. Or, les résul- 
tats obtenus en utilisant des microphones «protégés» 
(placés dans une boite) sont généralement faux, 
probablement a4 cause des résonances des parois de 
la boite, des vibrations parasites du support du 
microphone relié a cette boite protectrice, etc... . 
Il est préférable dutiliser pour ces mesures un 
microphone trés robuste, peu sensible, «nu». 

La mesure de la composante ultra-sonore se 
heurte a une autre difficulté. L’expérience montre 
que les microphones dits «ultra-sonores» étalonnés 
donnent des résultats peu sirs® et seul le procédé 
de la pression de radiation sonore [9], [10] convient 
pour ces mesures. Or, le son émis par des sifflets et 
par des sirénes s’accompagne toujours d’un jet d’air 
et ce jet, méme lorsqu’il est trés faible (sifflets a 
séparation du jet et du son) risque de fausser les 
mesures. Mais il est possible de séparer complete- 
ment les ultrasons et ce jet d’air au moyen d’un 
miroir (plaque métallique lourde) 4 45°, comme le 
montre la Fig. 21. En effet, la longueur d’onde de 
Pultra-son étant trés inférieure aux dimensions du 
«miroir», ils se réfléchissent 4 90°, tandis que le jet 
d’air suit la surface du miroir. On vérifie facilement 
en déréglant le sifflet (ou en arrétant la siréne) que 
le jet d’air seul n’influe nullement sur la balance. 
On utilise toujours un obstacle conique, pour éviter 


8 Des mesures faites «a quelques décibels prés» ne 
présentent aucun interét lorsqu’il s’agit d’évaluer la 
puissance acoustique en Watt. 
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Fig. 21. Mesure de l’intensité du son émis par un sifflet 
par le procédé de la pression de radiation 
sonore. 


tout risque d’ondes stationnaires et de résonance 
entre le générateur du son et |’obstacle. La formule 
a utiliser est 


(9) 
ou P est le poids en grammes mesuré par la ba- 
lance, S la surface de la base du céne en cm? et 


© Vangle que font les génératrices du céne avec 
sa base. 


8. Sirénes industrielles [8], [11] 


Contrairement a une opinion assez repandue, il ne 
suffit pas d’interrompre périodiquement un jet d’air 


pour le transformer intégralement en énergie sonore. 


Les troncons d’un tel jet, freinés par l’air ambiant 
se recollent et reconstituent un jet d’air continu 
silencieux. C’est pourquoi, le rendement d’une siréne 
est déplorable lorsqu’elle n’est pas munie d’un pa- 
villon suffisamment allongé. C’est dans le pavillon 
que prennent naissance, sous l’effet de l’inertie du 
jet brusquement interrompu, les dépressions périodi- 
ques engendrant les ondes sonores. 

Pour cette raison, les sirénes pour les alertes anti- 
aériennes a pavillons multiple allongés, sont trés ef- 
ficaces. Par contre les sirénes industrielles, surtout 
les sirénes dites «ultra-sonores» a pavillon trés évasé 
[12], [13], ont généralement un trés mauvais rende- 
ment et une faible puissance ®. 

Des sirénes bien meilleures, a pavillon allongé, 
ont été cependant réaliées aussi bien aux Etats-Unis 
{15] qu’en Allemagne [16] ces derniéres ont un 
rendement entre 25 et 30% a 5 kHz. 

En méme temps que des dépressions périodiques 
dans le pavillon, ’inertie du jet d’air brusquement 
interrompu, produit des compressions périodiques 
a l’intérieur méme de la siréne et ces compressions 
contribuent également a |’émission sonore, mais 
seulement a certaines fréquences: les sirénes pré- 


® Le calcul du rendement [14] basé sur l’hypothése 
inexacte d’une transformation totale, de 100% de 
lénergie du jet interrompu en énergie sonore, ne cor- 


_ respond a aucune réalité physique. 


V. GAVREAU: GENERATEURS PNEUMATIQUES D’ULTRA-SONS 


129 


sentent des fréquences de résonance trés prononcées, 
pour lesquelles le rendement est meilleur que pour 
les autres fréquences. 

Une difficulté particuliére aux sirénes est la 
nécessité d’assurer une étanchéité parfaite entre le 
rotor et le stator, sans provoquer un frottement ex- 
cessif. 


9. Avantages et inconvénients des sifflets 
et des sirénes 


Les sirénes permettent une émission de _ trés 
grande puissance acoustique, surtout aux basses fré- 
quences. Elles peuvent étre facilement synchronisées 
et utilisées en batteries. Mais elles exigent un aju- 
stage de haute précision et le fait qu’elles compor- 
tent des organes tournants et exigent un moteur 
d’entrainement, rend pratiquement impossible leur 
emploi a hautes températures et sous fortes pressions 
(bruleurs de charbon pulvérisé et de mazout, syn- 
thése de l’ammoniac). 

Les sifflets permettent d’obtenir trés facilement 
des ultrasons de trés haute fréquence (plus exacte- 
ment, des sons complexes a forte proportion de com- 
posantes ultra-sonores), convenant bien pour la coa- 
gulation de fumées et de poussiéres. On peut utiliser 
pour leur construction n’importe quel métal (acier 
inoxydable, bronze, etc... 
emploi pour différents traitements chimiques, méme 


.) ce qui permet leur 


a hautes temperatures et fortes pressions, en les ali- 
mentant au besoin directement par le mélange gazeux 
a traiter. 

Le principal inconvénient des installations a siré- 
nes et a sifflets est une consommation assez im- 
portante d’énergie, nécessaire pour assurer leur 
fonctionnement. Mais des progrés importants ont 
été réalisés réecemment et ont permis d’accroitre con- 
sidérablement leur rendement, donc de réduire la 
consommation de |’énergie. 


10. Applications pratiques 


Les premiéres applications envisagées ont été la 
précipitation des fumées, et poussiéres industrielles 
[17], [18] et la précipitation des brouillards. Mais 
on est rapidement arrivé a la conclusion que seuls 
certains types de poussiéres et de fumées pouvaient 
étre précipitées par ce procédé: les poussiéres humi- 
des ou grasses, collantes, et les fumées goudronneu- 
ses (telles que la fumée de tabac, la fumée de caout- 
chouc, la fumée des usines fabriquant l’oxyde de 
zinc, etc....). Dans d’autres cas, on est obligé d’uti- 
liser une pulvérisation d’eau pour permettre le col- 
lage et la précipitation de la poussiére par les ultra- 
sons. 

Une application assez intéressante parait étre la 
précipitation des goudrons «non-coagulables» par 
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tous les autres procédés, notamment des goudrons 
trés finement pulvérisés contenus dans le gaz des 


gazogenes. Dans ce cas, le sifflet est alimenté en gaz ' 


comprimé purifié, pour éviter la formation d’un 
mélange explosif avec lair. Nous avons également 
utilisé dans ce but un sifflet alimenté en vapeur sous 
pression fournie par une chaudiére ordinaire. 

La précipitation des brouillards naturels a été 
réalisée au moyen de sirénes, mais le son audible 
de 3 a 5 kHz émis par de telles sirénes est absolument 
insupportable et on a été obligé d’abandonner ce 
procédé. I] sera peut-étre possible de précipiter le 
brouillard par des batteries de sifflets, a fréquence 
plus élevée, beaucoup moins pénible et d’ailleurs 
fortement absorbée par le brouillard en déhors de la 
zone «nettoyée»: La fréquence de 10 kHz parait étre 
particuliérement efficace. (Plus exactement, il s’agit 
d’un mélange de 10 kHz et de 150 kHz émis par un 
sifflet Levavasseur; les essais ont été faits sur un 
brouillard d’eau produit artificiellement, au labora- 
toire.) 

Mais les applications les plus intéressantes des sif- 
flets (les sirénes ne pouvant pas étre utilisées dans 
ce but) sont la combustion accélérée et la catalyse 
accélérée. Les essais préliminaires ont donné des 
résultats absolument remarquables et des sifflets- 
bruleurs industriels ont été construits. 


11. Conclusion 


Une trés longue et difficile mise au point des 
émetteurs pneumatiques d’ultra-sons intenses (siré- 
nes et sifflets) a permis d’accroitre considérablement 
leur rendement et les puissances acoustiques émises. 
Leur emploi industriel a ainsi été rendu possible et 
un certain nombre d’installations sont en cours de 
réalisation dans différents pays (en France, en Po- 
logne, en Hongrie). 

Il s’agit de dépoussiérage et de précipitation de 
fumeées, pour lesquels on peut utiliser indifféremment 
des sirénes ou des sifflets, mais il faut que la pous- 
siére ou la fumée ne soit pas «séche», car seule la 
poussiere collante et la fumée goudronneuse peuvent 
ainsi étre précipiteées. 

Les applications de loin les plus importantes et 
les plus intéressantes concernent l’accélération de 
certaines réactions chimiques: combustion, catalyse. 
Pour ces application seuls les sifflets conviennent, 
les sirénes tournantes ne pouvant pas étre utilisées 
aux temperatures élevées et sous trés grandes pres- 
sions. 
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Une autre application, -peut-étre possible, est la 
précipitation des brouillards. La encore, seuls les 
sifflets peuvent étre utilisés, car les sirénes de grande 
puissance émettant des sons de fréquence peu €élevee, 
sont insupportables pour le personnel travaillant 
dans les ports et sur les aérodromes. 
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EULSE TECHNIQUES 
MEASURING. ULTRASONIC ABSORPTION IN LIQUIDS 
by J. H. Anpreaz®, R. Bass*, E. L. Heaserx+ and J. Lams* 
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Summary 


The paper is concerned with the construction and operation of equipment for the 
measurement of the ultrasonic absorption in liquids by the pulse method over a range of 
frequencies from 1 to 200Mc/s. General considerations of the design of experimental 
systems are discussed together with the practical limitations. 

Details are given of the electrical and mechanical parts of three systems which are in 
current use, and permit the measurement of ultrasonic absorption within an accuracy of 
+2%. The first two systems, A and B, are designed to operate up to a frequency of 
200 Mc/s and illustrate alternative lines of approach to the problem. The third system C 
can be used for measurements on liquids under pressures up to 100 atm. It incorporates a 
common transmitter-detector system and provides a longer ultrasonic path than do the 
first two systems: thus smaller absorption coefficients can be measured or measurements 
made at correspondingly lower frequencies, down to a minimum of 1 Mc/s. 

Attenuation meastrements are referred to various piston attenuators which have been 
constructed and checked against each other. 


Sommaire 


On traite de la réalisation et de l’utilisation d’un appareillage destiné 4 mesurer l’absorp- 
tion des ultra-sons dans les liquides, par la méthode des impulsions, entre 1 et 200 MHz. 
On en discute les lignes générales, ainsi que les limites pratiques. 

On décrit en détail les parties électriques et mécaniques de trois systemes couramment 
utilisés, qui permettent de mesurer l’absorption des ultra-sons 4 +2% prés. Les deux pre- 
miers dénotés A et B, sont destinés 4 fonctionner jusqu’a 200 MHz et représentent deux 
facons différentes d’aborder le probleme. Le troisiéme systéme dénoté C peut étre utilisé 
pour des mesures dans des liquides soumis a des pressions allant jusqu’a 100 atm; il 
comporte un systeéme transmetteur-détecteur commun et un trajet d’ultra-sons plus long 
que dans les deux premiers; ceci permet de mesurer des coefficients d’absorption plus 
faibles, ou. encore de descendre a des fréquences plus basses (jusqu’a 1 MHz). 

Les mesures d’affaiblissement sont données par comparaison avec différents affaiblisseurs 
a piston, qui ont été construits a cet effet, et que l’on a comparés entre eux. 


Zusammenfassung 


Die Arbeit befaBt sich mit Aufbau und Wirkungsweise einer Apparatur zur Messung 
der Ultraschallabsorption in Fliissigkeiten nach der Impulsmethode fiir Frequenzen zwi- 
schen 1 und 200 MHz. Es werden allgemeine Gesichtspunkte zur Konstruktion der Appa- 
ratur und die praktisch.'gegebenen Grenzen diskutiert. Einzelheiten des elektrischen und 
mechanischen Aufbaues dreier Systeme werden angegeben, die Absorptionsmessungen mit 
einer Genauigkeit von 2% gestatten. Die ersten beiden Systeme arbeiten bis zu Frequenzen 
yon ca, 200 MHz und benutzen verschiedene Methoden. Das dritte System kann fiir Drucke 
bis 100 atm benutzt werden. Es enthialt ein gemeinsames Sende- und Empfangssystem und 
benutzt eine lingere Ultraschall-Wegstrecke als die ersten beiden Systeme. Somit kénnen 
kleine Absorptionskoeffizienten gemessen oder Frequenzen bis hinunter zu 1 MHz benutzt 
werden. Die Diimpfungswerte wurden auf verschiedene geeichte Dampfungsglieder bezogen. 


1. General considerations oscillator and is fed to a piezoelectric transducer. 


The resulting ultrasonic pulse from the transducer 


1.1. Features of the different possible experimental 
arrangements 


In principle, a short radio-frequency pulse con- 
taining an appreciable number of oscillations at the 
frequency f, is generated by a suitably modulated 
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passes through the liquid to be measured and is 
reconverted into an electrical pulse by a detecting 
transducer. After suitable amplification and de- 
modulation the envelope of the received pulse is dis- 
played on an oscilloscope. By varying the distance 
between the launching and detecting transducers and 
by measuring the change in amplitude of the detected 
pulse, a plot of attenuation versus path length is 
obtained which gives the absorption coefficient a 
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directly. From this and the known frequency, /, the 
quantity a/f? may be calculated. 

Measurements of the attenuation in the liquid are 
referred to a precision attenuator. The attenuator may 
be placed (i) in series with the detecting transducer, 
or (ii) in series with the oscillator. In both these 
cases variations in the attenuation in the ultrasonic 
path are compensated for by equivalent changes in 
the setting of the attenuator, maintaining a constant 
signal level on the oscilloscope. A disadvantage of 
these arrangements is that the accuracy depends on 
the gain of the electrical receiver remaining constant 
during the period of a measurement (10 —15 minu- 
tes). The amplitude of the received signal before it 
reaches the receiver input is reduced by losses in the 
acoustic system and by the insertion loss and added 
attenuation of the attenuator. Thus, for a given 
oscillator output and receiver noise level, the signal 
to noise ratio is determined by the sum of these 
attenuations and remains constant throughout a 
measurement. Considerable improvement may be 
obtained with arrangement (i) if a preamplifier can 
be placed between the receiving transducer and the 
attenuator. The signal to noise ratio is then de- 
termined by the attenuation in the acoustic system 
and the noise level in the preamplifier only, provided 
that the gain of the latter is large compared with the 
insertion loss of the attenuator. It is, however, es- 
sential that the pre-amplifier should have a linear 
gain characteristic over the required dynamic range 
of input signals and this is not easily combined with 
a high gain. On the other hand arrangement (ii) has 
obvious advantages over (i) if the permissible output 
from the oscillator is limited by “‘spark-over” at the 
transducer. 

A third arrangement which is particularly satis- 
factory at high frequencies requires the addition of 
an electronic delay circuit and a second (com- 
parison) oscillator tuned to the same frequency as 
the main oscillator. The pulse from the comparison 
oscillator is delayed with respect to that from the 
main oscillator so that the envelopes of the two 
pulses are displayed adjacently when they reach the 
oscilloscope. Both pulses pass through the same 
amplification path. The receiver is operated at maxi- 
mum gain only for the highest attenuator settings. 
the gain being reduced successively as the atten- 
uation decreases, with a consequent improvement 
in the signal/noise ratio. There is no need for the 
receiver to have a constant gain over the period of 
time required for a measurement. 

The detecting transducer can be dispensed with by 
using a reflector, in which case the ultrasonic wave 
is returned to the transmitting transducer for de- 
tection. In this common transmitter-detector system 
the receiver must be capable of recovering from the 
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application of the oscillator output. The double path 
length provided by this arrangement can be advant- 
ageous at low frequencies when the absorption is 


‘‘Jow. Also at high frequencies the problem of produc- 


ing transducers with identical harmonic resonant 
frequencies is avoided. 


1.2. Practical limitations to the design of experi- 
mental systems 


A variety of practical limitations must be con- 
sidered in the design of experimental systems utilis- 
ing the pulse technique: 

(i) At any one frequency there is a wide variation 
in the absorption coefficient of different liquids. For 
example, in water (which has a low absorption coef- 
ficient), /f?=25°10-* s* cm-*,, and Vateag re 
quency of 50 Mc/s a liquid path of 5.5cm would 
be required to reduce the amplitude of an ultra- 
sonic wave by 30 dB; in carbon tetrachloride (a/f? 
=550:10~'7s?cm~*) the absorption is high and 
the corresponding path length would be 0.25 cm. 
This means that separate mechanical systems may be 
needed for measurements on different liquids at the 
same frequency. 

(ii) For a specific liquid, in the absence of re- 
laxation processes, the absorption coefficient is pro- 
portional to the square of the frequency: thus, to 
reduce the amplitude by, say, 30 dB, the path length 
required at 10 Mc/s would be one hundred times 
that at 100 Mc/s. In the case of carbon tetrachloride 
cited above the path length for 30 dB attenuation 
would be 6.25 cm at 10 Mc/s and only 0.0625 cm at 
100 Mc/s. These two conditions could not easily be 
dealt with in the same mechanical system. 

(iii) At low frequencies diffraction effects due to 
the finite radius, a, of the transmitting transducer 
may have to be taken into account. The sound field 
can be divided approximately into the Fresnel region 
near the transducer and the more distant Fraunhofer 
region. The boundary between these two lies at a 
distance of about a?/A from the transmitting trans- 
ducer. In the Fresnel region, the total ultrasonic 
energy radiated from the transducer is still confined 
to a narrow beam, but the pressure amplitude 
fluctuates rapidly both across the beam and along 
its axis. If, however, the receiving transducer is made 
somewhat larger than the transmitting transducer 
and so covers the whole cross section of the beam, 
constant signal is obtained throughout the Fresnel 
region in the absence of absorption. This is not 
exactly true for transmitting and receiving trans- 
ducers of equal radius, but even then diffraction ef- 
fects are in practice usually negligible at frequencies 
above 5 Mc/s and with ultrasonic pathlengths not 
exceeding a?/22[1]. On the other hand, in the 
Fraunhofer, or diverging, region of the sound field 
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the aperture of the detecting transducer may be re- 
stricted so as to record as nearly as possible the 
pressure on the axis; it is then necessary to apply a 
correction to the measured attenuation to allow for 
the divergence of the sound beam in order to obtain 
the true ultrasonic absorption. A discussion of these 
corrections has been given by Pinkerton [1]. 

(iv) When the time taken for the ultrasonic pulse 
to traverse the liquid path is less than the pulse 
duration an additional time delay must be provided 
in the ultrasonic system. This serves to separate on 
a time scale the required signal from the transmitted 
pulse which, if two transducers are used, by-passes 
the ultrasonic system by radiation, or if a common- 
transmitter-detector system is employed, is fed 
directly to the receiver. Owing to the high power 
of the oscillator and the sensitivity of the receiver. 
it is not generally feasible to eliminate this ‘pick-up’ 
pulse. The necessary time delay may be provided 
conveniently by introducing a length of fused quartz 
rod into the ultrasonic path. Fused quartz possesses 
a low attenuation coefficient and gives a delay of 
about 1.7 us/em. In practice further advantages ac- 
crue from the use of such delay rods and these will 
be discussed in chapter 2. 

(v) The temperature of the test liquid must be 
kept constant to within +0.1° C in order to prevent 
local thermal fluctuations from causing variations in 
the received signal. Heating of the liquid by the dis- 
sipation of ultrasonic energy does not present a pro- 
blem because the duty cycle (on/off ratio) may be 
kept as small as necessary, say, 107. 


1.3. Standards of attenuation 


The overall accuracy of the absorption measure- 
ments depends upon the accuracy of the attenuator 
employed. Piston attenuators have been used 
throughout this work [2], [3] and in order to 
verify the theoretical attenuation constants various 
attenuators have been constructed and compared 
with each other. These include types operating in the 
E,, (cireular), H,, (circular) and Ho, (rectangular) 
modes. Two attenuators of different transverse di- 
mensions have been built to operate in each of these 
modes. The design of the H,, attenuators followed 
the description by Garnssorovcu [2]. Diagrams il- 
lustrating the construction of the Ey, and Ho, at- 
tenuators are shown in Fig. 1. In all cases the cal- 
culated attenuation constants were found to agree 
to within better than +0.5%, which was the limit 
of accuracy of the method of comparison. 


1.4. Transducers 

X-cut quartz crystals are normally used as trans- 
ducer elements at all frequencies from 1 to 200 Mc/s, 
but polarised barium titanate discs are suitable 
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below about 10 Mc/s. The electrical impedance pre- 
sented to an oscillator by a quartz transducer be- 
comes excessive at the lower frequencies and power 
is more readily transferred from the oscillator to the 
lower impedance offered by a barium titanate trans- 
ducer. However, this advantage of the barium titan- 
ate transducer is not particularly significant since 
there is no need for large amounts of ultrasonic 
power at the lower frequencies. In fact, there is little 
difficulty in obtaining a wide range of attenuation 
(e. g. 50 dB) for measurements taken with apparatus 
operating below 100 Mc/s. 
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electrode 


electrode 


IS 
(100 divs. =1dB) 


(0-130 dB) 


Fig. la. Eo, (circular) piston attenuator. 
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Fig. lb. Hy, (rectangular) piston attenuator. 


The thickness-vibration of an X-cut quartz crystal 
disc is resonant when the thickness of the disc is 
equal to one half wavelength or an odd number of 
integral half wavelengths. One half wavelength is 
about 0.014in. at 20 Mc/s in quartz so this fre- 
quency is about the highest at which a crystal can be 
operated conveniently in a fundamental frequency, 
or half wavelength, mode. Higher frequencies are 
reached by using the odd harmonic resonances of 
the crystals and, by selecting a crystal with a con- 
venient fundamental frequency, a complete range of 
discrete frequencies can be covered. The bandwidth 
of a crystal increases with the frequency of its fun- 
damental vibration but is the same at all harmonic 
frequencies for the same fundamental: therefore it 
is advisable to select crystals with the highest 
possible fundamental frequency if a reasonably 
broad frequency response is required at harmonic 
frequencies above 100 Mc/s. This is important in 
two-crystal systems which depend on the optimum 
behaviour of two crystals at the same frequency. 
Fundamental frequencies of 10 and 20 Mc/s have 
proved to be suitable for operation up to 200 Mc/s. 
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In order to apply an alternating voltage across a 
quartz crystal, an electrode is provided on each side 
of the disc and the area of at least one of the elec- 
trodes is limited to the area of the required ultra- ‘ 
sonic beam. The electrodes may be metallic layers 
evaporated on to the surfaces of the crystal, or 
they may consist of conducting surfaces between 
which the crystal is clamped. When fused quartz 
delay rods are employed, one electrode is con- 
veniently provided by evaporating a layer of 
chromium or gold on to the optically flat end of 
the rod before pressing the crystal disc against it. 
Also, to permit ultrasonic propagation from the 
crystal to the rod, it is necessary to indroduce a drop 
of non-volatile liquid (e.g. vacuum oil) between the 
crystal and the rod. Methods of mounting quartz 
crystals are illustrated in Figs. 2, 7 and 8. 

Experience with the ultrasonic systems A and B 
of section 2 has shown that it becomes increasingly 
difficult above 100 Mc/s to obtain the desired at- 
tenuation range of 30 dB for measurements. Neither 
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2. Detailed description of experimental systems 


Three of the experimental systems used by the 
authors have been selected to illustrate the salient 
features of apparatus employing the pulse technique. 
The first two systems, A and B, were designed for 
the purpose of measuring ultrasonic absorption in 
small quantities of liquids up to 200 Mc/s; both are 
included so that alternative lines of appreach to the 
design of equipment can be emphasized. The third 
system, C, is used for measuring the absorption in 
liquefied gases at pressures up to 100 atm. It incor- 
porates a common transmitter-detector system and 
provides a longer usable ultrasonic path than do the 
first two systems. These three sets of apparatus are 
described in turn. 


2.1. System A 
2.1.1. Mechanical design 


Two considerations arise in the design of the 
mechanical system, namely, the production of small 


Depth gauge head 
/ 


Side elevation 
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Fig. 2. Mechanical design of system A. 


of the systems allows a useful variation in the ampli- 
tude of the received signal for more than 25 dB at 
frequencies above 200 Mc/s. The reason for this 
decrease in the available attenuation range is not 
clear at present, but it is likely to be connected with 
the operation of the crystals, the efficiency of the 
coupling between the crystals and the delay rods 
and the attenuation in the latter [4], [5], [6]. 


displacements and the alignment of the two trans- 
ducers. It is desirable to obtain ten, or more, points 
on the final graph of attenuation plotted against 
distance and the points should be accurate to ©1%, 
so it is necessary to produce measurable displace- 
ments of, say, 0.2% of the total path length covered 
in the measurement. Also, to obtain the maximum 
ultrasonic signal through the system, parallelism of 
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the transducers must be achieved and maintained to 
within a fraction of an acoustic wavelength during 
displacement. These problems can be illustrated by 
the case of carbon tetrachloride at a frequency of 
200 Mc/s. For a measuring range of 30 dB, the total 
change in path length is only 0.016cm and so 
measurable displacements of the order of a micron 
must be possible. Also the ultrasonic wavelength is 
only 6 microns so parallelism should be maintained 
to within about one wavelength of visible light. At 
lower frequencies the problem is proportionately 
simpler. 

The mechanical design of system A is illustrated 
in Fig. 2. The detecting transducer and its associated 
fused quartz delay rod are supported in a rigid steel 
framework which slides upon two steel pillars, the 
surfaces of which are hardened and ground to be 
accurately straight and cylindrical. The framework 
is kinematically mounted, being supported by five 
hardened steel balls, which are held in firm contact 
with the bearing pillars by two flat phosphor-bronze 
springs: as can be seen from Fig. 2, there are four 
point contacts on one pillar and one oneihe other. 
The whole moving system is counterbalanced by a 
brass weight and is driven by a micrometer depth 
gauge. The transmitting crystal, quartz rod and 
matching lines are rigidly fixed to the base of the 
apparatus and alignment is achieved by tilting the 
two bearing pillars with respect to the base. These 
pillars are fixed into accurately matched holes in the 
two end plates, the lower plate being spring-loaded 
against the base of the apparatus and tilted by means 
of three fine-pitched screws (40t.p.i.). The dis- 
placement of the upper fused quartz rod is measured 
by a dial gauge, the sensitivity of which is chosen to 

_ suit the absorption of the test liquid. 

The design of the vessel containing the liquid is 
influenced by the need to isolate the temperature- 
controlling liquid from the test sample. The con- 
trolling liquid flows between the inner and outer 
walls of a double wall container (Fig.3), which 
fits on to the lower fused quartz rod by a tapered 
ground joint. This seal is found to be liquid tight 
without the use of any grease, which might con- 
taminate the test sample and has proved to be 
satisfactory over the temperature range —60° to 
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Temperature 
controlling 


liquid 


Tapered ground joint 
Fig. 3. Liquid vessel for system A. 
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+ 100° C. The temperature of the liquid is measured 
by a thermocouple attached to the upper fused 
quartz delay rod. The diameter of the fused quartz 
rods is 9/16” and, when these are touching, less 
than 10cm? of test liquid are needed to fill the 
vessel. 


2.1.2. Electrical design 


In order to achieve simplicity of design and con- 
struction, self-excited oscillators are employed. By 
using the push-pull oscillator circuits of Fig. 4 only 
a single tuned-circuit is required. Normal induct- 
ances and condensers are satisfactory for the tuned- 
circuit up to about 120 Mc/s and above this fre- 
quency Lecher lines are used. The output coupling 
consists of a coil or loop placed close to the tuned 
circuit. The valve chosen for these oscillators is the 
double tetrode type QQV 06/40. The anode voltage 
for normal operation should not exceed 600 V, but, 
since modulation is applied to the anodes in pulses 
of 5s duration and at a repetition rate of 500 per 
second, voltages up to 5kV have been employed. 

The circuit of the high power modulator built to 
generate the 5 kV pulses in shown in Fig. 5. A delay 
cable (Telcon Type Z3M) of characteristic im- 


@) O0V 06/40 


Fig. 4. (a) Tuned-line oscillator, for use above 120 Mc/s 
(LL Twin Lecher line, S Movable short cir- 
cuit) ; 

(b) Lumped circuit oscillator, for use up to 


120 Me/s. 
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Fig. 5. Modulator circuit system ‘A’. 
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pedance equal to 330 ohms is used for the storage 
element and this is discharged in series with the load 
resistance through an hydrogen thyratron type 


BT 83. The impedance of the combined load formed ' 


by the load resistance in parallel with the oscillator 
must be matched to the characteristic impedance of 
the delay line for correct operation [7]. In order to 
overcome the effect of fluctuations in the impedance 
presented by the various oscillators it is desirable 
that the characteristic impedance of the delay line 
should be small compared with the impedance of the 
oscillator. It may be noted that the use of a pulse 
transformer which would require more careful im- 
pedance matching is avoided in this circuit. 
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Fig. 6. Schematic diagram of system ‘A’. 


A block diagram of system A is given in Fig. 6. 
The pulse generator, which triggers the modulator 
and the comparison oscillator, is of conventional 
design and this is also true of the delay circuit. The 
matching indicated in Fig. 6 is achieved by the use 
of coaxial lines of variable length. The transducers 
are X-cut quartz crystals having a fundamental fre- 
quency of 20 Mc/s and a capacity of 45 pF. In the 
case of the transmitting transducer an additional 
‘stub’ line is provided to tune out this capacity. At 
frequencies below about 100 Mc/s matching to the 
attenuator is more easily accomplished by the use of 
tapped tuning coils across the input and output 
capacities. . 

Modified Radar receivers have been employed to 
provide the necessary amplification for the detected 
pulse. The GEE receiver type R 1355 has been used 
up to 130 Mc/s and the receiver R 1359 for the 
higher frequencies. In order to obtain a better signal- 
to-noise ratio in the latter, the initial stage of inter- 
mediate frequency amplification has been replaced 
by a Cascode low-noise amplifier [8]. 
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Use of a comparison pulse for attenuation mea- 
surement implies that both comparison and _ultra- 
sonic signals shall follow identical paths through the 
amplifier. To achieve this a small coupling coil or 
probe is introduced near the input radio-frequency 
tuning circuit of the receivers for injection of the 
comparison signal. This circuit is very loosely 
coupled and little interaction occurs between the two 
signal paths. 

To reduce the effects of receiver paralysis result- 
ing from overloading by the transmitter ‘pick-up’ 
pulse, ‘blanking’ techniques may be employed. While 
the oscillator is generating power, a pulse is applied 
to the grids or cathodes of the initial amplifying 
stages of the receiver to reduce their gain to zero. 
As a further precaution a clipping network may be 
inserted after the second detector to prevent pulses 
of more than a predetermined height from entering 
the video-frequency amplifier. 


2.1.3. Measuring procedure 


The piston attenuator is set to minimum attenuation 
and, using its known calibration, the comparison 
oscillator is tuned to the required frequency (an odd 
harmonic of the fundamental frequency of the crystal 
transducer). No matching is attempted at this stage 
but sufficient signal reaches the receiver to permit 
it to be tuned. Next an absorption wavemeter is 
tuned to the output of the signal generator and this 
is now used to check the tuning of the transmitting 
oscillator. 

The transmitter and receiver are connected to the 
ultrasonic system and the tuning of the former is 
checked. A small quantity of a liquid of low ab- 
sorption, for example, alcohol, is put in the vessel 
and the two fused quartz rods are brought into 
contact. Under these conditions it is usually pos- 
sible to observe a small signal which has passed 
through the ultrasonic system. The electrical system 
is now tuned and the matching lines adjusted to 
provide maximum signal. The final frequency does 
not necessarily correspond to an exact harmonic of 
the fundamental frequency of the transducer. The 
fused quartz rods are then separated to increase the 
ultrasonic, path length and to reduce the signal 
amplitude to a little above noise level and the 
mechanical alignment is carried out as follows. The 
levelling screws are adjusted in turn to maximise the 
signal passing through the system: when an apparent 
maximum is obtained the apparatus is tilted grossly 
out of alignment in two directions at right angles to 
ensure that this is the largest maximum. The nuts 
controlling the compression springs are then tigh- 
tened to prevent accidental alteration of the align- 
ment and, finally, the tuning of the electrical system 


is rechecked. 
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It now remains to adjust the comparison oscil- 
lator. This is set to its pulse modulation position and 
the comparison pulse is made to overlap the first 
received ultrasonic pulse on the oscilloscope. Under 
these conditions beating between the two pulse 
envelopes can be observed and this is reduced to a 
minimum by tuning the comparison oscillator. Al- 
though there is no phase coherence between the 
transmitting and comparison oscillators, their sta- 
bility permits the reduction of the beat to a very 
low rate, implying operation at identical frequencies. 
Now the attenuation in the comparison path is 
gradually increased and the attenuator matching is 
adjusted simultaneously to an optimum position. 

The frequency of the oscillations is measured by 
means of a heterodyne wavemeter [9], the reference 
erystal of which is periodically checked against the 
Droitwich transmitter as a standard. Generally the 
frequency of the comparison oscillator is measured 
directly: this produces a narrower frequency spec- 
trum having a longer pulse length than the main 
oscillator or it can be switched to continuous wave 
operation. 

It is found that mechanical alignment at the higher 
frequencies above 100 Mc/s is complicated by the 
change in path length produced during the angular 
adjustment. For this reason the alignment procedure 
is performed at about 100 Mc/s and only residual 
corrections are made at higher frequencies. After 
careful alignment at any one frequency, note is made 
of all electrical settings so that the conditions can 
be repeated readily at a later date. 

In conclusion it may be remarked that the system 
has been used successfully as a continuous wave 
interferometer for velocity measurements. The com- 
parison oscillator is then employed to drive the 
ultrasonic system, internal sine-wave modulation 
being used to render the received signal suitable for 
display on the oscilloscope. 


2.2. System B 
2.2.1. Mechanical system 


The precise mechanical movement (Fig.7) is 
obtained by the sliding of one Invar steel cylinder A 
inside a second Invar cylinder B. One of the two 
fused quartz rods, C, , is held in a conical hole along 
the axis of the cylinder, A, and the top end of this 
rod is also ground to take a small glass container. 
D, for the liquid to he measured. The other quartz 
rod C, is fixed into the Invar block E, which can be 
tilted about the centre point of the lower face of this 
rod; this is achieved by giving part of the block E 
a spherical surface F, which is held against the inside 
conical surface of the upper part of the cylinder B 
by four screws G. The adjustment of these screws 
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Fig. 7. Mechanical design of system ‘B’. 


tilts the axis of the fused quartz rod C, without 
varying the distance of it from the rod C, . 

The displacement of the cylinder A is transmitted 
to the dial gauge H through a vertical rod J, the 
two sections of which are thermally insulated from 
each other by a Bakelite connector K. For very short 
ultrasonic paths in the liquid between the fused 
quartz rods an accurate Swiss dial gauge is used 
having graduations down to 0.001 mm. For longer 
ultrasonic paths up to a maximum of about 1 cm a 
less accurate gauge is employed. The displacement 
of the cylinder A is controlled by the screw L which 
supports the cylinder through two thin steel wires, 
M. These wires act as thermal insulators so that the 
driving screw L can project outside the thermal 
insulation. 

The generating and detecting X-cut quartz crystals 
of 10 Mc/s fundamental frequency are pressed 
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against the outer ends of the fused quartz rods by 
means of the crystal holders N, and Nj. When the 


crystal holders and crystals are removed and the 


tilt screws G are loosened, it is found that optical 
interference fringes can be observed through the 
rods as the lower quartz rod is raised to touch the 
upper one. However, this method of alignment is 
insufficiently precise for operation above 100 Mc/s 
when it is necessary to adjust the tilt screws for 
maximum ultrasonic signal. 

About 12 cc of the liquid to be measured is intro- 
duced into the container D. A heavy copper wire coil 
is wound around the outside of the container and a 
periodically interrupted current of about 15 amperes 
is fed to the coil by a thermostat. This on-off cycle 
of current provides heating for proportional tem- 
perature control and also generates an intermittent 
magnetic field to drive a small annular magnetic 
stirrer Q. 

Three tapered holes P are ground to take standard 
B7 glass cones and these are used, respectively, for 
the insertion of two thermistors (Standard Tele- 
phones and Cables, type F) — one for measuring 
the temperature and the other for controlling the 
thermostat bridge — and for the introduction and 
removal of liquid from the container D. Since some 
liquids cannot be exposed to the atmosphere, the 
apparatus can be fitted with chromium-plated brass 
bellows which seal off the space around the tank 
between the upper and lower steel cylinders, A and 
B. When the bellows are in position the apparatus 
can be connected directly to a distillation plant via 
a B7 cone in one of the holes P, and air can be 
completely excluded from the space above the tank. 
Although the liquid is in direct contact with glass 
and fused quartz only, the vapour of acids will 
attack the Invar steel between the bellows and the 
quartz rods and it is necessary to dismantle the 
apparatus afterwards to clean off the slight cor- 
rosion. This involves the loss of the alignment bet- 
ween the fused quartz rods so, if corrosive liquids 
are to be used regularly, the steel should be coated 
or plated for protection. 

The apparatus was designed to operate over the 
temperature range from +100° to —150°C and 
low temperatures are attained by means of liquid 
nitrogen. The outer Invar cylinder B is surrounded 
by a copper jacket R, which is partially filled with 
a low melting-point liquid, and liquid nitrogen is 
forced into it by a cooling system similar to that 
described by Vaspa and Harr [10]. Around the 
copper jacket is a Dewar vessel, S, having a hole in 
the bottom to take the coaxial cable from the detect- 
ing crystal and the heater wire from the liquid 
container. The hole in the bottom of the Dewar 
vessel is fitted with glass bellows to allow for thermal 
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contraction of the inner wall. The thermal insulation 
of the apparatus is completed by two Paxolin discs, 
T, a Paxolin tube, U, which is hinged for removal, 
and a rubber tube which is rolled down over U. 
Before lowering the temperature of the system it is 
essential to dry the air inside and, in particular, to 
insure that the graphite lubricant between the slid- 
ing cylinders is quite dry. A separate thermostat 
regulates the temperature of the cooling liquids in 
the copper jacket. For temperatures above room 
temperature the jacket R and the Dewar vessel S are 
unnecessary and an additional heater coil is wrapped 
around the outer cylinder B to raise the mean tem- 
perature of the apparatus to a few degrees below 
that required for the liquid in the container D. 


7.22 Eheetrrcak system 


The electrical arrangement is similar to that of 
system A apart from a few relatively unimportant 
details. 


2.3. System C 


2.3.1. General considerations 


A hermetically sealed system has been constructed 
for measurements in liquids and liquefied gases at 
pressures up to 100 atm. over the frequency range 
1 to 100M/s. As discussed in Section 1.2. com- 
paratively large ultrasonic path lengths are required 
for absorption measurements in the region of 1 Mc/s 
while at the same time a simple and compact design 
is required for pressure applications. By employing 
a common transmitting and receiving transducer the 
effective ultrasonic path is doubled and the mechani- 
cal layout is considerably simplified. In addition 
this avoids the difficulty of producing electrical 
connections of constant impedance to a moving 
transducer. 

A disadvantage of the common T-R ‘system, 
however, is the unavoidable paralysis of the electri- 
cal receiver for a short time following the trans- 
mitter pulse. Further, the transducer itself will con- 
tinue to vibrate after the termination of the trans- 
mitter pulse and will therefore be unsuitable for the 
detection of weak echoes during a certain recovery 
period. 

An ordinary plane reflector is unsuitable for the 
measurement of relatively high absorption coefficients 
since if the reflector is set at sufficiently short range 
the reflected signal arrives at the transducer before 
the system has recovered from the transmitter pulse. 
This difficulty is overcome by the use of a fused 
quartz delay line as a reflector [11]. On the other 
hand, with low absorption coefficients, the delay 
line is unsatisfactory because pulses reflected several 
times within the liquid remain strong enough to 
interfere with the delayed signal; further the length 
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of the delay line restricts the available liquid path. 
The two reflectors are therefore used in turn. 

Pulse systems working as common T-R arange- 
ments at 7.5 Mc/s and above have been described 
by Prevtam and Gatr [12] and by Pinkerton [1]. In 
the present discussion attention is particularly focus- 
sed on the problems associated with the application 
of pressure and the measurement of low absorption 
coefficients in liquids at relatively low frequencies. 


2.3.2, The mechanical system 


A cross-section of the apparatus is shown in Fig. 8. 
The thick-walled cylindrical tube A, in the following 
_ referred to as the cylinder, has an inside diameter 
of 2”, an outside diameter of 5” and is 8” long. The 
top-plate B and baseplate C are both 3/4” thick and 
held to the cylinder by means of 5/8” diameter high 
tensile steel bolts. All metal parts, except for the 
bolts, are made of FMB grade stainless steel. The 
bore of the cylinder and the surface of the piston D 
are ground to secure smooth and strictly parallel 
movement of the latter. The reflector E has an ac- 
curately machined and polished surface and is 
screwed into the piston D. The piston is moved by 
turning the screwed rod F. The thread has 26 turns 
per inch and the position of the reflector can be read 
from an externally mounted revolution counter, 
which counts 100 units per turn of the drive. 
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Fig. 8. Mechanical design of system ‘C’ (for description 
see text). 


The drive is supported by the ball G and the rod 
H, the height of the latter being adjusted by means 
of the screw J in order to maintain a slight pressure 
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onthe ball bearing K and thus fix the position of 
the drive. The crystal holder (with the transducer 
M) is supported in such a way that it can be tilted 
about the centre of the top surface of the transducer. 
The tilting is carried out by means of two screws 
operated from outside which do not appear in the 
cross-section shown. By the eccentric arrangement 
of both ultrasonic path and drive with respect to the 
cylinder axis the maximum possible ultrasonic 
pathlength — in this case approximately 23 cm — 
is obtained for a given length of the cylinder. For 
the measurements of high absorptions the reflector E 
in Fig.8 is replaced by a fused quartz delay line. 
The quartz rod has a diameter of 17.5 mm and a 
length of 128 mm providing a time delay of 43 us 
per double journey. The end faces are ground parallel 
to each other to within 1’ and perpendicular to the 
axis to within 2’. 

The cylinder is surrounded by a liquid bath which 
is kept at the desired temperature by circulating the 
liquid through a thermostat. The resistance element 
N is connected to a Sunvic resistance thermometer 
controller type RT. 2 operating in conjunction with 
the thermostat. In this way the temperature of the 
liquid is maintained constant to within +0.1°C at 
temperatures between —60°C and +60°C. The 
temperature of the liquid bath is measured with 
mercury or alcohol in glass thermometers and is 
assumed -- after about one hour has been allowed 
for equilibrium to be established — to be identical 
with the temperature inside the cylinder. 

The whole assembly is well insulated thermally 
and is supported by two pieces of angle steel which 
in turn are mounted on a rigid wooden framework. 
O-rings are placed in annular grooves in the base 
and top plate to accomplish a pressure tight 
seal. The O-rings are made of synthetic rubbers 
especially selected for the substances to be measured. 
Alternatively, polythene or P.T.F.E. rings can be 
used. All other parts are sealed by suitably com- 
pressed washers made of either P.T.F.E. or poly- 
thene. It should be noted that due to the relatively 
high thermal expansion of the materials involved, 
all these gaskets are only effective over a limited 
temperature range unless they are readjusted at 
different temperatures (which is impracticable in the 
present arrangement) or are spring loaded as shown 
in Fig. 8 where the drive is led through the top plate. 

Two stainless steel tubes enter the apparatus 
through the baseplate. One tube ends inside the 
baseplate and is normally used for feeding in or 
extracting the substances to be measured. The other 
projects inside the cylinder through a hole in the 
piston and is terminated immediately below the top 
plate. This second pipe line is connected to a 
pressure gauge and can be used alternatively for 
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evacuating the apparatus by means of a rotary pump 
or for pressurising the system from a high pressure 
nitrogen supply via a reducing valve. This provides 
a convenient means for producing a constant work- 
ing pressure well above the saturation point, and over 
a temperature range, without having to fill the entire 
system with the liquid under investigation. It is ne- 
cessary, however, to take a number of test measure- 
ments before and after the admission of the nitrogen. 
in order to ensure that any nitrogen, which may be 
dissolved in the liquid, does not unduly influence 
the measured absorption. Needle valves are incor- 
porated in both pipelines to disconnect the apparatus 
from the external supply system. 


2.3-5. Che eleetrical system 


A block diagram of the electrical system employed 
in the frequency range 1 to 20 Mc/s is shown in 
Fig. 9. The trigger pulse generator has a repetition 
rate of 500 per second and produces trigger pulses 
for the time base of the oscilloscope and for the 
pulse generator incorporated in the transmitter. The 
latter provides switching pulses of 8 to 40 us dura- 
tion and of variable delay. The main part of the 
transmitter is formed by an EF 91 Colpitts oscillator 
and a QV 04/7 power amplifier. Both these stages 
can be pulsed independently and are, in normal 
operation, cut off during the intervals between trans- 
mitter pulses. The two stage transmitter is chosen in 
preference to a single stage pulsed oscillator because 
it enables accurate measurements of ultrasonic 
velocity to be carried out at low frequencies in the 
following manner: The oscillator is operated con- 
tinuously, while the power stage is pulsed. A certain 
amount of CW breaks through the power stage and 
enters the receiver, there interfering with the signal 
pulse which is received in the usual way. If a measure- 
ment is taken over an ultrasonic path of, say, 10 cm 
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and the corresponding number of beats (about 300 
at 3 Mc/s) is counted then, in principle, the ultrasonic 
wavelength can be measured with an accuracy of 


' +3 parts in 104. In practice, however, the accuracy 


of a velocity measurement is limited by the frequency 
stability of the oscillator and is estimated to be +2 
parts in 10%. This procedure for measuring velocity 
is not applicable at frequencies above 5 Mc/s, be- 
cause sufficient CW then breaks through the power 
stage to overload the receiver permanently. 

Radio frequency pulses of approximately 100 V 
peak amplitude appear at the output of the power 
stage and are fed to the transducer via a matching 
network. The latter is mounted in a small box as 
near as possible to the mechanical system without 
interfering with controls or thermal insulation. An 


H,, (circular) piston attenuator is placed between 


the transducer and the receiver input. This arrange- 
ment is chosen, because at low frequencies pulses 
with a duration of up to 40 us have to be used and 
sufficiently long time intervals between transmitter 
pulse, received echo and subsequent reflections are 
not normally available for the accommodation of a 
comparison pulse. 

Two stages of pre-amplification are used which 
provide a net gain of 25 dB, sufficient to compensate 
for the insertion loss of the attenuator. No deviations 
from linearity can be detected, within the experimen- 
tal accuracy, for input voltages up to about 0.1 V. 
Thus if 100 uV is considered as a minimum signal 
voltage, a range of 60 dB is available for measure- 
ments. The GEE receiver type R1355 is used as a 
receiver for frequencies above 10 Mc/s. This receiver 
is slightly modified to improve its overload perfor- 
mance and recovery time and plug-in RF-units have 
been adapted to cover the required frequency ranges. 
For frequencies below 10 Mc/s a T.R.F. receiver is 
used which consists of three tuned stages and three 
stages of wide band amplification. The tuned stages 
operate at the five spot frequencies 1, 3, 5, 7 and 
9 Mc/s. The low frequency receiver is completed by 
a diode detector and the video amplifier in the 
R1355 receiver. 

At frequencies above 20 Mc/s the system remains 
essentially the same, but the transmitter is replaced 
by a self-excited push-pull oscillator driven from a 


high power pulse modulator as described in section 
2.1.2. and illustrated by Figs. 4 and 5. 


2.34. Recovery after the transmitter 
pulse 


Using a transmitter output voltage of 100 V the 
receivers described in section 2.3.3. recover in 
about 5 us, which at low frequencies is much less 
than the mechanical recovery time of the transducer. 
No further improvement is obtained by ‘blanking’ 
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the first two receiver stages and this is only per- 
formed as an occasional check. In connection with 
the mechanical recovery time of the transducer, three 
different effects have to be distinguished: 


(i) The decay of the thickness vibration of the 
transducer. 
(ii) The decay of cross modes having resonance 
frequencies near the resonance frequency of 
the thickness mode and excited simultaneously 
with the latter. 
Interaction between the transducer and the 
crystal holder. 


(iii) 


The first of these effects is always present and the 
resulting recovery time depends mainly on the 
fundamental frequency, /), of the transducer, but 
very little on the actual harmonic frequency, at 
which it is operated. Practical values lie around 
50 us for fy=1Me/s and 10 us for fy=5 Mc/s, in 
fair agreement with theoretical expectations. Ab- 
normally long and irregularly shaped decays are 
produced by the second and third effects which are 
generally more troublesome than the first one and 
may, under certain conditions, make accurate meas- 
urements impossible. In order to study these effects 
and to find ways for their elimination, a number of 
different transducers have been used under different 
mounting conditions, and the following points have 
emerged from practical experience: 

1. At any given frequency of operation a trans- 
ducer with the highest practicable fundamental fre- 
quency is used because the recovery time due to both 
the thickness mode and to spurious modes decreases 
with increasing fundamental frequency, fy. of the 
transducer, 

2. When a fundamental frequency of 1 Mc/s is 
considered essential, barium titanate ceramic trans- 
ducers are used in. preference to X-cut quartz 
crystals, because the former have been found to ex- 
hibit less tendency for excitation in cross modes, and 
shorter decay times of the thickness mode. In ad- 
dition, electrical matching is much more easily 
achieved at frequencies up to about 11 Mc/s. The 
high capacity of the transducers presents some dif- 
ficulty at higher frequencies but, if special care is 
taken to avoid unnecessary line inductance, satis- 
factory signal strength can be obtained up to 15 Mc/s. 
No attempt has been made to use barium titanate 
transducers with fundamental frequencies higher 
than 1 Mc/s because of the inconveniently high 
capacities involved. 

3. Three different arrangements of the exciting 
electrodes have been used: 

a) The two opposite surfaces of the transducer are 
plated in full and each surface forms one elec- 
trode; 
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b) One of the electrodes covers only a central area 
of one of the surfaces, the other surface being 
plated in full; 

c) As with (b) except that on the surface which is 
completely plated the plating is extended around 
the edge of the transducer on to the outer part 
of the opposite surface occupied by the smaller 
electrode. 

Arrangement (a) produces a negligible electric 
field in the radial direction and is preferable for 
transducers having a low value of fundamental fre- 
quency (e. g. barium titanate at 1 Mc/s, X-cut quartz 
at 1 and 3 Mc/s), where there is strong excitation of 
cross-modes owing to the comparatively large ratio 
of thickness to diameter. Any of these arrangements 
can be employed with transducers of higher funda- 
mental frequency but (c) is preferred since it re- 
duces the influence of the clamping. 


2.3.5. Limitations of the system 


A low frequency limit is imposed on measurements 
of the absorption in liquids by the pulse method as 
a result of diffraction effects associated with the 
source of finite diameter (= 2a). Pinkerton |1] has 
found experimentally that with a common T—R 
system operating above 7.5 Mc/s no correction to 
the measured absorption is required provided the 
path length is confined to the Fresnel region of the 
sound beam and does not exceed a?/2 1. Recent cal- 
culations of Sexr et al. [13] show that the recorded 
output of the transducer should fall slightly and also 
fluctuate inside the Fresnel region, the last of a 
series of maxima occurring at a path length equal to 
2.4a?/2. The calculations show that these effects 
would have little influence at the higher frequencies 
of Pinxerron’s experiments but below about 5 Mc/s 
they become increasingly important. Preliminary 
investigations on the present system indicate that 
diffraction effects can be taken into account at fre- 
quencies above 1Mc/s and do not constitute a 
fundamental limitation to the measurements. 

A maximum change in ultrasonic path-length of 
20cm is possible with system C and this requires 
a volume of liquid of 250 cm*. If a minimum signal 
range of 10 dB is considered necessary the ab- 
sorption in the liquid must exceed 0.5 dB/cm. 

Finally, an upper frequency limit is imposed on 
the measurements. This depends upon the precision 
of the reflector movement and its accurate mea- 
surement and is further determined by the ab- 
sorption in the liquid. In practice the maximum 
absorption which can be measured with reasonable 
accuracy is about 250 dB/cm. 

The above considerations can be illustrated by 


considering as a typical example carbon tetra- 
chloride in which a/f?=550-10717s?cm™! at 
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FORMATION OF AIR BUBBLES IN AIR-SATURATED WATER 
AT REDUCED, PRESSURE AND THEIR INDICATION 


BY AN ACOUSTICAL MEASURING PROCEDURE 
by H.J. Naaxe, K. Tamm, P. Dimmic and H. W. Hersere 


III. Physikalisches Institut der Universitat Gottingen 


Summary 


With the aid of sound-absorption measurements the formation of bubbles in water 
saturated with air is investigated. The bubbles occur when the static air pressure in the 
sample holder is reduced gradually or in steps. The formation of bubbles is reduced when 
the water is subjected to a pressure greater than atmospheric before commencing the 
experiment. 

The behaviour of air pockets adhering to solids and of air-filled pores in the walls of the 
container or in solid particles under conditions of increased or of decreased pressure is 
considered. Calculations of the growth and ascent of bubbles and the sound absorption 
caused thereby are compared with measurements. 


Sommaire 


On a étudié la formation de bulles dans l’eau saturée d’air, en recourant 4 des mesures 
d’absorption du son. Les gouttes se forment lorsqu’on réduit, graduellement ou de facon 
discontinue, la pression statique de l’air sur l’éprouvette; il y en a moins si, avant ]’expé- 
rience, l’eau a été soumise a une pression supérieure a la pression atmosphérique. 

On s’est penché sur le comportement des poches d’air accolées aux solides, et a celui 
des pores remplis d’air dans les parois de l’éprouvette ou des particules solides, lorsque la 
pression est augmentée, ou bien diminuée. On a calculé la croissance et |’ascension des 
bulles, ainsi que l’absorption du son qui en résulte, et l’on a comparé les résultats de ces 
calculs avec ceux des mesures. 


Zusammenfassung 


Mit Hilfe von Messungen der Schallabsorption wird die Blasenbildung in luftgesattig- 
tem Wasser untersucht, die auftritt, wenn der statische Druck tiber dem MeBgefa8 allmah- 
lich oder in Stufen reduziert wird. Die Blasenbildung wird verringert, wenn man yor 
Beginn des Versuches einen Uberdruck auf das Wasser einwirken 1]aBt. 

Weiterhin werden Uberlegungen iiber das Verhalten von Luftansammlungen an Fest- 
kérpern und von luftgefiillten Poren in der GefaBwandung oder in festen Partikeln bei 
Anwendung von Uber- oder Unterdruck angestellt. Berechnungen iiber das Anwachsen und 
Aufsteigen von Luftblasen und die dabei zu erwartende Schallabsorption werden mit den 
Messungen verglichen. 
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1. Introduction 


Under certain physical conditions, e.g. a rapid 
decrease of pressure, vapour or gas bubbles are form- 
ed in liquids and may give ‘rise to considerable 
changes in their physical behaviour. A change in the 
absorption and velocity of sound is an example of 
this. If the physical properties of the liquid such as 
pressure, temperature, surface tension and degree 
of saturation with gas — all of which affect the 
formation of bubbles — are known, it is relatively 
easy to describe qualitatively the behaviour of 
bubbles already existing. Within certain limits it 
is possible, in addition, to calculate the rate at which 
the size of the bubbles changes. On the other hand 
no satisfactory explanation of the mechanism of 
origination of gas or vapour bubbles in a liquid has 
yet been found. 


This last mentioned problem has already been 
investigated from various directions. Becker and 
Dorine [2], [4], Turnsurz and Fisuer [21] and 
others have attempted to give a theory explaining 
the origination of a bubble or the cavitation of a 
liquid from a purely thermodynamical viewpoint. 
The conclusion of all these considerations is that a 
gas or vapour bubble which appears spontaneously 
will again vanish immediately unless considerable 
heating or high tensional stresses occur or unless the 
liquid is highly supersaturated with gas. In practice, 
however, formation of bubbles may already occur at 
much lower temperatures, without the application of 
tension and with slight supersaturation. 


The new point of view is obtained by taking into 
consideration the properties of solid surfaces which 
are in contact with the liquid either as a wall of the 
containing vessel or as a solid impurity. For every 
combination of solid and liquid an angle of adhesion 
exists. S. Taxacr [18] has shown that the negative 
pressure required for the formation of a bubble on 
a wall or on a solid particle depends on the adhesion 
angle. In the case of an ideal non-wetting liquid 
(angle of adhesion 180°) spontaneous accumula- 
tions of gas occur at the boundary between liquid 
and solid, according to his theory. Assuming an angle 
of 179°,.a tensile stress of 0.2 kg/cm? is sufficient 
for the formation of bubbles. The angles of adhesion 
observed in practice are, however, without exception 
smaller and call for appreciably higher tensile stresses. 


In another group of investigations it is assumed 
that accumulations of gas are already present in the 
liquid before the negative pressure is applied. This 
conception was already arrived at by C. Tomiinson 
[19], [20] whose publication dating back to 1867 
is probably the oldest systematic investigation of the 
problem. Tomurnson observed that bubbles were 
formed on an immersed sample when the latter was 
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in an ,,active“ condition for the formation of bub- 
bles. In the case of some materials treatment with 
alcohol or heating until red-hot destroyed the acti- 
vity, whereas in the case of others (e.g. waxes, 
resins) it did not. Tomirnson concluded that the 
active condition was probably caused by the pre- 
sence of adsorbed gas films. 

The view that accumulations of gas on solids 
(walls of the container, impurities) contribute to 
the formation of bubbles was also held at a later 
date by E. Meyer and E. Sxuprzyx [13]. SxuprzyK 
regards, in particular, small air-filled cavities in 
the bounding surfaces as being the nuclei for the 
formation of bubbles. This possibility will later be 
examined in greater detail. 

Taxacr also points out that the wetting effect, 
which he examined very thoroughly, is not by itself 
sufficient to explain all the phenomena that occur. 
In addition, effects due to the micro-structure of the 
walls have to be taken into account. The effect that 
pores in the walls of the vessel have on the resistance 
of ihe liquid against cavitation is examined by 
Fisuer [7]. He considers two idealized cross-sections 
of pores and calculates the negative pressure re- 
quired for the formation of a gas bubble. 

Another possibility, suggested by Fox and Herz- 
FELD [8], is that accumulations of air in liquid may 
be stable over long periods of time. According to this 
suggestion organic impurities in water, in particular 
proteins, whose presence can never be fully avoided, 
may form a film around an air bubble and thereby 
prevent its complete solution. According to the 
authors there are therefore always sufficient nuclei 
present in the liquid that can grow to macroscopic 
bubbles when a negative pressure is applied. 

R. B. Dean [3] developed a theory for the forma- 
tion of bubbles that does not require the presence 
of nuclei foreign to the liquid. He suspects that 
turbulence phenomena caused by impacts or sound 
waves play a decisive role and provide the forces 
necessary to overcome the surface tension. According 
to the theory rotation-free eddies are capable, at 
their moment of origination, of creating large tensile 
stresses in the liquid since they are accompanied by 
a negative pressure proportional to 1/r? (r = distance 
from eddy centre). Sound waves of relatively low 
intensity may, according to Dzan, give rise to ap- 
preciable tensile stresses in the liquid if viscosity ef- 
fects cause turbulence in the boundary layers. 

The object of the investigation to be described 
was to contribute to the explanation of the formation 
of bubbles and to investigate the further behaviour 
of the bubbles in the liquid. Certain considerations 
of G. Kurrze [10], which are based on conceptions 
similar to those of SxuprzyK, served as basis for 


this work. They will be briefly described here. 
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2. The behaviour of air bubbles in water 
containing air 


If an air bubble is floating freely in water, a con- 
tinuous exchange between the air of the bubble and 
the air in the water occurs due to diffusion. The 
preferential direction of the diffusion process de- 
pends on the partial pressure p; of the air in the 
bubble and the partial pressure p, of the air dis- 
solved in the water. The partial pressure of the air 
inside the bubble is equal to the total pressure in the 
bubble if the vapour pressure of the liquid can be 
neglected, and is made up by the hydrostatic pres- 
sure Po and the additional pressure 2-a/r due to the 
surface tension 4: p;=py+2a/r. The partial pres- 
sure p, of the dissolved air is dependent on the de- 
gree of saturation. 

In undersaturated water bubbles of all sizes dis- 
solve, whereas in oversaturated water (pa>po) an 
equilibrium radius 


2a 


Re 
Pa Po 


exists, for which the partial pressures of the air 
inside and outside the bubble are equal. This equili- 
brium is, however, unstable, i.e. should the bubble 
radius increase, the internal pressure decreases, with 
the result that air diffuses into the bubble and the 
radius is increased still further (OV/dr>0 where 
V is the volume of the bubble). A limit to this in- 
crease exists in practice in that the bubble rapidly 
rises to the surface of the liquid, and bursts. Should 
the radius of a bubble that is in equilibrium de- 
crease slightly, the internal pressure increases, air 
diffuses out of the bubble and the latter dissolves 
completely within a few seconds. One may say, there- 
fore, that free bubbles cannot exist for long either 
in under- or in supersaturated water. 


Fig. 1. Behaviour of an air-filled pore with the radius 
r of the opening; (a) pp pa, (b) Po > Pa- 


Different relations exist for bubbles clinging to 
solid surfaces. Consider an air-filled pore in the con- 
tainer wall, with an opening assumed to be circular 
with radius r’ (Fig. 1). The spherically shaped air- 
water boundary that stretches across the pore open- 
ing can have positive or negative radii of curvature. 
The absolute values of these cannot, however, be less 


H. J. NAAKE et al.: FORMATION OF AIR BUBBLES IN. WATER 


ACUSTICA 
Vol. 8 (1958) 


than the radius r’ of the pore 
Rees eee 


Equilibrium of the air diffusion exists as above, 


when the boundary surface has a radius of curvature 
ples 


In contrast to the case of a free bubble, the 
equilibrium is now stable over the range 


[|R|>r 


since here the radius decreases when the volume 
increases (OV /Or <0). An equilibrium condition can 
therefore exist for undersaturated water also, when 
the boundary surface curves slightly into the pore 
(see Fig. 2). The reason for the stability is that 
with decreasing air volume the additional pressure 
due to surface tension also decreases, thereby de- 
creasing the diffusion of air into the water. How- 
ever, if the radius of the pore opening is greater 
than the absolute value of the equilibrium radius, 
no equilibrium can exist. 


Fig. 2. Schematic representation of the course of the 
total pressure p’ in a free air bubble — — — 
and in an air-filled pore as functions of 
the volume. Stable equilibrium is possible for 
dp’/dV > 0. a surface tension, r’ radius of the 


pore opening. 


In the case of large external pressures (pp > Pa + 
+2a/r’) the entire pore is filled with water and is 
lost as nucleus, provided that no part of the wall 
inside the pore has a radius of curvature less than r’. 

With sufficient reduction of the external pressure 
(Po<Pa—2a/r’) the bubble swells outwords, be- 
comes more than hemispherical and than grows to 
such an extent — since now the radius of curvature 
is increasing — that a large part of the air volume 
is detached and rises to the surface as a free bubble. 
The remaining part continues to grow, i. e. according 
to this mechanism the nucleus is continually pro- 
ducing bubbles and is continuously extracting air 
from the water through diffusion. This phenomenon 
may be observed very well when a bottle of soda- 
water is opened, it frequently occurs that a con- 
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tinuous stream of bubbles rises from a definite point 
in the wall or in the liquid. 

In view of the behaviour of the pores under con- 
ditions of excess and reduced pressures as described 
above, it is to be expected that the application of an 
excess pressure will render the water insusceptible to 
degassing by pressure reductions up to the value of 
the initial pressure excess. This idea has already 
been mentioned by Sxuprzyx [13]. 

A further possibility for the stabilization of air 
accumulations may be thought of as being provided 
by structures with smooth surfaces and cohering to- 
gether due to adhesive or electrostatic forces. If one 
assumes that at the point of contact both particles 
have a convex curvature and that an air cushion is 
present, no pressure increase whatsoever will theo- 
retically cause complete solution of the air cushion. 
Such a model could possibly be used to explain the 
micro-photos published by Lizsermann [11], in 
which the pore-model cannot be described due to the 
minuteness of the solid particles visible in the 
picture. 

One may also, in addition, imagine a solid particle 
having a rough surface as being surrounded by a 
layer of air. In this case also, the small, negative 
radii of curvature that arise when the air-water 
boundary is stretched across the surface roughnesses 
result in a stabilization of the surrounding air layer. 


3. Methods of measurement 


Two methods were employed to investigate the ef- 
fects occurring during the degassing of water: an 
acoustical and an optical method. Only the acoustical! 
method will be described here. A description of the 
optical method will be given in a later paper [22]. 

In the acoustical method, reverberation 
measurements were carried out to measure the sound 
absorption, thereby enabling the presence of air 
bubbles to be determined. Whereas dissolved air 
causes no measurable increase in sound absorption 
[14] the formation of bubbles immediately causes 
an increase in absorption. The increase depends on 
the number and the size of the bubbles in the vessel 
as well as on the freqency of the sound wave. The 
acoustical method has the advantages that even a 
few bubbles produce an effect that is easily measured 
and that the degree of damping can easily be re- 
corded over fairly long periods of time. Its dis- 
advantage is that the damping is determined by two 
independant quantities — number of bubbles and 
bubble radius — and that a further measurement is 
required to determine one of these before an un- 
ambiguous result is obtained. 

To carry out the absorption-measurements by the 
reverberation time method, the liquid sample under 


time 
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test was poured into a glass vessel suspended inside 
a pressure-flask. A piezo-electric crystal, attached to 
the outer side of one wall of the vessel and driven 
by pulsed oscillations, was used to excite numerous 
characteristic modes of vibration in the liquid. The 
decaying reverberation that commences when the 
driving signal is switched off, was picked up by a 
second crystal similarly attached to the wall. This 
receiver voltage was given a logarithmic characteris- 
tic, rectified and applied to the vertical deflection 
amplifier of the display C. R.O. which had a cali- 
brated constant-velocity sweep. The decay constant 
in dB/second of the entire system, consisting of con- 
tainer plus liquid, could be read off directly from 
the slope of the straight line drawn on the C. R. O. 
screen. Similar measurements with water that was 
completely devoid of bubbles were always carried 
out in order to eliminate the constant losses due to 
the containing vessel, thereby enabling the effect of 
oscillating bubbles to be determined. 


4, Preparatory treatment of the water 


A large number of preparatory experiments were 
necessary before reproducible conditions could be 
obtained. Distilled water was used, which was drawn 
through a very fine glass filter. Subsequently, a large 
quantity of the water (1801) was poured into a 
vessel and saturated with air. The air used for this 
purpose was filtered through cotton-wool to remove 
dust particles etc., and then allowed to bubble 
through a distributing filter. 

Two methods were used to determine the air con- 
tent of the water. In the first method, the air con- 
tained in a given quantity of water was displaced 
by the same quantity of a saturated potassium 
hydroxide solution, and its volume was measured. In 
the second method, the air contained in a known 
volume of water was driven off by carbon dioxide 
produced within the water from calcium carbonate 
and hydrochloric acid. The volume of air was meas- 
ured after the carbon dioxide had been removed by 
a concentrated potassium hydroxide solution. 

The air content was always adjusted so that the 
water was just saturated before the experiments 
began. The temperature of the water was kept at 
about 20°C. After being saturated with air, the 
water required was transferred through a glass tube 
to the testing vessel. To the water, hitherto cleaned 
of or kept free from contaminating substances, a 
weighed quantity (100 mg) of defined ,,contaminat- 
ing“ particles was added. For this purpose jeweller’s 
rouge was used at first, and later a powder was 
employed consisting of pure aluminium oxide and 
having a mean particle diameter of approximately 
4. um. 
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5. Calculation of the expected damping in the 
measuring vessel 


5.1. Damping caused by one bubble 


The calculation of the damping caused by a 
bubble is based on a theory due to Prriem [15]. 
This theory has been confirmed by measurements 
carried out on single bubbles by M. L. Exner [5], 
W. Hamer [6] and H. Hassxe [9]. 

An air bubble in water represents a system capable 
of performing oscillations. The compliance is due 
to the air enclosed in the bubble, the mass due to 
the medium taking part in the oscillation. The re- 
sonant frequency is given by the relation 


i 1 seen, (1) 


~ Oar Ow 


where y =c,)/c, is the ratio of specific heats and @, 
the density of water. 

For a constant radius, the resonant frequency is 
therefore proportional to Vpy+2a/r. The damping 
due to the bubble has two components: a thermal 
component, caused by the non-adiabatic compres- 
sion of the air in the bubble, and a radiation com- 
ponent due to the spherically symmetrical radiation 
of energy into the surrounding medium. Both com- 
ponents in combination essentially determine the 
amplitude of oscillation of the bubble near its re- 
sonant frequency. An irreversible loss of acoustical 
energy in the container is, however, caused only by 
the thermal component, since in the diffuse sound 
field existing in the vessel the reradiated energy is 
to all intents and purposes equivalent to that pre- 
viously present. 

The energy absorbed per unit time by thermal 
processes by the bubble has a value of 


where vy is the R. M.S. velocity at surface of bubble, 
F = dr/S op the compliance constant of bubble, S the 
surface area of bubble and 7, the loss factor of 
thermal damping. 

The relation between the velocity at the bubble 
surface vg and the free-field sound pressure p is 
given by the impedance of the bubble: 


1 eee 
Z= Favs (Mth + Nr) +i(om- =] (2) 


(where 7, is the damping factor due to radiation and 
M=4x 0,7 the oscillating mass) so that the re- 
lation between the thermally absorbed energy and 
the sound pressure in the container is (3) 


Nth SPs 


pease = —_—— ——- 2 
 (/o F)? (nn +1)? + (@ M—1/0 F)? wo FP 
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If one now assumes that per unit volume one bubble 
of radius r is present and one introduces the acousti- 


_cal energy present per unit volume at a sound pres- 


sure of po, the decay constant of the entire system 
is 
Nth 
4 — a 
(Nth + Mr)? + (@o?/ a9” — 1)? 
x S? w F (Oy Cy?) dB/s, (4) 


where (@y Cw") is the compressibility of water. This 
assumes that the losses in the water itself are negli- 
gible, and that the sound energy in the bubbles is 
negligible in comparison to the sound energy in the 
water. 

That the first condition is fulfilled, is immediately 
evident when the damping of distilled water quoted 
in the literature (0.8 dB/s at 50 kc/s) is compared 
to the damping due to air bubbles that was measured. 
(Order of magnitude 100 to 1000 dB/s.) Since 
furthermore a bubble by itself, with a resonant fre- 
quency of 50 kc/s, should have a decay constant of 
150,000 dB/s the only explanation that can be of- 
fered for the much lower values of damping con- 
stant actually measured is that by far the greater 
part of the energy is stored outside the bubble, i. e. 
in the water. Hence the second condition is also ful- 
filled. 

To apply the above equation for the decay con- 
stant to bubbles that are undergoing a slow growth, 
it is also necessary to know the damping factors 7p 
and 7, for bubbles excited at a frequency other than 
the resonant frequency. Prriem [15] gives the value 
for the thermal damping factor 7, as function of 


6=4.,3 


the dimensionless quantity y=2rVxf/a, where 
a=A/cy es the thermal conductivity of air. The mass 
of the medium participating in the oscillation does 
not enter into Prriem’s calculation, so that the results 
are also applicable to bubbles excited outside their 
resonant frequency. 

The loss factor due to radiation may easily be 
calculated from the compressional energy in the 
bubble and from the power radiated by the bubble 
as a spherical radiator [17]. We obtain 


_ 823 ev(foR)? F 
3 Cy Oa Ps fo® ; 
where 0,, Cw are the density and sound velocity 


respectively for water, @, and c, the density and 
sound velocity for air, and R the resonance radius. 


Nr (5) 


5.2. Damping caused by a given number of bubbles 


The optical investigations of the degassing process 
have shown that under the test conditions that ex- 
isted the concentration of bubbles was always less 
than 1.5 per cc. The absorptional cross-section of a 
bubble at resonance, on the other hand, is only 
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0.5 cm? at 50 kc/s and for bubbles not in resonance 
it is even less. Therefore, and also because of the 
great difficulties encountered when the problem is 
treated theoretically, the mutual reaction of bubbles 
on one another was not considered, even though 
under unfavourable circumstances the statistical 
distribution of the distances between bubbles may 
result in an effect. The increase in the decay constant 

of the system, caused by the bubbles, is therefore 
given by the expression (4) in the case of one 
bubble per cm’, and is proportional to the bubble 
concentration. 


5.3. Rate of growth of the air bubbles 


If the air pressure over a given volume of water 
saturated with air is suddenly or gradually reduced, 
the water becomes supersaturated with air. The de- 
crease in the hydrostatic pressure results, on the one 
hand, in an increase in the volume of the gas accu- 
mulations existing as thin layers, or in the pores of 
solid particles or stabilised in some other way by 
their shape. In many cases this results in an increase 
in the radius of curvature of the air-water bound- 
aries. On the other hand, the critical radius 


2a 2a 


R= Lis 
Pa—Po -4p 


(6) 


is decreased because of the increase in Ap. Both 
effects have as result that some air accumulations 
may increase beyond the equilibrium point, and 
from there on grow still further due to diffusion. 
If during this growth the radius does not again be- 
come less than the critical radius (e.g. when the 
opening of an air-filled pore is sufficiently large) 
the growth will continue steadily. The rate of growth 
may be calculated if the diffusion constant of air in 
water is known. 

Measurements of the diffusion constant have 
hitherto only been carried by observation of bubbles 
dissolving in undersaturated water. The bubbles con- 
cerned were either clinging to a wall, or rising free- 
ly in an observation tube, dissolving as they rose. 
Lizesermann [11] observed that the rate of growth 
was larger for freely rising bubbles than for bubbles 
attached to some surface. He obtains a diffusion con- 
stant apparently 2 to 3 times larger for freely rising 
bubbles using in both cases the formula (7) given 
below for resting bubbles. Similar results, i.e. an 
apparent diffusion constant (D’ =5.4-10~* cm?/s) 
2.7 times as large as the known value (D = 2.0-10~° 
em?/s) for diffusion of air in resting water, have 
recently been observed by Marrini [12] who also 
gives a simple model for explaining this fact. 

If one assumes that the diffusion constant observed 
when a bubble dissolves also holds good for the case 
of a bubble growing in supersaturated water, the 
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size of the bubble can be expressed as a function of 
time. Disregarding the effect of surface tension for 
the present, we find for resting bubbles 


Pale eGo (7) 


where rg is the initial radius, D the diffusion con- 
stant, s the solubility of air in water and 0 the 
degree of supersaturation. If surface tension is taken 
into account, one obtains (see ref. [10]) (8) 


rare R(G6-+2)(r—r9+ Rin) = S Oks 


where & is the critical radius. 

It is important that for small bubbles the surface 
tension be taken into consideration, since it reduces 
the rate of growth of small bubbles. In both equa- 
tions the initial radius of the bubbles at the time 
t=O appears. This is, to begin with, an unknown 
quantity. If, however, one reduces the external air 
pressure in small steps (e.g. by 100 mm Hg at a 
time) the range of values assumed by the initial 
radii may be confined. Additionally one must make 
the plausible assumption that in the time between 
two pressure steps the radii of nuclei with r<R do 
not alter appreciably. The bubble size as function of 
time for the individual pressure steps is given in 
Fig. 3. For the time t=0 the value r=2R was 
assigned arbitrarily. Deviations of the initial radius 
from this value result merely in a shift of the 
abscissa. From'these curves it is possible for example, 
to determine the time required for the bubble to 
become just visible (r=30um) or to attain the 
resonance radius (e.g. 60 um at 50 ke/s and 600 mm 
Hg, see Fig. 4). 


10? | 


r/R 


Fig.3. Growth of bubbles by diffusion. Radius r of 
the bubbles, referred to the critical radius 
R=2a/(pa—po), as a function of the time t, 
for different pressures pp in water saturated for 
Pa=750 mm Hg, calculated for D’=5,4- 107° 
em?/s. ro/R=2 has been chosen arbitrarily. 


5.4, Velocity of ascent of the bubbles 


The frictional force acting on a solid sphere mov- 
ing in a viscous medium has been calculated by 
Stokes as 


K=6nxvrv (9) 


(y: viscosity). 
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Fig. 4. Time t necessary for the growth of a bubble up 
to its resonance radius starting with a radius 
between two limiting values corresponding to 
r/R(po) =1.05 and r/R(p)+100 mm Hg) =1 
resp. calculated for the normal and for the lar- 
ger apparent diffusion constant i.e. for resting 
and for freely rising bubbles. x measured at 
50 ke/s, © measured at 100 ke/s. 


According to Osren and Gotpstern this formula is 
the first term of a series of powers of the Reynolds 
number. Deviations for larger Reynolds numbers, 
i.e. for higher velocities (onset of turbulence) may 
therefore be expected. 


Auten [1] confirmed the validity of Stoke’s for- 
mula for ascending bubbles of small radius, and 
found deviations tending in the direction of OsrEEN 
and Go.pstTEIn’s corrections for larger bubbles (r 
between 50 and 200 um): for bubbles having these 
radii the measured velocity of ascent is less than that 
calculated from Stokes’s formula. The applicability 
of Stokes’s law to bubbles is not immediately obvious, 
since it is not known if the layer of liquid im- 
mediately at the boundary of the bubble moves with 
the latter so that it acts as a rigid sphere, or whether 
it flows around the bubble and is constantly renewed. 
Rysczynski [16] therefore proposed a new formula 
for the force acting on a moving bubble. The 
measurements of ALLEN, which have recently been 
confirmed by Exner and Hamper [6] and by Lizser- 
MANN [11] appear, however, to justify the use of 
Stokes’s formula. 


5.5. Concentration of bubbles in the measuring vessel 


Because bubbles slowly rise during measurements 
and reach the surface of the liquid, the number of 
bubbles in the vessel decreases. The decrease in the 
bubble concentration can be found, since the bubble 
radius as function of time, the relationship between 
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bubble radius and rate of ascent and the dimensions 
of the containing vessel are known. 
If qo is the number of bubble nuclei that form 


growing bubbles at the application of the stepwise 


pressure reduction, v(t) is the velocity of ascent of 
the bubbles as a function of time and A is the mean 
height of the container, then the mean concentration 
of bubbles referred to the total volume of liquid is 
given by 


q(t) =49)1- 5 [tose . 


0 


(10) 


The expected variation of the concentration with 
time is shown in Fig. 5 for the two different values 
of the apparent diffusion constant (resting and ris- 
ing bubbles resp.). In each case two curves were 
computed for the two limiting values (see below) 


of the bubble nuclei. It frequently happens that the’ 


concentration does not decrease according to the 
above function, because a number of bubbles cling 
to the walls and do not rise. A degree of damping is 
then measured that does not return, even after a 
considerable period of time, to the value originally 
existing before commencement of the experiment. 


p= 500mm Hg 


0 20 40 60 80s 


Fig. 5. Relative bubble concentration in the measuring 
vessel as a function of time at the 500mm Hg 


pressure step. Initial radii of the nuclei: 
ro =1.05 R or 1.67R resp., calculated for the 
normal (— — —) and for the larger ( ) 


apparent diffusion constant corresponding to 
resting and to freely rising bubbles. 


Starting with the known initial bubble concen- 
tration c; (determined by optical methods [22] and 


given in the following table) the course of the damp- — 


ing in the measuring vessel then can be calculated 
by the aid of the curves given in Fig. 5: 


P 600 | 500 
Cee wehbe 


- 400 | 300 200 | 100 
0.64 | 0.31 | 0.36 | 0.36 


mm Hg 


1/em$ 


6. Measurements 


Two methods of pressure reduction were employed 
for the acoustical measurements. In the first, the air 
pressure was continuously reduced with an oil pump 
from normal (approx. 750 mm Hg) to about 100 mm 
Hg within a time of approximately 120 seconds. 
The rate of decrease of pressure was, in this case, 


: 


© aie 
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determined by the capabilities of the pump and the 
volume of the pressure flask. In a second series of 
measurements, the pressure was reduced in one 
sudden step first to 700mm Hg and then in six 
further steps of 100 mm Hg to a final pressure of 
100 mm Hg. By connecting the pressure flask briefly 
to an evacuated reservoir of adequate volume, the 
pressure reductions could be effected in less than 
one second, between each of the individual pressure 
steps a pause of a few minutes allowed completion 
of the degassing process that occurred. Before degass- 
ing begins, the damping measured in the vessel is 
determined principally by the losses in the vessel 
_ itself and in the suspension. With the set-up used, 
this damping was about 150 dB/s. When bubbles 
begin to form, this damping increases more or less 
strongly. With the apparatus available it was pos- 
sible to measure rates of decay up to 2000 dB/s in 
very rapid succession. 


6.1. Measurement of the pressure-dependence of the 
onset of degassing 


The onset of .degassing as a function of the 
previously applied pressure excess was measured 
under conditions of a steady reduction of pressure. 
The excess pressure could be varied from 0 to 1.2 
atmospheres, the upper limit being set by the 
maximum permissible loading of the pressure flask. 
The amount of air that diffused into the water 
through its undisturbed free surface during appli- 
cation of the excess pressure was so small that the 
degree of saturation of the water remained practi- 
cally unaltered. The defined impurities used in these 
measurements consisted of jeweller’s rouge and 
aluminium oxide. Measurements were carried out 


100 


50 


0 0.2 04 atm 06 


Fig. 6. Shift of the onset pressure of degassing against 
that of the trial measurements as a function of 
the applied overpressure. 
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with frequencies between 10 and 100 kc/s. The 
results obtained with aluminium oxide are shown in 
Fig. 6. In this figure the shift in the point of initia- 
tion is shown that occurs when the results of a test 
are compared to the results obtained with a trial 
measurement carried out on the same sample without 
an initial excess pressure being applied. This method 
of measurement was necessary because large fluc- 
tuations in the pressures at which damping com- 
menced were observed between different fillings of 
the measuring vessel. In spite of this precaution, 
Fig. 6 still shows a large spread in the measured 
values. Systematic differences between the different 
frequencies at which measurements were performed 
cannot be found. The direction in which the point 
of initiation shifts agrees with the prediction that 
the higher the initial excess pressure, the lower the 
pressure to which the water may be subjected before 
degassing begins. The view that excess pressure and 
pressure reduction are equal could not, however, be 
confirmed. 


6.2. Measurement of the maximum value of damping 


During evacuation the damping in the vessel 
could be measured continuously. The values of 
greatest damping occurring with jeweller’s rouge 
and aluminium oxide are plotted as a function of 
the previously applied excess pressure in Fig. 7. The 
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Fig. 7. Maximum of the measured additional damping 
as a function of the previously applied excess 
pressure. Impurity: (a) jeweller’s rouge, (b) 
aluminium oxide powder. 


measurements were carried out with frequencies 
between 10 and 200 ke/s. The effect was found to 
be strongly frequency-dependent. It is noteworthy 
that at high frequencies (200 kc/s) no measurable 
damping occurs after application of a slight pressure 
excess, whereas for low frequencies an increase in 
damping is still found even for high initial excess 
pressures. The explanation for this is that the few 
nuclei remaining after application of the excess 
pressure undergo a steady change in size and thereby 
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in resonance frequency and are less effective as small 
bubbles than as large bubbles. A larger fraction of 
the sound energy is therefore absorbed the larger 


the bubbles are, and the lower therefore their fre-' 


quency of resonance. After a high initial pressure 
the effect of the remaining bubbles is only noticeable 
when they have a large diameter, i.e. at low fre- 
quencies. At low frequencies and after the appli- 
cation of a small initial excess pressure the maxi- 
mum damping is appreciably greater than could be 
measured with the apparatus (2000 dB/s). 


6.3. Effect of the duration of the pressure excess 


One series of tests was performed to investigate 
the effect of the duration of the previous excess 
pressure on the degassing process. The degree of 
degassing was again determined from the value of 
the observed damping maximum. In every test the 
excess pressure used was 0.6 atm. The results are 
shown in Fig. 8. It will be seen that with increas- 
ing duration of subjection to the excess pressure, the 
degassing at first is reduced rapidly and then attains 
a steady value depending only on the value of the 
excess pressure. This steady value may be measured 
after a time of about 10 min, and therefore all the 
experiments in which the point of initiation and the 
maximum value of the damping were measured, 
were performed after application of an_ initial 
pressure excess of 10 min duration. 
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dB/s 
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Fig. 8. Maximum attenuation measured as a function 
of the duration of the application of an excess 
pressure of 0.6 atm. 


A possible explanation may be advanced to ac- 
count for the relatively long time required by the 
initial pressure before it becomes fully effective: It 
is possible that wetting effects inside a pore play a 
part in the process, i.e. that for complete inhibition 
of a nucleus it is necessary that certain parts of a 
pore wall be gradually wetted by the water, which 
would not reach those parts during the course of 
the deformation of the air-water boundary deter- 
mined by the ratios of the partial pressures. One 
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may well imagine this gradual progress of a film of 
water requiring times of the order magnitude of 
minutes. 


6.4. Effect of elapsed time between initial pressure 
excess and evacuation 


In general, experiments were commenced with im- 
mediately after application of the initial pressure ex- 
cess. In one series, however, the water was subjected 
to an excess pressure of 1 atm and then stored for 
15 days. It was found that nuclei, once inhibited, 
were not capable of producing bubbles even after such 
a long recovery period, i.e. no greater increase in 


damping could be observed. This result is in agree- 


ment with the idea of Dean [3] that no subsequent 
adsorption of air occurs at a surface once it has been 
wetted. 


6.5. Measurements with degassing carried out in 
steps 

To obtain closer information regarding the sizes 
of the nuclei present, a few series of measurements 
were carried out in which the pressure was not 
reduced gradually but in steps of 100 mm Hg at a 
time. For this purpose the pressure flask was 
connected to a previously evacuated container 
through a tube of large cross-section. In this way 
the pressure could be reduced almost instantaneously 
by the required amount, and then held constant 
during degassing. In addition this method possesses 
the advantage that clear distinction is made possible 
between pressure and time effects occurring during 
degassing. It then becomes possible to compare the 
measured damping curves with the calculated values. 

Very little air is extracted from the water by the 
formation of bubbles, in spite of the large effect on 
the absorption of sound. For each pressure step one 
may therefore be certain to have a nearly unaltered 
air content of the water, corresponding to a degree 
of saturation of 1 at 750 mm Hg due to the pre- 
paratory treatment. Fig. 9 shows how the additional 
damping obtained from the experiment varies with 
the time elapsed after the application of the cor- 
responding pressure step. Noteworthy in these ,,cur- 
ves“ is the time ¢, which elapses between the initia- 
tion of the pressure reduction and the onset of an 
increase in the damping. The higher the pressure in 
the vessel, the longer this time becomes. This re- 
lation is in agreement with the theoretical expecta- 
tions and it can easily be seen that the lower the 
pressure in the vessel is made, the higher the degree 
of supersaturation of the water becomes and the 
faster therefore the growth of the bubbles. The 
following table shows the times that were measured: 


P| 600 | 500 | 400 | 300 | 200 | 100 | mm Hg 


# | 20 | 10 . eee 
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Fig. 9. Course of the damping as a function of time during stepwise degassing (f =100 kc/s); 
calculated for the normal diffusion constant (D=2.0-10~* cm?/s), 
——— calculated for the larger apparent diffusion constant (D’ =5.4: 1075 cm?/s), 


© measured values. 


(—-— For the pressure step to 500 mm Hg an additional theoretical curve is shown, 
calculated with D=2.0-10~* cm?/s below and with D’ =5.4-10~* cm?/s above the 


resonance size.) 


The maximum value of the additional damping is 
determined by the number of bubbles formed in the 
vessel. When the degassing occurs in steps, the damp- 
ing is always much less than the maximum value which 
can be measured with the apparatus (2000 dB/s) so 
that the measured values may be compared to the cal- 
culated values. The number of bubbles formed during 
each step was obtained from the optical measurements. 
The resulting ,,theoretical“ curves, as calculated, are 
also plotted in Fig. 9. Two pairs of curves were com- 
puted for each pressure step: One pair using the nor- 
mal diffusion constant D, the other pair using the 
higher apparent diffusion constant D’. One of the two 
curves, each pair consists of, is calculated for the 
largest and one for the smallest nuclei in question. 
The largest nucleus radius was taken to be the 
critical radius of the previous step i. e. these bubbles 
could just not be induced to grow during the 


previous pressure level. The smaller value of radius 
was assumed to be 7/R=1.05, i.e. a radius 5% 
greater than the critical radius for the existing 
pressure level. Nuclei with a radius only just above 
the critical value (r/R between. 1 and 1.05) grow 
extremely slowly and contribute a negligible amount 
to the absorption. 

These curves would have to be applied if all the 
bubbles had either the one or the other of the two 
initial radii considered. In actual fact the initial 
radii have values lying between these limits, and 
the curve to be expected will therefore lie somewhere 
between the computed boundary curves. 

A comparison between the experimental and the 
theoretical curves shows that surprisingly the course 
of the experimental values first follows much better 
the theoretical curves calculated for resting bubbles 


(i.e. with the diffusion constant D = 2.0 - 10~*° cm?/s) 


J52 


than those for freely rising bubbles (with a larger 
apparent diffusion constant D’ =5.4-10~> cm?/s). 
This fact agrees with the results of the optical 
method [22] which show a slow growth for bubbles 
smaller than the limit of visibility (;=30 um) and 
a fast one for larger bubbles. A theoretical damping 
curve (see Fig.9, py =500 mm Hg) calculated ac- 
cording to this rule indeed shows best relative 
agreement with the experimental results. While there 
is a quite good agreement with respect to the relative 
course of the curves, a discrepancy must be stated 
between their absolute values; the theoretical ones 
being too high. A reason for this deviation cannot 
be advanced at present with certainty. It may be 
assumed the theoretical values are too high because 
of the neglection of the interaction of the bubbles 
between each other (see section 5.1). Further 
theoretical investigations seem desirable here. 


The work described in this paper was carried out 
under arrangements with the Department of Scien- 
tific and Industrial Research, London, whose as- 
sistance and interest are gratefully acknowledged. It 
is published with the permission of the Department. 
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BEITRAGE ZUR DEUTUNG. DER SCHALLGESCHWINDIGKEIT 
IN BINAREN FLUSSIGKEITSMISCHUNGEN 


von R. Grenaur 


I. Physikalisches Institut der Technischen Hochschule, Stuttgart 


Zusammenfassung 

Fu8end auf den bekannten, von M.R. Rao und W. Scuaarrs aufgestellten halbempiri- 
schen Beziehungen R=Vu-M/o und u=u.‘sB/V wird an Hand eigener und fremder 
Messungen die Abhangigkeit der Schallgeschwindigkeit von der Konzentration in binaren 
Fliissigkeitsmischungen diskutiert. Es zeigt sich, da} der Verlauf der Schallgeschwindigkeit 
iiber der Konzentration qualitativ gedeutet werden kann. Bei Mischungen, deren Kompo- 
nenten keine schweren Atome enthalten bzw. keine extrem hohe Schallgeschwindigkeit 
besitzen, la8t sich auch ein quantitativer Zusammenhang zwischen der Schallgeschwindig- 
keit und der Dichte in Abhangigkeit von der Konzentration angeben. Zum Schlu8 wird auf 
eine Moglichkeit hingewiesen, bei schalldispergierenden Fliissigkeiten die Dispersion bei 
konstanter Frequenz in Abhangigkeit yon der Temperatur zu messen. 


Summary 

The dependence of the velocity of sound on the concentration in binary fluid mixtures 
is discussed in terms of the well-known semi-empirical relations R=Vu -M/o and 
U=Uo*SB/V established by M. R. Rao and W. Scuaarrs. When the velocity of sound is 
plotted against the concentration, the result permits a qualitative interpretation. In cases of 
mixtures with components not having heavy atoms or extremely high velocities of sound, it 
is also possible to establish a quantitive interrelation between the velocity of sound and 
the density as a function of the concentration. For sound-dispersing fluids, finally, the 
possibility is pointed out of measuring the dispersion at constant frequency as a function 
of the temperature. 


Sommaire 


Le rapport entre la vitesse du son et la concentration de mélanges de liquides binaires 
est discuté sur la base des relations semi-empiriques connues, établies par M. R. Rao et 
W. Scuaarrs: R=Vu-M/o et u=U.*sB/V, et a Vaide de mesures effectuées par les 
auteurs et de mesures de sources étrangéres. I] est démonstré que la variation de la vitesse 
du son peut étre interprétée qualitativement en fonction de la concentration. Pour les 
mélanges dont les composantes ne contiennent pas d’atomes lourds et dans lesquels la 
vitesse du son n’est pas extrémement élevée, on peut également établir un rapport quan- 
titatif entre la vitesse du son et la densité en fonction de la concentration. A la fin, on 
indique la possibilité de mesurer, pour les liquides qui dispersent le son, la dispersion 
en fonction de la température, la fréquence restant constante. 


iiber 180 organischen Fliissigkeiten die Brauchbar- 
keit von Gl. (1) gepriift. 

1940 fand M. R. Rao [2] folgende, ebenfalls 
empirische Beziehung zwischen Schallgeschwindig- 
keit u, Dichte 9 und Molekulargewicht M in reinen, 


S 


1. Einleitung 


Fir die Schallgeschwindigkeit in reinen organi- 
schen Flissigkeiten hat W. Scuaarrs [1] die empiri- 
sche Formel 

U=U.'s B/V (1) 


nicht assoziierten Flissigkeiten: 


angegeben. In ihr bedeuten uw, eine Konstante fiir 
alle organischen Fliissigkeiten mit der Dimension 


R=Vu'M/o=Vu-V. (2) 


einer Geschwindigkeit, s ein die Elastizitat der Mole- 
kulst6Be bestimmender StoBfaktor, der den fiir ela- 
stische Kugeln geltenden Wert s=4 nicht iiberstei- 
gen kann, B das Eigenvolumen der Molekiile und V 
das Molvolumen. Die zunachst unbekannten Gréfen 
s und B setzen sich, wie Scuaarrs zeigt, additiv aus 
den entsprechenden, die Fliissigkeitsmolekiile auf- 
bauenden Atomgré$en zusammen. Scuaarrs hat an 


R ist eine temperaturunabhangige Stoffkonstante, die 
von Lacemann und Corry [3] den ungliicklichen 
Namen ,,molekulare Schallgeschwindigkeit* erhielt. 
Ihre Dimension ist cm!/?-s~‘*, Wie schon Rao [4] 
und auch Lacemann und Corry [3] zeigten, ist R 
eine additive GroBe, d.h., die molekulare Schall- 
geschwindigkeit einer Verbindung setzt sich additiv 
aus den sogenannten Atominkrementen (bei Lace- 
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MANN und Corry: Bindungsinkrementen) zusammen. 
Um dieser Additivitat eine Deutung geben zu kon- 
nen, geht Scuaarrs [5] von Gl. (1) aus und erhalt 
aus ihr eine GroBe mit der Dimension eines Volu- 
mens: 


S=Vulu..:-V =VsBV?. (3) 


Sie hangt mit der molekularen Schallgeschwindigkeit 
R in folgender Weise zusammen: 


S=R/Vitoo - 


In der vorliegenden Arbeit wird S als Schallvolumen 
bezeichnet. Da sich R und S nur durch einen konstan- 
ten Faktor unterscheiden, gelten alle Uberlegungen, 
die sich auf die molekulare Schallgeschwindigkeit R 
beziehen, auch fiir das Schallvolumen S, womit die er- 
wahnte Additivitat etwas plausibler gemacht ist. 
Ersetzt man in Gl. (2) mittels der Beziehung 
u® = (0 f)~+ die Schallgeschwindigkeit wu durch die 
adiabatische Kompressibilitat 6, so wird 


R=07-" B-"M 
oder 


(4) 
W = Re Mh = Bo M/o =u M/¢"". 


Die GroBe W, von Y. Wava [6] als ,,molekulare 
Kompressibilitat“ eingeftihrt, ist ebenfalls tempera- 
turunabhangig und additiv. Ihre Dimension ist 
e@-cm20+s~2, 


2. Messungen 


Den nachstehenden Uberlegungen liegt folgendes 
experimentelle Material zugrunde: 


R. GENAHR: SCHALLGESCHWINDIGKEIT IN FLUSSIGKEITSMISCHUNGEN 


ACUSTICA 
Vol. 8 (1958) 


1. Vom Verfasser wurden von den vier Mischun- 
gen 


Benzol — t-Butylalkohol (1) 
Cyclohexan — Essigsdure (2) 
Dekalin — n-Propylalkohol (3) 
Chlorbenzol — n-Propylalkohol (4) 


Dichte und Schallgeschwindigkeit als Funktion der 
Konzentration bei 20° C gemessen; die erstere mit 
einem Satz Ardometer (MeSgenauigkeit: +0,1%), 
die letztere mit einem Ultraschallinterferometer der 
Fa. Steeg & Reuter (MeBfrequenz 2 MHz, Mef- 
genauigkeit ebenfalls +0,1%). Diese MeBwerte sind 
in Tabelle | zusammengestellt. 

Die Fliissigkeiten wurden jeweils in der reinst er- 
haltlichen Form von der Fa. Merck, Darmstadt, bezo- 
gen und anschlieBend wie folgt gereinigt: 

Benzol und Cyclohexan: 

tiber Na fraktioniert destilliert, 
Dekalin: 

getrocknet mit Na, dann zweimal einfach destil- 

liert, 
n-Propylalkohol und t-Butylalkohol: 

einen Tag iiber CaO gekocht, 
Essigsaure: 

gekocht mit Essigsaureanhydrit, dann fraktioniert 

destilliert, 
Chlorbenzol: 

fraktioniert destilliert. 

2. Ferner wurden Messungen anderer Autoren 
ausgewertet, und zwar an folgenden Mischungen: 


Heptan— Dichlormethan 
Tetrachlorkohlenstoff — Athylalkohol 


(5) 
(6) 


Tabelle I 
Benzol-t-Butylalkohol 


Mol-% 0 4,10 9,87.) 21,25 2) 
w 1318,0 | 1308,7 | 1290,2 
@Q 0,8790 | 0,8738 | 0,8666 


45,80 60,20 78,98 100 


1256,0 1201,9 1175,3 1155,2 1141,3 m/s 
0,8538 0,8285 0,8143 0,8010 0,7868 g/em? 


Cyclohexan-Essigsaure 


23,60 | 36,50 


1267,0 
0,7829 


1195,4 | 1161,7 
0,8293 | 0,8754 


0,8077 


1223,7 | 


63,00 | 75,70 | 86,20 92,20 100 
1149,3 | 1136,0 | 1133,5 | 1142,0 | 1164,3 | m/s 
0,8937 | 0,9337 | 0,9804 | 1,0083 | 1,0550 | g/cm? 


Dekalin-n-Prophylalkokol 


1413,6 
0,8778 


1280,1 
| 1,1017 


1286,7 
1,1112 


1273,8 


1267,2 
1,0890 


1,0737 


24,56 | 


40,30 | 59,40 | 


1254,2 
1,0029 


1248,1 | 1240,1 


1226,2 
0,9514 | 0,8933 


0,8204 


1221,2 
0,8034 


g/cm? 


alla eal 
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Methylalkohol — Cyclohexanol (7) 
n-Propylalkohol — Cyclohexanol (8) 
iso-Propylalkohol — Cyclohexanol (9) 
n-Butylalkohol — Cyclohexanol (10) 
iso-Butylalkohol — Cyclohexanol (11) 
Oktylalkohol — Cyclohexanol (12) 
Glykol — Cyclohexanol (13) 
Glyzerin — Cyclohexanol (14) 
Diese MefBwerte stammen von Tvomikosxi und 


Nurmi [7] (Mischung 5, entnommen einer Arbeit 
von D. Serre [8], MeBtemperatur 25° C), K. Sacuer 
[9] (Mischung 6, MeBtemperatur 18°C) und F. 
_ Danusso [10] (Mischung 7 bis 14, MeBtemperatur 
30° C). Die Genauigkeit dieser Messungen ist nicht 
bekannt. 


3. Ergebnisse 


3.1. Zusammenhang zwischen u und @ in bindren 
Flissigheitsmischungen 


In bindren Mischungen gehorchen im allgemeinen 
weder die Dichte bzw. das Molvolumen noch die 
Schallgeschwindigkeit der Mischungsregel *. Irgend- 
ein Zusammenhang im Verlauf der Dichte und der 
Schallgeschwindigkeit in Mischungen konnte bisher 
nicht gefunden werden. Ist z.B. bei einer bestimm- 
ten Gruppe von Mischungen die Dichtekennlinie 
konvex nach oben gekriimmt, so konnen die Schall- 
kennlinien dieser Mischungen konvex nach oben 
oder unten gekriimmt sein, sie konnen aber auch 
einen Wendepunkt besitzen. Entsprechendes gilt, 
wenn die Dichtekennlinie konvex nach unten ver- 
lauft bzw. einen Wendepunkt besitzt ?. 

Ohne eine quantitative Annahme zu machen, kann 
man tiber den Verlauf der Schallkennlinie folgendes 
aussagen: Aus Gl. (1) sieht man, dafi die Schall- 
geschwindigkeit in Mischungen auch dann _ nicht 
linear mit der Konzentration verlauft, wenn sowohl 
das Molvolumen V als auch das Volumen s BP der 
Mischungsregel gehorchen. Bei bekanntem Verlauf 
des Molvolumens wird der Charakter der Schall- 
kennlinie nur noch durch den Verlauf des Volumens 
s B in der Mischung bestimmt. Uber den StoBfaktor 
s kann man von Fall zu Fall qualitative Aussagen 
machen. Er wird im allgemeinen kleiner sein als es 
der Mischungsregel entspricht, und zwar um so klei- 
ner, je mehr sich die Molekiile der reinen Fliissig- 


1 Vel. Marxerar und Nikurapse [11], die die Ab- 
weichungen des Molvolumens yon der Mischungsregel 
auf das Zusammenwirken eines rein geometrischen 
Packungseffektes mit zwischenmolekularen Kraftwir- 
kungen zuriickfihren. 

2 In einzelnen Fallen versuchten Parsnap [12], 
Sacuer [10] u. a. den Verlauf der Schallkennlinie mit 
Hilfe von Assoziationen, Dipolwirkungen u. 4. zu deu- 
ten. 
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als)5, 


keiten in bezug auf GréBe und Form voneinander 
unterscheiden. Da man annehmen kann, da sich 
das Molekilvolumen B der reinen Substanzen in der 
Mischung nicht sehr andert, bestimmt im wesent- 
lichen der StoBfaktor s den Verlauf des Volumens 
sB tuber der Konzentration. 

In den Bildern 1a bis 14a sind tiber der Konzen- 
tration 


WV =Ve Va 
A(s B) = (sB),— (s B) m 


und 


aufgetragen. Dabei ist Vg das gemessene Molvolu- 
men, V,, bzw. (sB), sind die der Mischungsregel 
gehorchenden Werte, (s B), ist das mittels Gl. (1) 
aus der gemessenen Schallgeschwindigkeit berechnete 
Volumen sB. Man sieht, dafS mit Ausnahme der 
Mischung 9 (iso-Propylalkohol — Cyclohexanol, vel. 
Bild 9a) in allen Fallen das Volumen s B in der Mi- 
schung kleiner ist als es der Mischungsregel ent- 
spricht, und zwar vollig unabhangig vom Verlauf 
des Molvolumens?. 


3.2. Berechnung der Schallgeschwindigkeit in bina- 
ren Fliissighkeitsmischungen mittels S bzw. W 


Es erhebt sich die Frage, ob ahnliche halbempiri- 
sche Beziehungen wie sie die Gl. (2) bzw. (3) dar- 
stellen, auch zur Deutung der Konzentrationsabhan- 
gigkeit der Schallgeschwindigkeit binarer Mischun- 
gen herangezogen werden konnen, d. h., ob die zwei 
Komponenten der Mischung gewisse ,,Inkremente“ 
zu irgendeiner geeignet zu definierenden Kenngrofe, 
z. B. zum Schallvolumen S, additiv beisteuern. 

Diese Frage hat O. Nomoro [13] an 37 binaren 
Fliissigkeitsmischungen untersucht* und zwar der- 
art, daf} er aus den gemessenen Dichte- und Schall- 
geschwindigkeitswerten die KenngroBe (in diesem 
Fall die molekulare Schallgeschwindigkeit R) be- 
rechnet, tber der Konzentration auftragt und gra- 
phisch die Linearitat priift. Auf diesem Weg wird 
jedoch die zum Teil sehr beachtliche MeBgenauig- 
keit mancher Arbeiten nicht voll ausgenitzt. 

Eine viel scharfere Priifung einer solchen halb- 
empirischen Arbeitshypothese ergibt sich, wenn man 
umgekehrt verfahrt. Man berechnet zunichst die 
Kenngrofen der reinen Fliissigkeiten aus den be- 
kannten Daten (Dichte und Schallgeschwindigkeit) , 
dann daraus mittels der Mischungsregel den Wert 
fiir eine beliebige Konzentration. Es gilt also fur 
das Schallvolumen S;,) bei der Konzentration x: 


3 Bei der genannten Mischung kénnte die unterhalb 
20 Mol-% Cyclohexanol auftretende Anomalie vielleicht 
auf das Vorhandensein besonderer Molekilkomplexe 
zuriickgefiihrt werden, deren StoBfaktor gréBer ist als 
der des reinen Cyclohexanol. 

4 Er stiitzt sich dabei auf Messungen von Tuomrxkosk1 
und Nurmi, Gasrrevir, Porant, ParsHap u. a. 
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a Ts In. den Bildern 1b bis 14b ist von den 

Soq= Sil =e) +8 es b ahivede Miisehekcr sansa dips Gnae 

tis Be oben erwahnten Fliissigkeitsmischungen die Grof8e 


. ; : | Aus =Uiq)s —Um) aufgetragen, also die Differenz 
abet bed eee =P uno ote peallvolwams ee zwischen der nach GI.(5) berechneten Schallgeschwin- 
reinen Fliissigkeiten. Nun berechnet man aus dieser ‘ae alent ao andider wach der: Niischaseeee ee 
Gleichung die Schallgeschwindigkeit u;,) und ver- A as 38 1 en 
leicht sie mit der gemessenen. Es ist (5) rechneten, “Crone wim) —UnC = #) ia 
8 a : und uy die Schallgeschwindigkeiten der reinen Flis- 
ge ee Sn thco Qa)” Fg (L—x) +S)x}. sigkeiten bedeuten. Die gemessenen Werte der Schall- 
geschwindigkeit, bzw. die Abweichung dieser Werte 


Bilder 1 bis 14. Verlauf der GréBen AV, A(sB), Jus bzw. Aup iiber der Konzentration in Mol-% der an zweiter 
Stelle genannten Komponente. Die Volumina sind in em?, die Schallgeschwindigkeit ist in m/s 
angegeben. Bei den Bildern 1b bis 14b entsprechen die Kurven den berechneten Groen Aus 
bzw. Auy, die eingetragenen Punkte sind Mefwerte. 


AV, A sB) —= 


80 Mol-%1 


80Mol*% 100 


60 


Bild 1. Benzol—t-Butyl- Bild 2. Cyclohexan—Essig- Bild3.Dekalin—n-Propyl- Bild 4. Chlorbenzol —n- 
alkohol. sdure. alkohol. Propylalkohol. 


——J 


cm] (@) av @so wo _consaso 


20 =40 60 — 80Mol% 


0 Molt. 100 


AV, A SB) —= 


0 20 40 60 80Mol-%100 0 20 40-60 80Mol-“/. 100 


©) 20 40 60 80Mol-%100 


- 80 Be 


Bild 5. Heptan—Dichlor- Bild 6. Tetrachlorkohlen- Bild 7. Methylalkohol— __ Bild 8. n-Propylalkohol — 
methan. stoff-Athylalkohol. Cyclohexanol. Cyclohexanol. 
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AV, Als) —~ 


7 


20. 40 60 80Mo-%100 0 20-40 60 BOMot% 00 0 20. 40 60 _—_80Mol-7.100 
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03 
cm3 


02 


20. 40 60  80Mol7.100 
0 


60 80Mol%100 0 20. 40 60 —_80Mol-1,100 


AlsB) 


0 20. 40 60 _80Mal-7.100 


Bild9.iso-Propylalkohol—__ Bild 10. n-Butylalkohol — 


Cyclohexanol. Cyclohexanol. 


0 Mol 100: 


| 


20.40 60 _80Mol%100 0 20. 40 60 _—80Mol, 100 


0 Mol-7. 100 


Bild 13. Glykol—Cyclo- Bild 14. Glyzerin — Cyclo- 
hexanol. hexanol. 


OMol180 


-40 


Bild 11.iso-Butylalkohol— Bild 12. Octylalkohol — 
Cyclohexanol. Cyclohexanol. 


von der Mischungsregel, sind in diesen Bildern eben- 
falls eingetragen. Man sieht, da} mit Ausnahme der 
Mischungen 4, 5, 6, 13 und 14 die gemessenen 
Werte mit den berechneten innerhalb der Mef- 
genauigkeit tibereinstimmen (siehe Bilder 4b, 5b, 
6b, 13b und 14.b). 

Wendet man die vorstehenden Uberlegungen auch 
auf die molekulare Kompressibilitat W (siehe Gl. 
(4)) an, dann erhalt man 

W a) =W,1—2) + Wer =u Mo 


*Ia 


Q(z) 
und daraus analog zu Gl. (5) (6) 
0 1 0%) ¢ 
Ua)w = Aa Wa = hs [Wi ay, xl". 


In den Bildern 4b bis 6b ist auBer der Groke 
Aus = uU;,)5 —Uy,) auch die Differenz Aujy = ucy yw — 
— Un) aufgetragen, bei den Mischungen 13 und 14 
(Bilder 13b und 14b) fallen die Aww- und Aus- 
Kurven innerhalb der Zeichengenauigkeit zusam- 
men. Man sieht, daf} bei den Mischungen 4 und 5, 
bei denen die nach Gl. (5) berechnete Schallgeschwin- 
digkeit von der gemessenen abweicht, die nach 
Gl. (6) berechnete GroBe u;,)jy innerhalb der Mef- 
genauigkeit mit der gemessenen Schallgeschwindig- 
keit tbereinstimmt. Das bedeutet, daB bei diesen 
Mischungen nicht das Schallvolumen S, sondern die 
GroBe W linear mit der Konzentration verlauft. 


158 


Dies gilt jedoch nicht fur die Mischungen 6, 13 und 
14, bei denen auch die nach Gl. (6) berechnete 
Schallgeschwindigkeit uz) nicht mit der gemesse- 


nen ubereinstimmt (siehe Bilder 6b, 13b und 14b). ' 


O. Nomoto [13] kommt zu dem SchluB, daB bei 
der Mischung Heptan — Dichlormethan (in der vor- 
liegenden Arbeit Mischung 5) sowohl R als auch V 
linear mit der Konzentration verlaufen, bei der Mi- 
schung Tetrachlorkohlenstoff — Athylalkohol (Mi- 
schung 6) dagegen nur W. Das liegt, wie schon er- 
wahnt, daran, da8 Nomoro die Linearitat graphisch 
pruft. Geht man dagegen wie vorstehend beschrie- 
ben vor, dann sieht man aus den Bildern 5b und 
6b, daB bei der Mischung 13 nur W, bei der Mi- 
schung 14 dagegen weder S (bzw. R) noch W der 
Mischungsregel gehorchen. 

Bei den Mischungen, bei denen Gl. (5) bzw. (6) 
nicht erfullt ist, fallt auf, da eine Komponente der 
Mischung entweder schwere Atome (z. B. Chlor) ent- 
halt®, oder aber ihre Schallgeschwindigkeit extrem 
hoch liegt (Glykol: wu = 1643 m/s; Glyzerin: 
wu = 1911 m/s). Davon unberihrt bleibt jedoch 
das im vorigen Abschnitt dargestellte Verhalten des 
Volumens sB in der Mischung, d.h., auch bei den 
Mischungen, in denen das Schallvolumen S bzw. die 
molekulare Kompressibilitat 7 nicht der Mischungs- 
regel gehorcht, wird das Volumen sB in der Mi- 
schung kleiner. 


3.3. Die Schallgeschwindigkeit als Funktion der 
Temperatur bei schalldispergierenden Fliis- 
sigkeiten 


Die Schallgeschwindigkeit reiner organischer Fliis- 
sigkeiten und auch binadrer Mischungen nimmt linear 
mit der Temperatur ab, falls bei der MeBfrequenz, 
bei der die Schallgeschwindigkeit gemessen wird, 
keine Schalldispersion auftritt. Ist letzteres jedoch 
der Fall (z.B. bei reiner Essigsaure bei einer Mef- 
frequenz von 2 MHz), so besitzt die Gerade, die den 
Verlauf der Schallgeschwindigkeit uber der Tem- 
peratur beschreibt, einen Knick, wie folgende Uber- 
legung zeigt. 

In Bild 15 sind die Dispersionsstufen fiir mehrere 
Temperaturen eingezeichnet. Mift man nun die 
Schallgeschwindigkeit im Dispersionsgebiet mit kon- 
stanter Frequenz bei steigenden Temperaturen T, 
bis T's, so erhalt man die MeSpunkte 1 bis 5, oder, 
wenn man die Schallgeschwindigkeit iber der Tem- 
peratur auftragt (siehe Bild 16), die doppelt ge- 
knickte Gerade J — II —III. Da mit steigender Tem- 
peratur die Dispersionsstufen groer werden, wer- 
den die Abstande der Parallelen oberhalb der Dis- 
persionsstufen kleiner, d.h. die Geraden I und III 
in Bild 16 divergieren nach hoheren Temperaturen. 


5 Hierauf hat Nomoto [13] bereits hingewiesen. 


R. GENAHR: SCHALLGESCHWINDIGKEIT IN FLUSSIGKEITSMISCHUNGEN 


Vol. 8 (1958) 


1 

{ 13 Ts, 
3 
| 
4 
5 


| 
| 
sale 
ig f— 
Bild 15. Verlauf der Schallgeschwindigkeit tiber der 
Frequenz im Dispersionsgebiet. 
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Bild 16. Die Schallgeschwindigkeit von Fliissigkeiten in 
Abhangigkeit von der Temperatur bei Schall- 


dispersion. 


Die Mischung 2 (Cyclohexan — Essigsaure) zeigt — 


auffallend starke Schallabsorption [14], aus der 
man auf eine erhebliche Dispersion bei 2 MHz 
schlieBen kann. Deshalb wurde vom Verfasser in 
dieser Mischung die Schallgeschwindigkeit tempe- 
raturabhangig gemessen. Dabei wurde der zweite 


»Knick* (Ubergang von I nach II in Bild 16) tat-— 
sachlich festgestellt, wie Bild 17 zeigt. Da mit der | 
Schalldispersion zusammen auch eine erhohte Schall- | 


absorption auftritt, die mit steigender Essigsaure- 
konzentration und steigender Temperatur stark zu- 
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Bild 17. Verlauf der Schallgeschwindigkeit iiber der 
Temperatur in der Mischung Cyclohexan— 
Essigsiure. Kurvea: 4,68 Mol-%, Kurveb: 
23,6 Mol-%, Kurve c: 56,3 Mol-% Essigsaure; 
MeBfrequenz f=2 MHz. 
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nimmt, war es mit dem zur Verfiigung stehenden 
Ultraschall-Interferometer nicht moglich, auch bei 
héheren Essigsdurekonzentrationen bzw. bei reiner 
Essigsaure diesen Knick zu messen. 

Lamp, AnpREAE und Birp [15] haben bei reiner 
Essigsdure die Schallgeschwindigkeit bei 0,5, 4,0 
und 5,0 MHz temperaturabhangig gemessen (siehe 
Fig. 3 in [15]). Da ihre Meffrequenzen unter- bzw. 
oberhalb der Dispersionsfrequenz liegen, haben sie 
innerhalb ihres Temperaturbereiches nur die Gera- 
den I und III (vgl. Bild 16) gemessen, nicht aber 
die Verbindungsstiicke II. Man sieht aber deutlich, 
daf} diese Geraden mit steigender Temperatur diver- 
gieren. 


Herrn Professor Dr. H. O. Kyeser mochte ich fiir 
die Anregung zu dieser Arbeit und seine Unter- 
stiitzung besonders danken. Der Deutschen For- 
schungsgemeinschaft verdanken wir einige der be- 


nutzten Apparate. (Eingegangen am 3. August 1957.) 
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ZWECKMASSIGE EINTEILUNG DER ELEKTROMECHANISCHEN 
WANDLER MIT RUCKSICHT AUF IHRE BERECHNUNG 
AUS IHREN ELEKTRISCHEN UND MECHANISCHEN DATEN 


von A. Lenk 


Mitteilung aus dem Institut fiir Elektro- und Bauakustik der Technischen Hochschule Dresden 


Zusammenfassung 


In dieser Arbeit wird eine Systematik der elektromechanischen Wandler angegeben, die 
sich in verschiedenen Punkten von den in der Literatur vorhandenen Einteilungen unter- 
scheidet. 

Fiir die vier Gruppen, in die sich die elektromechanischen Wandler unterteilen lassen, 
wird ein sehr einfaches Einteilungskriterium angegeben, das sich nur auf einfache konstruk- 
tive Eigenschaften des Wandlers bezieht. Jeder Gruppe ist dann ein Ersatzschaltbild zuge- 
ordnet, dessen Elemente sich aus dem Aufbau des Wandlers einfach bestimmen lassen. 


Summary 

In this paper a systematic is given for the electro-mechanical transducer, which differs 
in several points from that commonly given in the literature. 

A very simple criterion for allotting such a transducer to one of four groups is given 
which depends only on simple properties of the construction. Each group is then classified 
into a substitute circuit, the elements of which can be derived from the construction of 
the transducer. , 


Sommaire 

On donne dans le présent travail un classement systématique des transformateurs électro- 
mécaniques; ce classement différe 4 certains égards de ceux déja publiés ailleurs. 

On indique un critére trés simple de classement dans les 4 groupes proposés pour les 
transformateurs électromécaniques, ce critére ne faisant appel qu’a des caractéristiques 
de construction simples. On établit ensuite pour chaque groupe un schéma équivalent dont 
les éléments constitutifs peuvent étre déterminés de fagon simple a partir de la constitution 


du transformateur. 


1. Einleitung 


Die Frage der systematischen Einteilung der elek- 
tromechanischen Wandler wurde zuletzt in einer Ar- 


beit von F. A. Fiscuer [1] angeschnitten. Einige 
Punkte dieser Arbeit waren der AnlafB zu einer Dis- 
kussion des Verfassers mit Herrn Dr. F. A. Fiscuer, 


160 


in deren Verlauf ein grofer Teil der hier wieder- 
gegebenen Ergebnisse geklart und formuliert wurde. 
Es besteht im tbrigen tiber alle Punkte der hier ver- 


éffentlichten Arbeit vollstandige Einstimmigkeit des - 


Verfassers mit Herrn Dr. F. A. Fiscuzr. 

Zur Beschreibung von elektromechanischen Wand- 
lern werden weitgehend elektrische Ersatzschaltbil- 
der oder unter Verwendung der elektromechanischen 
Analogien elektromechanische Vierpolersatzschalt- 
bilder verwendet [2], [3], [4]. 

Mit Riicksicht auf eine rationelle Beschreibung 
entsteht dabei die Frage, welche Gruppen von Wand- 
lern jeweils durch ein gemeinsames Ersatzschaltbild 
dargestellt werden konnen, und dariber hinaus die 
weitere Frage, ob die so erhaltenen einzelnen 
Wandlerklassen durch irgendwelche gemeinsame 
konstruktive Eigenschaften oder gemeinsame physi- 
kalische Merkmale der Energieumsetzung gekenn- 
zeichnet sind. 

Solche Einteilungen sind in [2], [5], [6] ent- 
halten. Sie enthalten gemeinsam das Ergebnis, dah 
nach der Art des Ersatzschaltbildes vier verschiedene 
Klassen von Wandlern unterschieden werden miis- 
sen. Beziiglich der Frage, welche Wandler zu welcher 
Klasse gehoren und welche gemeinsamen konstruk- 
tiven oder physikalischen Eigenschaften eine Klasse 
kennzeichnen, unterscheiden sich die obengenannten 
Arbeiten in einigen Punkten. 

In der vorliegenden Arbeit sollen diese Unter- 
schiede geklart werden, und weiterhin wird ein 
neues, sehr einfaches Kriterium angegeben, nach 
dem die Zuordnung eines Wandlers zu einer der 
erwahnten vier Klassen ohne weiteres moglich ist. 


2. Ersatzschaltbilder der elektromechanischen 
Wandler 


Aus [2], [6], [7] kann man entnehmen, dal 
man im quasistatischen Zustand, der vorlaufig be- 
trachtet werden soll, samtliche elektrischen Wandler 
durch die beiden miteinander identischen Ersatz- 
schaltbilder nach Bildla und 1b, und samtliche 
magnetischen Wandler durch die beiden ebenfalls 
miteinander identischen Schaltbilder nach Bild 2a 
und 2b darstellen kann. Unter elektrischen Wand- 
lern sind solche zu verstehen, bei denen die Energie- 
umwandlung mit Hilfe von elektrischen Feldern vor 
sich geht. Analog dazu geschieht bei magnetischen 
Wandlern die Energieumwandlung mit magnetischen 
Feldern. 

In den Bildern 1 und 2 haben C; und L; die Be- 
deutung der Kapazitat bzw. Induktivitat im mecha- 
nisch festgebremsten Zustand. n, und n, sind die 
mechanischen Nachgiebigkeiten fiir konstante elek- 
trische Feldstarke (elektrischer Kurzschlu8) und kon- 
stante dielektrische Verschiebung (elektrischer Leer- 
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Bild 2. Ersatzschaltbild der magnetischen Wandler, 
(a) mit Reiheninduktivitat, (b) mit Parallel- 
induktivitat. 


lauf). Entsprechend bedeuten im magnetischen Fall 
n die Nachgiebigkeit bei konstanter magnetischer 
Feldstarke (elektrischer Leerlauf) und n, die Nach- 


giebigkeit bei konstanter Induktion (elektrischer — 


KurzschluB). & ist der elektromechanische Kopp- 
lungsfaktor. 

Auf der mechanischen Seite wird die Analogie 
u—>v,i—F verwendet. 


3. Unterteilung der Wandler in vier Klassen 


Interessiert man sich nur fiir Messungen an be- 
reits fertigen Wandlern, so waren damit die aufge- 
stellten Fragestellungen schon beantwortet. Es gibt 
danach zwei verschiedene Sorten von Wandlern, 
namlich elektrische und magnetische, die durch die 
Ersatzschaltbilder nach den Bildern 1 und 2 darge- 
stellt werden konnen. Die darin enthaltenen Kon- 
stanten C;, L}, X, Y, nj, ny konnen samtlich durch 
geeignete Messungen an den Klemmenpaaren des 
Wandlers bestimmt werden. Es besteht also bis jetzt 
kein Anlaf fiir eine weitergehende Unterteilung. 

Das andert sich jedoch, wenn man sich die Auf- 
gabe stellt, die Eigenschaften eines Wandlers aus 
seinem konstruktiven Aufbau zu bestimmen. Die 
Verhaltnisse sollen zunachst einmal am Beispiel der 
elektrischen Wandler betrachtet werden. 

Zu dieser Berechnung mu man den Zusammen- 
hang der GroBen C;, X bzw. Y und n, bzw. ny, mit 
den Konstruktionselementen kennen. Die Kapazitat 
C; laBt sich aus dem Aufbau des Wandlers in jedem 
Fall einfach bestimmen, ohne daB der Mechanismus 
der Energieumwandlung naher bekannt sein muf. 
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Es handelt sich dabei um ein rein elektrisches Pro- 
blem. Nicht so einfach ist es auf der mechanischen 
Seite. Im System ist eine rein mechanisch bedingte 
Nachgiebigkeit n, etwa in Form von Federn, ge- 
spannten Membranen o. 4. vorhanden, deren Zu- 
sammenhang mit nm und n, erst noch festzustellen 
ist. 

An dieser Stelle ergibt sich nun eine nochmalige 


_ Unterteilung. Es zeigt sich naémlich, da bei einem 


Teil der elektrischen Wandler die rein mechanisch 
bedingte Nachgiebigkeit n, bei elektrischem Kurz- 
schluB und bei einem anderen Teil im elektrischen 
Leerlauf an den mechanischen Klemmen gemessen 
werden kann. Im ersten Fall ist also ngs=n,, und 
im zweiten Fall ist n,=n,. Der Zusammenhang zwi- 
schen n; und n, ist durch den elektromechanischen 
Kopplungsfaktor k bestimmt: m/n,=1—k?. Die 
Kopplungsfaktoren tblicher Wandler liegen zwischen 
10 und 30%, so dafs die beiden Nachgiebigkeiten 
um 1 bis 10% voneinander verschieden sind, was 
oft vernachlassigt werden kann. 

Es zeigt sich nun weiter, daB bei allen Wandlern, 
bei denen die elektrischen FeldgroBen auf den me- 
chanischen Gro8en v und F senkrecht stehen n, = ni 
ist (tangential-dielektrische Wandler, piezoelektri- 
sche Langsschwinger). Bei allen Wandlern, bei denen 
elektrische und mechanische Grofen parallel zuein- 
ander liegen, gilt ns = (normal-dielektrische Wand- 
ler; piezoelektrische Dickenschwinger) . 

Die Richtigkeit dieser Behauptung lat sich bei 
den dielektrischen Wandlern leicht nachpriifen, weil 
dort die Groe von n; aus der mechanischen Kon- 
struktion zu entnehmen ist. Man kann auch dann an 
den mechanischen Klemmen messen, wenn man die 
PolarisationsgroBe abschaltet. 

Bei den piezoelektrischen Wandlern ist das nicht 
so einfach einzusehen. Elektrischer Energiespeicher 
(Kapazitat) und mechanischer Energiespeicher (Ela- 
stizitat) sind hier untrennbar miteinander verkniipft, 
und es hat zunachst offenbar keinen Sinn, von einer 
rein mechanisch bedingten Nachgiebigkeit oder einer 
rein elektrisch bedingten Kapazitat zu sprechen. Die 
Abschaltung der inneren Polarisation ist auch ohne 
gleichzeitige erhebliche Anderung aller anderen 
Eigenschaften nicht méglich, so da nicht festgestellt 
werden kann, welche der beiden Nachgiebigkeiten 
und n,, als die des mechanischen Systems zu betrach- 
ten ist. 

Diese Entscheidung wird méglich, wenn man die 
stetige Massen- und Federungsverteilung langs der 
Schwingungsrichtung beriicksichtigt, also die Be- 
schrankung auf den quasistationaéren Zustand fallen 
lat. Aus den Schaltungen des Bildes 1 werden dann 
die Schaltungen des Bildes 3 a,b. Man erkennt, da 
sich Dickenschwinger und Lingsschwinger verschie- 
den verhalten. Beim Langsschwinger wird aus der 
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Nachgiebigkeit n, der mechanische Mitgang, den 
man nach der mechanischen Anordnung erwarten 
miiBte. Die Wellengeschwindigkeit ¢ ist hier durch 
die elastische Konstante bei fester Feldstarke be- 
stimmt. 
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Bild 3. (a) Ersatzschaltbilder des einseitig eingespann- 
ten piezoelektrischen Dickenschwingers, 
(b) Ersatzschaltbild des einseitig eingespannten 
piezoelektrischen Langsschwingers. 


Beim Dickenschwinger tritt der gleiche Mitgang, 
der hier nur durch die elastische Konstante bei fester 
Verschiebung bestimmt ist, an die Stelle der Nach- 
giebigkeit nj. Der Grund dafiir ist, da} beim Langs- 
schwinger die Feldstarke E langs der Schwingungs- 
richtung konstant ist, wohingegen beim Dicken- 
schwinger dasselbe fiir die Verschiebung D der Fall 
ist. 

Vergleicht man nun die entsprechenden dielektri- 
schen Wandler (Bild 4a, b), bei denen das Dielektri- 
kum gleichzeitig die Rolle des mechanischen Systems 
iibernimmt, so erkennt man, daB der piezoelektri- 
sche Dickenschwinger das vollstandige Analogon zum 
normal-dielektrischen Wandler ist. Dasselbe gilt fiir 
den piezoelektrischen Langsschwinger und den tan- 
gential-dielektrischen Wandler. Der Unterschied be- 
steht lediglich darin, daf man beim dielektrischen 
Wandler unmittelbar erkennt, was unter dem mecha- 
nischen System zu verstehen ist, wohingegen das bei 
den piezoelektrischen Wandlern erst aus dem dyna- 
mischen Verhalten zu sehen ist. 

Bei den magnetischen Wandlern liegen die Ver- 
haltnisse ganz ahnlich. Infolge der Dualitat zwischen 
elektrischen und magnetischen Vorgiéngen vertau- 
schen sich jedoch die Zuordnungen. 

Man mibt die mechanische Nachgiebigkeit im elek- 
trischen Kurzschlu8 (n,;=n,), wenn magnetische 
und mechanische Gro8en parallel zueinander sind 
(normal-elektromagnetischer Wandler). Stehen ma- 
gnetische und mechanische GroBen senkrecht aufein- 
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Bild 4. (a) Ersatzschaltbilder des normal-dielektrischen 
Wandlers, 
(b) Ersatzschaltbild des 
schen Wandlers. 
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ander, so miBt man die mechanische Nachgiebigkeit 
ns im elektrischen Leerlauf (n,= 7), tangential-elek- 
tromagnetischer Wandler, elektrodynamischer Wand- 
ler). 

Die gleiche Unterteilung, wie sie eben beschrieben 
wurde, findet man durch die Betrachtung des Kraft- 
gesetzes [2]. Sucht man eine Beziehung zwischen 
der erzeugten Kraft F und einer elektrischen Grofe, 
in der die Schnelle nicht vorkommt, so erhalt man 
sowohl bei den elektrischen als auch bei den magne- 
tischen Wandlern die erwahnte Unterteilung in zwei 
Gruppen. 

Bei der einen Gruppe hangt F auf die beschrie- 
bene Weise von der Spannung w, und bei der ande- 
ren Gruppe von dem Strom i ab. 

Es gelten dabei im einzelnen die Zuordnungen: 


ng =n > Fant 5 


ng=n—>Frnu. 


Der Nachteil dieses Kriteriums liegt darin, da 
man den physikalischen Mechanismus der Energie- 
umsetzung schon kennen mu, was bei dem zuerst 
angefihrten Kriterium (Richtung der Verkopplung) 
nicht notig ist. 

In der Literatur [2], [5] ist noch eine Reihe 
anderer Kriterien angegeben, die ebenfalls eine 
Unterteilung sowohl der elektrischen als auch der 
magnetischen Wandler in je zwei Gruppen liefern. 
Diese Unterteilungen stimmen jedoch nicht mit der 
oben behandelten, aus der Art der Ersatzschaltbilder 
folgenden Unterteilung tiberein. 

Es handelt sich dabei erstens um die Frage nach 
der wesentlichen Beteiligung oder Nichtbeteiligung 
des elektrischen :Energietragers an der Energie- 
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umwandlung [5] und zweitens um die Frage, ob 
lineare oder quadratische Kraftgesetze vorliegen. Im 
ersten Fall ergibt sich, da8 bei allen dielektrischen 
und elektromagnetischen Wandlern der Energietra- 
ger an der Energieumwandlung wesentlich beteiligt 
ist, wohingegen das bei den piezoelektrischen und 
elektrodynamischen Wandlern nicht der Fall ist. Im 
zweiten Fall, bei der Frage nach dem Grad des 
Kraftgesetzes, kann man keine so einfache Aussage 
machen. Obwohl bei den dielektrischen und elektro- 
magnetischen Wandlern immer quadratische Kraft- 
gesetze vorliegen, gibt es piezoelektrische Wandler 
sowohl mit linearem als auch mit quadratischem 
Kraftgesetz. 


4. Zusammenfassung 


Bei der LEinteilung der  elektromechanischen 
Wandler nach der Art ihres Ersatzschaltbildes erge- 
ben sich vier Klassen. In zwei Klassen sind alle 
magnetischen Wandler und in zwei anderen Klassen 
alle elektrischen Wandler enthalten. 

Die Unterteilung der elektrischen und magneti- 
schen Wandler unter sich in zwei Klassen hat nur 
dann eine Bedeutung, wenn man nach dem Zusam- 
menhang der im Ersatzschaltbild vorhandenen Nach- 
giebigkeiten n, (elektrischer Kurzschlu8) und m 
(elektrischer Leerlauf) mit der rein mechanisch be- 
dingten Nachgiebigkeit n,, die auch das dynamische 
Verhalten wesentlich bestimmt, fragt. Diese Frage- 
stellung tritt auf, wenn ein Wandler berechnet wird, 
wohingegen sie bei der Messung eines vorhandenen 
Wandlers ohne Bedeutung ist. 

Fur die erwahnte Unterteilung kann man eine ein- 
fache Regel angeben, die darauf Bezug nimmt, ob 
die mechanischen Groen v, F auf den elektrischen 
bzw. magnetischen FeldgroSen senkrecht stehen oder 
mit ihnen parallel laufen. 


Es gilt das Schema entsprechend Bild 5. 


elektrische Wandler 
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Bild 5. Richtungskriterium fiir die Einteilung von elek- 
tromechanischen Wandlern. 


ns = MW 


Da der Unterschied zwischen n, und n, bei tb- 
lichen Wandlern zwischen 1 und 10% liegt, kann 
er im allgemeinen, wenn es sich nicht um die Be- 
rechnung abgestimmter Wandler handelt, vernach- 
lassigt werden, und man hat lediglich zwischen elek- 
trischen und magnetischen Wandlern mit den Er- 
satzschaltbildern nach den Bildern la und 2a zu 
unterscheiden. 
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Bild 6. Zusammenfassung. 


Fur ihr forderndes Interesse an dieser Arbeit und 
viele wertvolle Ratschlage bin ich den Herren Dr. 
F. A. Fiscuer und Prof. Dr. W. Rercuarpr zu Dank 
verpflichtet. 


(Eingegangen am 1. April 1957.) 
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Der Verfasser behandelt einleitend die Notwendigkeit 
der Sprachaudiometrie und schildert dann die Gesichts- 
punkte, die zur Aufstellung seiner Sprachpriifteste ge- 
fiihrt haben. Diese Teste kénnen durch Umgruppierung 
der Worte sowohl in sinnvoller als auch in sinnloser 
Form yerwendet werden. Anschliefiend an die Beschrei- 
bung der MeSapparatur und der Hoérpriifréume folgt 
eine Diskussion der Methodik der Sprachhorpriifung 
und der Ergebnisse von Untersuchungen des Verfassers. 
Die verschiedenen Einfliisse auf die Lautverstindlich- 
keit werden analysiert; die Kenntnis dieser Zusammen- 
hange kann mitunter dazu dienen, Simulanten zu ent- 
larven. 

Eine neue Definition des Sprachhorverlustes und der 
dadurch bedingten Erwerbsminderung wird erlautert 
und den Definitionen yon Hantsrock, Mrisrer u. a. ge- 
geniibergestellt. Daran schlieBen sich kritische Bemer- 
kungen iiber Priifmethoden und Sprachteste anderer 
Autoren. 


Besonders eingehend wird das _hallaudiometrische 
Verfahren des Verfassers dargestellt, bei welchem ver- 
schiedene Mischungsverhaltnisse von Direktschall zu 
Indirektschall eingestellt werden kénnen. Die Auswer- 
tung dieser verhaltnismafig schlecht zu verstehenden 
verhallten Teste ermdglicht eine Aussage iiber die Fa- 
higkeit, aus einem Frequenzgemisch Sprache heraus- 
héren zu kénnen, und damit iiber den vermutlichen 
Nutzen einer Horhilfe. AbschlieBend wird u.a die An- 
passung von Horgeraten behandelt. 

Der Berichter méchte anregen, bei einer Neuauflage 
des Buches die in der Einleitung aufgefiihrten Grund- 
begriffe und Formeln zu tiberarbeiten und auch kon- 
sequent im Text anzuwenden; auferdem sollte eine 
Reihe von Druckfehlern berichtigt werden. Das Verwen- 
den verschiedener Symbole (z.B. grofe und kleine 
Buchstaben) fiir die gleiche physikalische GroBe er- 
schwert unnétig das Verstaéndnis der Formeln. Die etwas 
uneinheitliche Darstellung ist vermutlich darauf zuriick- 
zuftihren, daB das Buch — wie im Vorwort mitgeteilt 
wird — aus einzelnen Teilen eines langjahrig gelesenen 
audiometrischen Kollegs des Verfassers entstanden ist. 

Das Buch mit seinem sehr ausfiihrlichen Literatur- 
yerzeichnis gibt einen interessanten Uberblick iiber den 
gegenwartigen Stand der Sprachaudiometrie und ihrer 
Priifmethoden. H. G. Diesrer 


FORCED BENDING AND EXTENSIONAL VIBRATIONS OF A 
TWO-LAYER COMPOUND LINEAR VISCOELASTIC BEAM 


by F. Scuwarz 


Centraal Laboratorium TNO, Delft, Netherlands 


Summary 


In the forced vibrations of a two-layer viscoelastic beam four basic cases are to be 


distinguished: 


A pure tensile vibration excited by a tensile force and a bending moment. 
A pure bending vibration excited by a bending moment and a tensile force. 
A superposition of a bending vibration and a longitudinal vibration under a pure bend- 


ing moment. 


A superposition of a longitudinal vibration and a bending vibration under a pure tensile 
force. 

The apparent stiffness against bending or tension, the damping and the energy dissipa- 
tion in these four cases are discussed as functions of the elastic constants and the ratio of 
the layers. 


Sommaire 


Il y a lieu de distinguer quatre cas fondamentaux de vibrations forcées d’une poutre 
viscoélastique a deux couches: 

Une vibration de traction pure, excitée par une force de traction et un moment de 
flexion. 

Une vibration de flexion pure, excitée par un moment de flexion et une force de traction. 

La superposition d’une vibration de flexion et d’une vibration longitudinale sous l’effet 
d’un moment de flexion pur. 

La superposition d’une vibration longitudinale et d’une vibration de flexion sous leffet 
d’une force de traction pure. 

La rigidite apparente a la flexion ou 4 la traction, l’amortissement et la dissipation de 
l’énergie, sont étudiés dans ces quatre cas en fonction des constantes élastiques et du rap- 
port des épaisseurs des couches. 


Zusammenfassung 


Bei den erzwungenen Schwingungen eines aus zwei Schichten zusammengesetzten visco- 
elastischen Stabes sind vier grundlegende Falle zu unterscheiden: 1. Reine Dehnschwin- 
gung, erregt durch Dehnungskraft und Biegemoment. 2. Reine Biegeschwingung, erregt 
durch Biegemoment und Dehnungskraft. 3. Uberlagerung von Biege- und longitudinaler 
Schwingung durch reines Biegemoment. 4. Uberlagerung von Biege- und longitudinaler 
Schwingung durch reine Dehnungskraft. 

Die Steife gegen Biegung oder Dehnung sowie die Dampfung und die Energieverluste 
werden fiir die vier Falle in Abhangigkeit von den elastischen Konstanten und dem Ver- 


haltnis der Schichtdicken diskutiert. 


1. Introduction, Discussion of the concepts of 
“pure bending” and “neutral fibre”, 
Relation to previous work 


It is well known that the vibrational properties of 
linear viscoelastic materials may be described by 
means of a complex elasticity modulus £* = E’ +i E” 
which depends on the frequency ~ of the vibration. 
The storage modulus E’ and the loss modulus E£” 
represent the components of tensile stress in phase 
and out of phase with unit tensile strain. The abso- 
lute modulus E=/VE?+£E? and the loss tangent 
tan d= E/E’ give the ratio and the phase lag be- 
tween stress and strain amplitudes. 

The purpose of this investigation is a study of the 
vibrational behaviour of a rectangular beam consist- 


ing of two layers of different viscoelastic materials 
with the complex elasticity moduli £,* and £,*. The 
breadth of the beam is taken equal to unity and all 
forces and moments will be given per unit breadth. 
A discussion of bending vibrations only will turn 
out to be too restricted. Therefore the superposition 
of bending vibrations and longitudinal tensile vibra- 
tions will be considered. 

The definitions of the physical concepts of “pure 
bending” and “neutral fibre” will meet with un- 
expected difficulties in our problem. Therefore we 
have to start with a careful investigation of the defi- 
nitions concerning pure bending, 

Up to now it has been usual to define pure bend- 
ing of a rectangular beam as a state of strain which 
is obtained under the action of a pure bending 
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moment. In addition the principle of St. Venant? is 
assumed to be valid: The amplitude of the longitu- 
dinal strain is linearly distributed over the cross- 
section” 


&2=a(x-+1r). (1a) 


The tensile strains in the transverse directions, ¢;, 
and é,, are assumed to be equal to — v¢,, (v Pois- 
_son’s ratio) and the shear strains are assumed to 
vanish. This ensures that only one component of the 
stress tensor, the longitudinal stress, o,,, is different 
from zero. There exists one fibre x= —r (cf. eq. 
(1 a)), called the neutral fibre which remains un- 
_deformed and therefore unstressed. The position of 
this fibre, i. e. the value of r, is fixed by the condi- 
tion that only a bending moment but no resulting 
tensile force is acting on the bar. Consequently the 
mean tensile force of the stress distribution 0;,, 
must be equal to zero: 


F= [ o,,dx=0. (le) 


This equation determines r. The bending moment is 
calculated as the moment of the stress distribution 
relative to the neutral fibre 


M= f[ («+r) Ox, dx. 


The conditions (la), (1b) and (lc) together 
form the usual definition of pure bending. This defi- 
nition is a mixed one, partly consisting of condi- 
tions for the strain (eq. (1a)) and partly of con- 
ditions for the stress (eq. (1b) and (1c)). One 
must be prepared that in certain cases the conditions 
(1a) on the one hand and (1b) and (1c) on the 
other hand turn out to be no longer compatible. 
Anticipating the results of the following sections 
we mention three cases of increasing degree of in- 
tricacy. In two of these cases the definition (1) for 
pure bending is admissible, in the third case it fails. 
1. A homogeneous rectangular bar of a linear 

elastic or viscoelastic material is subjected to the 

action of a periodic bending moment after (1 b) 

and (lc). The bar undergoes periodic pure 

bending of the same frequency in the sense of 
(1a). There exists a neutral fibre in the middle 
of the cross-section, which remains unstrained 
during the vibration. 

2. Consider a compound beam consisting of two dif- 
ferent linear viscoelastic materials. The elasticity 


1 A discussion of the principle of Sr. Venanr may 
be found in N. I. Muskxueuisnvi1, Some basic problems 
of the mathematical theory of elasticity. Noordhoff Ltd, 
Holland 1953, pp. 76 and 561. 

2 The position of the coordinate system may be 
found in Fig. 1. The z and x coordinates are directed 
parallel to the longitudinal fibres of the bar and the 
cross-section of the bar in the direction of bending 
respectively. 


(1b) 
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moduli of both layers may be different, their, 
damping, however, is assumed to be equal, 
tan d,=tan0,. Under the action of a periodic 
bending moment, again a certain fibre remains 
undeformed during vibration. This neutral fibre 
is not situated in the middle of the beam. Its 
position is determined by the ratio of the thick- 
nesses and the elasticity moduli of both layers. 

Again eq. (1) are compatible and may be used 

for the definition of bending. 

3. A compound bar with layers of different damp- 
ing is considered, tand,#tan0d,. It will be 
shown that a periodic bending moment without 
longitudinal force results in a vibration which is 
not of the simple form (la). No single fibre 
exists which remains unstrained during the entire 
period of vibration, At each particular instant a 
certain fibre of the beam is unstrained, but this 
will be a different one every moment. The neu- 
tral fibre “migrates” twice through the cross- 
section of the beam during one period of vibra- 
tion. It even leaves the beam during a certain 
part of the period. It will be quite clear that in 
this case the conditions (la) and (1b,c) ap- 
pear to be incompatible. Consequently a redefi- 
nition of the concept of pure bending is neces- 
sary. 

This definition may be based on either (1b, 
c) dropping condition (la) or on (la) drop- 
ping (1b,c). As pointed out above, the first 
possibility results in a vibration without a fixed 
position of the neutral fibre. This would not be 
called a pure bending by an experimental physi- 
cist. We call it a “vibration under pure bending 
moment”; it will be shown to be composed of a 
pure bending and a pure tensile vibration. 

The second alternative is chosen in defining 
pure bending by means of eq. (1a). However, 
difficulties arise now in connection with the 
meaning of the concept of neutral fibre x= —r. 
No preferential position for this fibre exists as 
no value for r can be found, for which the lon- 
gitudinal force F vanishes. (It is impossible to 
excite a pure bending vibration by means of a 
bending moment without a longitudinal force). 
Therefore the meaning of the neutral fibre has to 
be redefined. 

The experimental physicist measuring the deflec- 
tion and the elongation of the beam attaches a poin- 
ter to the free end of it at a fibre a= —r, say. The 
deflection and extension of the beam is then measur- 
ed as the vertical and horizontal movement of the 
pointer, i.e. as the deflection and extension of the 
(arbitrarily chosen) fibre x= —r. This fibre will be 
called the neutral fibre. Pure bending of the beam 
will be defined as a state of strain where the neutral 
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fibre does not undergo an extension, pure elongation 
as a state of strain where the neutral fibre does not 
undergo a deflection. 


If defined in this way, superposed bending and ° 


tensile vibrations of the compound beam may be 

excited by the common action of a periodic bending 

moment and a periodic tensile force. However, a 

bending moment only does not suffice for a pure 

bending vibration and a tensile force only does not 

suffice for a pure tensile vibration. We may distin- 

guish four basic experiments: 

1. A pure tensile vibration excited by a tensile force 
and a bending moment. 

. A pure bending vibration excited by a bending 
moment and a tensile force. 

. A superposition of a bending vibration and a lon- 
gitudinal vibration under a pure bending moment. 

4, A superposition of a longitudinal vibration and a 

bending vibration under a pure tensile force. 


i) 


ow 
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The calculations made by van Oorr and by 
Oxserst may be classified as special cases of our 
problem 2, as mentioned above. Their calculations 
contain certain approximations which will be dis- 
cussed in section 5. 


2. Forced vibrations of a two-layer beam under 
the common action of a periodic bending mo- 
ment and a periodic tensile force 


Forced vibrations of a two-layer compound beam 
are considered. The beam has a length 2/ and a 
width unity. It is composed of two different mate- 
rials, the viscoelastic constants of which are sum- 
marized in Table J. Both materials are assumed to 
be joined rigidly together. 

In the undeformed state the origin of the coordi- 
nate system is placed in the middle of the plane se- 
parating the materials with the x direction perpendi- 
cular to this plane; the z direction is taken parallel 


Table I. Material constants of the two-layer beam 


| Absolute | 


* Storage Thickness Poisson’s 
2 : ‘ dul ‘ 
eee le Damir’ modulus hese moot of layer ratio* 
’ E | | | 
Material 1 | By} tan 0, By’ = B,cos 6, | B,/’ = Hy sin 6, d, | » 
| | | 
Material 2 | E, | tan d, EL,’ = HE, cos 6, E,’’ = EH, sin 6, d, v 
Material 2 tan 6 é 
Ratio ————— 6 = —— a=H,’/E,’ ab = Be’ |B! ~ |b dad = 
Material 1 | tan 0, 2 [By ais | 2ldy 


The apparent stiffness against bending or tension 
and the energy dissipation in these four experiments 
will be discussed in dependence of the elastic con- 
stants and the thickness of the two layers. 

These considerations may be applied to two prob- 
lems of technical interest. The measurement of the 
viscoelastic constants of materials too soft to sup- 
port their own weight may be established by glue- 
ing them to a thin layer of a stiff material. From the 
viscoelastic properties of the lamellated bar the pro- 
perties of the soft material may be deduced, as pro- 
posed by van Oorr?. 

A more important application is the antivibration 
damping of constructional parts in machines, ships 
or cars. The propagation of noise through the metal 
frame of those constructions is usually prevented by 
a cover of antivibration points or of an other highly 
damping material. The attenuation of transverse 
bending vibrations through the total system (metal 
+layer) will be governed by the aparent loss tan- 
gent of this system. Therefore optimum conditions 
are required in order to obtain high losses. This 
problem has been analyzed by Osersr‘. 


3 W. P. van Oort, Transactions of Instruments and 
Measurement Conference, Stockholm 1952. 

4H. Osersr, Acustica 2 [1952], Beiheft 4, AB 181 
— AB 194. 


to the length of the beam (see Fig. 1). The fibres 
of material 1 extend from r= —d, to x=0, those 
of material 2 from x=0 to x= +d). 


Material 2 


Material 1 


Fig. 1. Right half of the composite two-layer beam and 
position of the coordinate system. 


The deflection and the elongation of the beam 
will be defined by means of the movement of a 
pointer attached to the free end of the beam at a 
fibre x= —r. The distance r between this fibre and 
the boundary between the layers will not be specifi- 
ed for the present. In special cases the value of r 
may either be given by the experimental arrange- 
ment or be chosen by the desire for simple mathe- 
matical expressions. 

We shall try to solve the problem by means of the 
following system of displacements (u, v, w are the 


5 Poisson’s ratio is assumed to be independent of 
the frequency and to be equal for both materials. 
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displacements in the x, y, z direction respectively) : 


W(x, Y, 25-4) = — Fa[2+y(r+r)2—vy?] sin wt— 
—sv(x+r)sin(mt—%), (2) 
v(x, y,2z;t) = —ay y(x+r)sinwt— 


—svysin(wt—%), 
w(x, y, 2; t) =az(e4+r) snwt+szsin(wt—y). 


This system of displacements corresponds to the 
superposition of a pure bending with the amplitude 
a and the neutral fibre at x = —r, and an elongation 
in the longitudinal direction with the amplitude s. 
A phase lag 7 is assumed to exist between bending 
and elongation. 


The system of strains corresponding to (2) is 
found as 


Vey =Vy2= Yer =O, Exz = Eyy = —Y Ezz 5 (3) 


&,=—a(x+r) sinwt+ssin(wt—y7) = 


=a(x+r)sinwt+s sinwmt—s’ coswt. 


The tensile strain, ¢,,, consists of a bending part, 
which varies linearly with the x coordinate ®, and an 
extensional part, independent of x. The total state 
of strain seems to be determined by four constants, 
viz. a, r, s and 7. However, if the extension is sepa- 
rated in a component s’=scosy7 in phase and a 
component s’=ssin7 out of phase with the bend- 
ing vibration, the state of strain is seen to be com- 
pletely determined by three constants, viz. a, s’” and 
ar+s . Up to now the decomposition of (3) in a 
bending and an extension in phase with it is purely 
artificial. 

However, the fourth constant, r, obtains an inde- 
pendent physical significance, if we consider the de- 
flection and the elongation of the beam. The deflec- 
tion of the beam is measured asthe vertical dis- 
placement of the pointer, viz. the displacement of 


the fibre x= —r in the negative x direction (4) 
D= —u(z= —r,y=0,z=l; t) ae Pasinwt. 
Similarly the elongation of the beam is measured as 
the displacement of the fibre «= —r in the positive 
z direction (5) 
L=+w(e¢= —r,y=0,z=l; t) =Issin(wt—7). 


By the definitions (4) and (5) the meaning of r is 
fixed, x= —r being the fibre at which elongation 
and deflection measurements are to be performed. 


® This special form of the bending deformation is a 
consequence of the assumption of Sr. Venanr (vanish- 
ing shear stresses). The equality of the strain gradient 
in both materials is a consequence of the continuity of 
the displacements at the plane x=0. (The materials 
are joined rigidly together.) 
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Next the stress distribution is calculated from 
(3) by means of the equations for linear viscoelastic 
behaviour’. The only stresses different from zero 
turn out to be the tensile stresses in the z direction. 
These show a component in phase with the strain 
proportional to E£,’, and a component out of phase 
with the strain proportional to E,”’ 


Og = Oya Oey =O pn Oy = 0, (6) 
Oz, =a(Ey sinwt+E,” coswt) (x+r) + 
+s[E,’ sin(wt—7) +E,” cos(mt—7)] 
(for material 1). 


For material 2 a corresponding formula is valid 
with all the lower indices changed correspondingly. 


The stress system (6) satisfies the equilibrium 
conditions without inertia terms 


O0,,/Ox + O0z,/dy + 96,,/dz =0 


and therefore constitutes a possible solution of the 
viscoelastic equations®. Further it may easily be 
seen that (6) obeys the right boundary conditions: 
All lateral faces of the beam are free from external 
forces, with the exception of the planes perpendicu- 
lar to the z direction. On the cross-sections of the 
beam two types of resulting forces are found. The 
stresses are equivalent to a longitudinal tension 


F= [ o..dx (7) 


cross-section 


and a bending couple 


M= [ (e+r) o,,dx. (8) 
cross-section 
Here we have chosen x= —r as the fibre relative to 
which the bending moment is defined. This is in ac- 
cordance with the selection of this fibre in the de- 
finitions (4) and (5). 

If the integrations in (7) and (8) are performed, we 
get the tensile force and the bending moment each 
as a sum of four components: 

A component proportional to s and in phase with 
the elongation, 

a component proportional to s and out of phase 
with the elongation, 

a component proportional to a and in phase with 
the deflection, 


7 These equations valid for linear viscoelastic isotro- 
pic materials under sinusoidal strain under the assump- 
tion of a frequency-independent Poisson’s ratio » may 
be found in H. A. Sruart, Die Physik der Hochpoly- 
meren, Band IV, § 6c. 


8 Neglecting the inertia terms in the dynamic equa- 
tions means that our solution is restricted to forced 
nonresonance vibrations; however, we know from the 
theory of elasticity that the expressions for the stiff- 
ness obtained in this way, may also be used for the 
approximate description of resonance vibrations. 
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a component proportional to a and out of phase 
with the deflection; 


F =sR'sin(mt—™%) +s R’ cos(mt—%) + 


+am' snwt+am’ cost, (9) 
M=sm’' sin(wt—7) +sm’” cos(mt—%) + 
+aB’snwt+aB’ cos@t. (10) 


R’ and R” represent the tensile force due to elon- 
gation, B’ and B” represent the bending moment due 
to the deflection. m’ and m” are cross components 
which give the tensile force due to the deflection and 
the bending moment due to the elongation. All these 
characteristic quantities depend on the viscoelastic 
constants of both materials and on the thicknesses 
of both layers. m and B depend further on the 
choice of the fibre 7 where the elongation, the deflec- 
tion and the bending moment are determined. 


Rees? de Ba Rend che de eee (ON) 

MoatER = (kde y de) 

Mm SOR =e (hy dee de (12) 

Bos SE df 2 Ey dsr df hs ye 
+r R’, 

Bee (Edge ky dgy 1 Grd? 8, de) + 
472 RR”, (13) 


Finally we calculate the energy dissipation per 
period; the energy dissipated per unit of volume 
and per period of the vibration is equal to 


By using (3) and (6) this leads to 
w=nE,’[o?(x+r)?2+5?4+2as(x+r) cosy] 


for material 1 


and to a similar formula for material 2. The energy 
dissipation per period in a strip of the beam of unit 
length and unit width is obtained by integration 
over the cross-section 


_ [wdx= (a? BY + s? R” +2asm” cosy).(14) 


cross-section 


3. Discussion of a special case: Equal damping 
in both layers 


The most simple conditions are obtained if the 
damping in both layers is the same, viz. tan 0, = 


=tan6,=tan 0. It follows immediately that 
Reith ei elie} 


However, we may still choose the position of the 


VIBRATIONS 


- (15)) m” vanishes together with m’ and the tensile 
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neutral fibre r. There is only one choice which re- 
sults in m’ =0, namely 


r=ro = (Ey d,?— Ey’ dy?) /(Ey d, + Es dy). (16) 


Because of the equality of tand, and tan 0, (eq. 


force and the bending moment do not show cross- 
terms 
F =sR’ sin(wt—y) +s R’ cos(mt—y), 


M=aB' sinwt+aB’ coswt. (17) 


If this “natural” choice r= ro for the neutral fibre | 
is used, pure bending may be excited by a bending — 


moment only and pure elongation may be excited 
by a tensile force only. This is a fundamental pro- 
perty inherent to the special case tan 0, = tan 0,. 

If tan 0, # tan 0, as generally, there is no choice 
of r, by which tensile force and bending moment 
may be separated. Still we may cause m’ to vanish 
by a proper assumption about r, but m” will remain 
different fiom zero or vice versa. There is no longer 
a “natural” position for the neutral fibre and the 
choice of r becomes an arbitrary decision, either 
prescribed by the experimental conditions or by the 
desire for simple mathematical expressions ®. 

In the following discussion we shall still use 
r=ry as the position of the neutral fibre. By this 
choice m’ vanishes and m” takes a simple form. 


(18) 
ry — Vidi +d) By By de 


see tan 05 — tan 0,). 
2° Ed +E do 


4. Discussion of the four basic experiments 


We return to the general case of unequal damp- 
ing and discuss the four basic experiments, where 
successively the deflection, the elongation, the tensile 
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force and the bending moment are assumed to — 


vanish. Note that the definitions of the experiments 
sub 4.2, and 4.4, depend on the choice for the situa- 
tion of the neutral fibre. 


4.1. Pure extensional vibration; D=0 


Assuming a = 0 in the general equations and using 
(16) and (18), we get immediately the result: A 
pure extensional vibration without deflection, 


D=6 


L=Issinwt, 


may be excited by means of a periodic tensile force 
and a periodic bending moment of the same fre- 


: The choice of r will generally be restricted by the 
requirement 
tan d,/tan 6, . 
But this restriction is not sufficient to determine the 
expression for r in the general case tan d 3 tan 0, . 


limr—r,_ for 


[7 wn 
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quency. The periodic tensile force precedes the ex- 
tension by an angle 4: (19) 


F=sR sinwt+sR’ coswmt=sRsin(wt+ A). 


The storage and loss parts of the extensional stiff- 
ness are given by 


R =E, d,1+a6), R’ =E, d,1+ab&), (20) 


where the abbreviations of Table I have been used. 
The absolute stiffness R and the phase lag 4 are 
determined as usually from the equations 
R2-R? +R, tand=R’/R’. 
Apart from the tensile force, a periodic bending 
moment is required, which precedes the extension 
by 90 degrees. 


M=smsin(mt+n/2). (21) 


Its amplitude is given by eq. (18). 

These formulae are illustrated by a special 
example: Consider a two layer beam of total thick- 
ness unity (d; +d = 1), which consists of a material 
with low damping and high modulus (tan 0, =0.1; 
E,=1), and of a material with high damping and 
low modulus (tand,=10; £,=0.1). This corres- 
ponds to the following ratios 


EJ/E,=0.1, a=£,'/E,/=0.01, 
b=tan 6,/tand;=100, ab=E,”/E,” =1. 


Re 


12 


The energy dissipation per unit of deformation, the 
loss modulus, is the same for both materials. 


10' 


wl 
107 10" 10° 10' 10? 10° 10° 
£=d,/d,—> 


Fig.2. The characteristic functions in pure tensile 
vibration for a two-layer beam according to for- 
mulae (19), (20) and (21): E,=1; tano, 
=0.1; E,=0.1; tand,=10 (R’,R” storage 
and loss part of the extensional stiffness, R ratio 
between the amplitudes of longitudinal force 
and extension, 4 phase lag between longitudi- 
nal force and extension, m ratio between the 
amplitudes of bending moment and extension). 
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Fig. 2 shows the quantities of formulae (19), 
(20) and (21) as a function of the thickness ratio 
£=d,/d,. With increasing relative thickness of the 
damping layer, the absolute stiffness R and its sto- 
rage part R’ decrease monotonously and the phase 
lag between elongation and tensile force increases 
monotonously, The energy dissipation per unit of 
elongation, R’, is independent of the thickness ratio 
as to be expected. The amplitude of the additional 
bending moment, m, shows a maximum at the thick- 


ness ratio =a 1/210. 


4.2. Pure bending vibrations; L=0 


The assumption s=0 leads the following results: 
A pure bending vibration without extension 


D=Pasinwt, L=0 


may be realised by means of a periodic bending 
moment and a periodic tensile force of the same fre- 
quency. The periodic bending moment precedes the 
deflection by an angle 0 


M=aB' sinwt+aB’ coswmt=aBsin(wt+0). 


The storage and loss parts of the bending stiffness 
are given by 


Ey dP 1+2a&(24+36+2 &) +a? 64 


B= ee 
12 (1+a@&) 


5 


(22) 


B?=B?+B2, tand=B’/B’, (23) 


Apart from the bending moment, a periodic ten- 
sile force is required, which precedes the deflection 
by 90 degrees: 


F=amsin(wt+n/2). (24) 


The behaviour in pure bending is illustrated in 
Fig. 3. Again the stiffness is decreased and the phase 
lag between bending moment and deflection is in- 
creased by an increase of the relative thickness of the 
damping layer. However, this course is no longer mo- 
notonous as there is a relative maximum in tan 0 and 
a relative minimum in B at thickness ratios of about 
10 and 1 respectively. At these thickness ratios a rela- 
tive increase of the amount of the soft, damping 
material may lower the damping of the total system 
or raise its bending stiffness. It will be quite clear 
that the situation of these extreme values may be 
of importance for the design of practical systems. 
The energy dissipation per unit of deflection, B’ 
equals 0.1 for a beam of either material 1 or mate- 
rial 2. However, if both materials are joined rigidly 
together, the energy dissipation in bending may be 
up to three times as high. 


plz 
10° 10" 10° 10' 10° 10° 10" 
6-0/0. 


Fig.3. The characteristic functions in pure bending 
vibration for a two-layer beam according to for- 
mulae (22), (23) and (24); the material con- 
stants are chosen as in Fig.2 (B’,B” storage 
and loss part of the bending stiffness, B ratio 
between the amplitudes of bending moment and 
deflection, 6 phase lag between bending mo- 
ment and deflection, m ratio between the am- 
plitudes of extensional force and deflection). 


4.3, Vibration under pure bending moment; F =0 


A vibration under pure bending moment is cha- 
racterized by the condition F=0. This condition 
yields equations for ssin 7 and s cos 7 

ssin yn =amR’/R?, (25) 
scosy7= —amR”/R?. 
If the deflection is assumed to be a sinusoidal func- 
tion of time, 
D=;P asin wt, 


the longitudinal vibration is determined by means of 


(25) as follows 
S 4) 


Inserting (25) into (10) yields the bending moment 


pS la 7 sin for (26) 


to) a 


M=aB' sinwt+aB" coswt=aBsin(mt+ i 
2 
=aBsin(wt+6) +a, sin(w1— A), (27) 
The storage and loss parts of the flexural stiffness 
now become 

B’ =B’ im? RR, 

Pe Bei? R’/R2, (28) 

B? =B? +B’, tand=B’/B’. (29) 

The result of a periodic bending moment will 

therefore be a vibration composed of a bending and 

an extension of the same frequency. The phase lag 

between bending and extension equals 4 +2/2 and 

the ratio of their amplitudes equals m/R (determin- 


ed by eq. (18) and (20)). 
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As may be seen from (14) and (25) the energy 
dissipation is given by W =xo? B’. Eq. (28) con- 
stitutes a direct comparison between the flexural 
stiffness in pure bending and the flexural stiffness 
under pure bending moment. The storage part of the 
flexural stiffness is always larger in the case of bend- 
ing under pure bending moment, whilst the loss part 
(and therefore the energy dissipation) as well as the 
loss tangent are always larger in pure bending. 

The behaviour in vibration under a pure bending 
moment is illustrated in Fig. 4. 


10° 


= Ss 
cea Ba 
oS 


0.1(m/R) 


10? 10" 40° 10' 10° 10° 10" 
$= d2/d, —~ 


Fig. 4. The characteristic functions in vibration under 
pure bending moment for a two-layer beam 
according to formulae (26), (28) and (29); 
the material constants are chosen as in Fig. 2 
(B’, B” storage and loss parts of the bending 
stiffness, B ratio between amplitudes of bend- 


ing moment and deflection, 0 phase lag between » 


bending moment and deflection, m/R_ ratio 
between the amplitudes of extension and de- 
flection) . 


4.4. Vibration under pure tensile force; M=0 


This case is characterized by the condition M=0. 
Inserting this into (10), we get two equations which 
yield 

Dp swe ns ee (30) 

acosy= —smB"/B?. 
If the extension is assumed to be a sinusoidal func- 
tion of time 


L=Issinwt 


the deflection is obtained by means of (30) as 


Damp PoE sin Ge z -8). (31) 


rs 


The tensile force is found by inserting (20) into 
(9) 
F=sF sinwt+sF’ cosw =sFsin(wt+A) = 
2 
=sRsin(wt+A)+s 7 


sin(wt—6). (32) 


¥ 
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Table II. Description of the four basic experiments in terms of the characteristic quantities R, 4, B, 0, m 


Pure tensile Pure bending Vibration under Vibration under 
vibration vibration pure bending moment pure tensile force 
| 
] | ee 
D/P) 0 a sin wt asin wt sipsin (wt — 3 — 5 
| ‘ Hie ae son) ; 
Lit . s sin wt 0 & sin [wt —A— 4 ssin wt 
| | | 
| a | E m2 4 
F | sksin(wt-+ A) am sin (we + = 0 | s Rsin (wt + A) +8" sin (wt — 0) 
| 2 
M  smsin(wt + 7/2) a Bsin (wt + 0) a Bsin(wt +d) +o", sin (wt — A) 0 
7 | wt | 9 vet ( yt m> aa | 2 yr m? vs 
Wie | &R | 2B | B ae) | 8 (2 — 8 ) 


The storage and loss parts of the extensional stiff- 
ness are 


F =F + m? (B’/B?), 
Ye a, BB Z 
P2=F24F", tand=F’/F’. 


(33) 


(34) 


We conclude that a periodic tensile force causes 
a vibration composed of an extension and a bend- 
ing of the same frequency. The phase lag between 
these vibrations equals 6+2/2 and the ratio of 
their amplitudes equals m/B (determined by 
eq. (18) and (23)). 

From eq. (33) we may compare the stiffness 
under pure tensile force and the stiffness in pure 
extension. The storage part is larger in the case of 
vibration under pure tensile force, whilst the loss 


10' 


Se ee 


0.1(m/B) 
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Fig.5. The characteristic functions in vibration under 
pure tensile force for a two-layer beam accord- 
ing to formulae (31), (33) and (34); the mate- 
rial constants are chosen as in Fig.2 (R’, R” 
storage and loss part of the extensional stiff- 
ness, R ratio between the amplitudes of longi- 
tudinal force and extension, 4 phase lag be- 
tween longitudinal force and extension, m/B 
ratio between amplitudes of deflection and ex- 
tension). 


part, the energy dissipation and the loss tangent are 
larger in pure extension. 

The course of the characteristic functions is 
shown in Fig. 5). 

The results for the four basic experiments are 
compiled in Table II. Note that the expressions of 
column 3 and 4 may be obtained by a superposition 
of those of column | and 2. 


5. Comparison with the results of other authors; 
Application to the problem of anti-vibration 


damping 


As already mentioned in the introduction, the 
same problem has been treated earlier by van Oorr 
and Oserrsr. Both authors did not realize the exi- 
stence of cross components between bending and 
extension. Accordingly they did not distinguish 
between pure bending and vibration under pure 
bending moment. Both authors assumed the case 
of pure bending and calculated the corresponding 
bending moment without paying attention to the 
required tensile force. Therefore their treatments 
should conform to case 4.2. and the expressions 
(22) for the bending stiffness 1°. 


Van OortH'! calculated the absolute bending 
stiffness without taking into account the effect of 
damping. Accordingly he obtains an approximation 


B=A,+A,, 


10 Unfortunately they used their results partially for 
the description of the forced bending vibrations of a 
two layer cantilever beam. Such vibrations, however, 
will take place under conditions of zero longitudinal 
force rather than zero extension, and case 4.3 would 
be more appropriate for this problem. 

11 W.P. van Oort, Transactions of Instruments and 
Measurement Conference, Stockholm 1952, Appendix I, 
formula (11) and appendix II, formula (8). 


ee 
where Aj=d, E, E dy Fei) 
A, = dy, Ey E dq? +1r9(r9 + dp) - 


As may be seen from eq. (13) the exact formulae 
in the same notation would read 


B’ =A, cos 6, + Ay cos 05, 
B =/B?+B", 
Bis a A, sin 4 + As sin ds e 


(22 a) 


Van Oort’s expression is therefore valid, as long 
as both loss ratios tan 6, and tan 0, may be neglect- 
ed against unity. 

For the calculation of the damping van Oorr 
starts with a consideration of the energy dissipation. 
Consequently he obtains the exact formula in com- 
plete agreement with eq. (23): 


A, sin 0, + Ag sin 0g 


tan 0 £, 
A, cos 0, + Ag cos 09 


(23 a) 


The theory of Oxsersr’? refers chiefly to the 
problem of antivibration damping, where the lamel- 
lated bar or plate consists of a constructional mate- 
rial [1] of high modulus and low damping and of 
a damping layer [2] of low modulus and high 
damping. This implies a high damping ratio and a 
small modulus ratio #°. 


a=E,'/E,) <1, b=tanod,/tand,>1. (35) 


Under these conditions we may simplify formulae 
(22) considerably: 


12 H. Oserst, Acustica 2 [1952], Beiheft 4, AB 181 
—AB194; formulae (13) and (14). 

13 Note that these conditions do not claim anything 
about the order of magnitude of the product ab, which 
may be well comparable with unity in many cases. 
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_ Ey dPl+2aF(2+3E4+28) +a (36) 
12 (1+ aé) ; 

TE dP ltabe(3+6F+4P+ 208+? &) 
so ae2 (l+aé)2 


B’ 


(37) 
(38) 


(140bE(34+66+484208+07 &) 
(1+4@€) [1+206(2 +325 2) fae 


-evtan Oo; 


The application of these formulae to the problem of 
antivibration damping is postponed to a separate 
treatment. 

In his calculations Oserst assumed the damping 
of the constructional material to be zero, i.e. 
tan 0, =0. This does not effect formula (36) which 
agrees exactly with the expression found by Oserrsr. 
However, instead of (38) he obtained the approxi- 
mation (39) 

34+6644F4208 +a % 


L+aé 14+206(24+36 428) +0’ 


This expression may be derived from (38) by 
neglecting in the numerator | against 3ab &. Con- 
sequently eq. (39) for the damping is valid in cases 
where 


Baan ae 


i.e. where the thickness of the damping layer is 
sufficiently high. It fails for thin damping layers. 
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MEASUREMENT OF SOUND TRANSMISSION THROUGH 
Poon h IN A*DUCT WITH AN-APPLICATION 


TO A RESONATOR 
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Summary 


Magnitude and phase values of both the reflected and transmitted (plane) sound waves 
at an orifice in a brass plate terminated by a square input duct and a non-reflecting output 
duct were carefully measured at numerous frequencies and at an intensity level at which 
linearity is valid. 

A resonator was formed by terminating portion of the output duct with a brass piston, 
and a formula for the reflection of sound from this resonator is derived and proved to be 
valid by experiments. 


Sommaire 


On a mesuré l’amplitude et la phase des ondes sonores planes incidente et réfléchie, au 
droit d’un orifice pratiqué dans une plaque de bronze séparant un tuyau d’amenée 4 section 
carrée d’un tuyau de départ non réfléchissant. Les mesures ont été faites avec soin, pour 
de nombreuses valeurs de la fréquence, et a des niveaux d’intensité ot les phénoménes sont 
linéaires. 

Le tuyau de sortie a été ensuite fermé par un piston de bronze, de maniére A former 
un résonateur: on établit une formule donnant la réflexion du son par ce résonateur; la 
validité de cette formule est expérimentalement confirmée. 


Zusammenfassung 


Betrag und Phase der reflektierten und durchgelassenen (ebenen) Schallwellen an einem 
Loch in einer Messingplatte, die von einem quadratischen Eingangsrohr und einem nicht- 
reflektierenden Ausgangsrohr abgeschlossen ist, sind bei vielen Frequenzen und im linearen 
Intensitatsbereich sorgfaltig gemessen worden. 

Durch Abschluf eines Teiles des Ausgangsrohres mit einem Messingkolben wird ein 
Resonator gebildet; eine Formel fiir die Schallreflexion von diesem Resonator wird her- 
geleitet und ihre Giiltigkeit durch Messungen bestatigt. 


1. Introduction 


The classic theory of the propagation of sound 
waves through simple apertures in unlimited plane 
walls was presented by Raytercu [1] when he de- 
veloped a theory of acoustical (HetmHotTz) resona- 
tors. His well-known result was that an end correc- 
tion has to be added to the length of an orifice 
when the resonance frequency is computed from the 
volume and orifice dimensions of such resonators. 

The more general problem of the complex impe- 
dance of an orifice has been treated theoretically 
by Cranpatt [2]. Efforts of other workers to verify 
experimentally the theoretical results. obtained by 
Rayxercu and Cranpa.i were unsuccessful, and there 
is still no general agreement as to the extent of the 
differences between theory and experimental results. 
A critical discussion of these earlier experiments 
can be found in the introduction to Bour, LaBare 
and Incarp’s paper [3], which presents a precise 
experimental study of the acoustic reactance of 
small circular orifices in thin walls. These workers 


derived an empirical formula for the effective orifice 
radius provided that the intensities of the sound 
fields are such that linearity is valid, as will also be 
assumed throughout this paper. 

The resistive component of an orifice impedance 
has been determined by Srvian [4] and there are 
more recent investigations on related subjects by 
Incarp [5], Incerstev and Nietsen [6], Ferrero 
and Sacerpote [7] and others (see references in 
[3]). 

The sound transmission through an orifice is a 
fundamental problem of considerable practical in- 
terest, e.g. in the design of acoustical resonators 
or perforated sound absorbers, and the results of 
the various workers suggest that further precision 
measurements are very desirable. 

The purpose of this paper is to record the results 
obtained when applying an experimental technique, 
described by the author [8], to measure the propa- 
gation of sound through a circular orifice in a brass 
plate which is terminated on either side by a square 
input or output duct, the output duct being non- 
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reflecting. Magnitude as well as phase values of the 


input and output sound fields have been measured 


carefully and used to derive characteristic quantities 
for specifying the discontinuity in the duct formed 
by the orifice. Interesting new conclusions are drawn 
and a simplified theory of certain types of resona- 
tors! in a duct is derived; this is shown by experi- 
ment to be valid in a wide frequency range. 


2. Measuring technique 


The method used successfully and with good ac- 
curacy by the present author in other discontinuities 
and filters in acoustical ducts [9], [10] has been 
applied to a duct discontinuity consisting of a cir- 
cular orifice in a brass plate. 

Fig. 1 shows the position and measurements of 
the orifice used in the experiments. The orifice is a 
circular opening of diameter 2.07 cm in a brass 
plate 0.635 cm thick. These particular measurements 
of the orifice were chosen because they were con- 
sidered favourable for obtaining good experimental 
conditions and accuracy of final results. The brass 
plate is terminated at one side by the input duct and 
at the other side by the output duct, both of square 
cross-section 7.6 cm 7.6 cm and with their axes 
perpendicular to the brass plate and coincident with 
the axis of the cylindrical orifice. The plate and the 
terminating duct ends with angle pieces have been 
carefully machined to provide a tight fit when bolted 
together. The duct walls made of brass 0.95 cm thick 
are rigid and free of unwanted vibrations and re- 
sonances, the input duct being 170 cm long and the 
output duct being 160 em long, followed by the non- 
reflecting terminal 123 cm long. The non-reflecting 
terminal and its measured frequency characteristic 
of the reflection factor (about 2 to 4 per cent, in a 
wide frequency range) have been described in a 
previous paper [9]. 

A hole (0.6 cm diameter) is drilled through the 
lower part of the brass plate, and a brass tube 
4.1 cm long is pressed into this hole. The inner dia- 
meter of the tube is such that the probe tube can 
slide in it easily. The probe tube (steel, 3.5 mm_in- 
ner and 5 mm outer diameter) is about 300 cm long 
between microphone and probe tube orifice, and has 
an extension piece of methyl methacrylate (5 mm 
diameter) about 180 cm long. The function of this 
extension piece is to seal the passage tube when the 
sound field in the input is probed, and to guide the 
probe tube through the passage tube when the out- 
put sound field is explored. (In the latter case the 
extension piece can slide along the bottom of most 
of the length of the non-reflecting terminal.) 


1 These resonators consist of an orifice plate ter- 
minated by a highly reflecting duct of which the length 
is greater than the cross-sectional dimensions. 
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Fig 1. Orifice plate terminated by the square input and 
output ducts showing position of orifice, input 
and output reference planes, the passage tube 
and probe tube with extension piece. (Measure- 
ments are in cm.) 


A check was made to determine if any appreciable 
leakage of sound occurred through the passage tube, 
leading to a reduction in measuring accuracy. For 
this purpose the orifice in the brass plate was sealed 
by a well fitting cylindrical brass piece and the in- 
put and output sound fields were probed. The en- 
couraging result was that in the frequency range 
from about 300 to 2000 c/s the magnitude of the 
transmitted sound wave was very small. The magni- 
tude ratio of the transmitted sound wave in the non- 
reflecting ouput duct and of the incident sound wave 
in the input duct was 0.01 in most of the measured 
frequency range, with slightly higher ratios (up to 
0.03) at the lowest and highest frequencies. These 
measurements showed also that the pattern of the 
standing wave in the input duct was regular, i. e. 
showing maxima and minima envelope curves that 
can be expected from sound attenuation in the duct 
[11]. This means that unwanted vibrations of the 
probe tube and similar disturbing effects are very 
small, so that it is possible to correct the measure- 
ments for attenuation effects in the input duct. 

After the plugging disc had been removed from 
the orifice a number of tentative measurements of 
the characteristic reflection and transmission factors 
were made using the method described elsewhere 
[8]. These measurements were made at various fre- 
quencies in the range from 300 to 2000 c/s with the 
lowest sound intensity in the duct which was suf- 
ficient to ensure good experimental conditions, This 
level was specified by determining the current 
through the loudspeaker coil at each frequency. The 
measurements were then repeated by changing the 
coil current by factors 2, 4 and 10 corresponding 
to increases of 6, 12 and 20 dB. There were no 
changes of the measured characteristic factors at 
increases of 6 and 12 dB, and only slight changes 
at an increase of 20 dB. These results were taken 
as evidence that measurements at the lower intensi- 
ties were free from non-linear effects. 
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3. The measured characteristic factors 


In Fig. 2 the measured magnitude and phase values 
of the characteristic reflection factor 7 R = A exp (j 2) 
and of the characteristic transmission factor 


T = B exp(j /) 


are plotted against 0/x=2a//, where 2 is the wave- 
~ length in the input and output duct and a=7.6 cm 
the width of the square duct cross-section. The input 
and output reference planes are the surfaces of the 
orifice plate facing the input and output ducts res- 
pectively, as shown in Fig. 1. 


0 02 04 06 08 10 
Ola—¢ —— 


Fig.2. (a) Measured magnitude values 4(@) and B(x) 
and (b) phase values a(@) and f(x) of the 
characteristic reflection and transmission factors 
of the orifice plate of Fig. 1 plotted against 0/x. 


Several independent sets of measurements of the 
reflection and transmission factors at various fre- 
quencies were made by the method previously de- 
scribed [8], and the results of one set are illustrated 
in Fig. 2 as smooth curves of A, B, a and f plotted 
against 0/x. The room temperature during the meas- 
urements was 15 °C + 2°. 

From Fig. 2a it can be seen that A? + B® is close 
to unity (0.93 — 0.97), indicating that the orifice is 
nearly loss-free. In earlier investigations on sym- 
metrical and nearly loss-free discontinuities it was 
shown that a— B= +2x/2 [9], [12], [13]. In the 


2 The characteristic reflection factor is by definition 
the complex ratio of the sound pressure amplitudes of 
the reflected and incident waves at the input plane. 
The characteristic transmission factor is the complex 
ratio of the sound pressure amplitude of the transmitted 
wave at the output plane and of the incident wave at 
the input plane, when the discontinuity is terminated 
by a non-reflecting duct [8]. 
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present case is will be noted from Fig. 2b that the 
measured difference between the phase values of the 
characteristic factors is also very closely equal to 
m/2. In fact, the following equation is fulfilled with 
good accuracy: 


B=a-—(n/2). (1) 


4. The reflection and transmission factors for 
completely reflecting termination (resonators) 


Two interesting problems arise, viz. (i) how can 
measured values of A, B, « and ( be used to find an 
analogous impedance circuit for the orifice, and (ii) 
how can the sound reflection and transmission be 
described in cases of other than non-reflecting ter- 
mination? A special case of the latter problem will 
be considered first. 
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Fig. 3. Cross-section of the resonator with the orifice 
plate, portion of output duct and terminating 
piston. 


The output duct behind the orifice plate was ter- 
minated by a well fitting brass piston at a distance 
1) =30.4 cm from the output reference plane as 
shown in Fig. 3, the whole device representing a 
resonator. The probe tube with extension piece 
(Fig. 1) was replaced by a shorter one, the passage 
tube was removed and the hole was plugged by a 
cylindrical brass disc and the sound field in the 
input duct was probed at numerous frequencies. 
From the pattern of the standing waves in the input 
duct the reflection factor Ru = Ay exp(j %) was de- 
termined in the usual way [8], [10]. In Fig. 4 
the measured values of Ay and ay are plotted against 
@/x showing three successive resonances after the 
principal resonance which is below the measured 
frequency range. It has been shown previously [10] 
that these measured frequency characteristics can be 
computed if the values of the characteristic factors 
of the orifice plate are known. 

The distance J) =30.4.cm being great compared 
with the width of the duct to ensure a plane output 
wave at the piston, the characteristic reflection and 
transmission factors 


R=Aexp(ja) and T=Bexp(j /) 


of a symmetrical discontinuity determine the modi- 
fied reflection and transmission factors 


Ry = Ay exp(j %) and Ty = By exp(j Pu), 


as has been derived in the earlier paper [10]: 
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Fig. 4. (a) Measured magnitude values Aq (Xx) and 
computed values of Ay and By ( ) plotted 
against 0/x, (b) measured (xX) and computed 
( ) phase values ay for the resonator of 
Fig. 3 plotted against 0/n. 


eTrex (—janh/t) 

1—o Rexp(—j4al)/A) ’ 

a ZT 
1—oRexp(—j4xh/A)’ 


Ry=R+4 


(2) 


Ty (3) 


where Q is the reflection factor at the reflection re- 
ference plane (piston). In the special case consider- 
ed above where @=1 and B=a-—(n/2), eq. (2) 
yields: 


» B? 42 A? 2 Acos(a—4n1)/d) 
Ay= |/ 47+B 
c V <3 1+A4?—2 Acos(a—42])/A) 


and (4) 
pewigns Riise sin(a —4 7 1,/A) 
Paik cos(a—4 2 1)/A) S 
Breathe sin(a—4 7 1,/A) 
“a —cos(a—421)/d) 
and eq. (3) yields: 
By 2 (6) 


V14+42—2Acos(¢—4.21)/A) 
and 


sin(a—421)/A) 


baie (a —4 2 1)/2) 


pu =o- Stare tan 


(7) 


With eqs. (4) and (5) and the values of A, B and a 
from the measured curves of Fig. 2 (taking J) = 30.4 
cm) Ay and ay have been computed and plotted 
against 0/x in Fig. 4. These computed curves for 
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Ay and especially for ay agree very well with the 
measurements, thus justifying the procedure adopted 


_ above. Figure 4 indicates that eqs. (4) and (5) are 


also valid in the lower frequency range, e. g. at the 
principal resonance, which in the particular example 
can be found by extrapolation to occur at 0/1~0.07 
which is below the measured frequency range. The 
validity of eqs. (4) and (5) being established by 
the measurements suggests that the related eqs. (6) 
and (7) are valid with a similar degree of accuracy. 
These equations are therefore useful for finding the 
sound pressure amplitude within the resonator at 
some distance from the orifice plate. For the case 
illustrated in Fig. 3 the values of By (which is the 
ratio of the magnitude of the plane sound pressure 
wave at the output reference plane propagated along 
the axis of the resonator and in the direction of the 
piston to that of the incident plane wave in the 
input duct) have been computed with eq. (6) and 
plotted against 0/x (Fig. 4a). There are resonances 
of By and the sound pressure amplitude close to 
the brass piston being 2 By (because of pressure 
duplication at a hard surface) it can be seen that 
at the three resonances the pressure amplitude at 
the piston is from 6 to 11 times greater than the 
amplitude of the incident wave. 


5. Electrical circuit analogy for an orifice 
impedance in a duct 


The other problem mentioned at the beginning 
of the preceding section, viz., how to find an ana- 
logous impedance circuit for the orifice, will now 
be discussed. Using A, B and a with f=a— (x/2) 
for specifying the orifice discontinuity, it has been 
shown elsewhere [8], [10] that the symmetrical 
discontinuity consisting of the orifice plate between 
the input and output reference planes can be re- 
presented by an analogous four-terminal electrical 
circuit, e.g. by the lattice type circuit of Fig. 5 
between, transmission lines. The impedances of the 
two lattice arms are Z, and Z, and the characteristic 
impedance of the transmission lines is Zy). In the 
analogous acoustical case the characteristic impe- 
dance of the input and output duct is 


Zy = c/a? 2 (8) 


where Q) density of medium (air), c velocity of 
sound, @ width of square duct. 


Z 


Fig. 5. 

Lattice type circuit (ana- 
logous to the orifice plate 2) 2 

of Fig.1) with the im- 

pedances Z; and Z, of the 4 2 
two arms. 2; 


AE a i 
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It has also been shown previously [8] that the 
impedance ratios 


Z,/Zo=2=2ir+j23, 2Z2/Zo=22=2R+j 223 (9) 


are simply related to the characteristic factors R 
and T of a discontinuity 


_ 14+(R-T) 


: Sak 1+ (R+T) 
1 1—(R—T) ? 2 


eee (10) 


These formulae yield in the present case with 
B=a— (n/2) (11) 
; 1— A?— BP +j 2(Asina+B cosa) 
4=2%4R+)24y5= = 
eee eile 4 B= 8 (Acosa—Bsin a)’ 
1 — A? — B? +j 2(Asina—B cosa) 
1+ A? + B?-—2(Acosa+Bsina) ° 


22=ZR +] 225= 


—— 
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Fig.6. (a) and (b) Real and imaginary components 
of the impedance ratios z,=2Z,/Z)=z;p4+j 213 
and 2z.=Z./Z)=zop+j 2; of the lattice type 
circuit of Fig.5. (Zp is the characteristic im- 
pedance of input and output ducts.) 


With these formulae the real (zjp, z2R) and imagi- 
nary (213, Z23) components of the impedance ratios 
have been computed using A, B and a (Fig. 2). 
They are plotted against 0/x in Fig. 6. It can be 
seen that z, has a positive imaginary component 21) 
and a relatively small real component z;p, i.e. 


1 - . 53D 
ZR~ 56 21d - The impedance ratio z, has also positive 


imaginary (zgy) and real (z22) components but the 
magnitude of z, is great compared with unity and is 
about 20 times larger than that of z,. This means 
that in the lattice type circuit the two impedances 
Z, can be omitted as a first approximation and Z, 
can be simplified by Z;~j 21; Zy). Fig. 6 shows that 
z13 is roughly proportional to 0/x. As 6/x is pro- 
portional to the (circular) frequency w=c0/a, the 
value of z;; is also proportional to @ so that the 
analogous impedance of the orifice is approximately 
a mass reactance (j#M) with 


mM 221520_ 2215 Go 
@) Aa 


(12) 


From this formula and the data in Fig. 6 an averag- 
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ed mass value M,,.=0.64 x 10° g-cm 4 has been 
found, taking 0) =0.00123 g-cm~? (air density at 
15°C and 760mm Hg [14]). A calculated value 
M omp. = 9.68 x 107? g-cm™~ is obtained by apply- 
ing the formulae derived by Bott, Lasare and In- 
carp [3]. The difference between these two results 
is only 6 per cent and they comparse well as far as 
the averaged mass reactance of Z, is concerned. 
(The value of M from the classic theory is in this 
particular case Mass. = 0.87 x 107? g-cm™*4 which 
is about 36 per cent too high.) 


6. Discussion 


The impedance of an orifice discontinuity in a 
duct is usually represented by an analogous series 
impedance between transmission lines. In this case 
it is assumed that the reflected and transmitted 
waves are inter-dependent (R+7=1 see eq. (10) 
for z)= co) so that only one set has to be measured, 
for instance the standing wave in the input duct, 
and the series impedance can then be specified in 
terms of A and a. 

The experimental method adopted here has been 
extended to measuring the transmitted (plane) waves 
in the output duct as well. If the inter-dependence 
mentioned above is not taken for granted a priort, 
an evaluation of these measurements requires the 
series impedance (2Z,) to be replaced by a four- 
terminal circuit which can be regarded as a lattice 
type circuit with two arm impedances Z; and Z, or 
with the corresponding impedance ratios z, and z.. 

Z, and Z, of the lattice type circuit analogous to 
the orifice discontinuity were specified in terms of 
A, B, a (and f=a— (x/2)), yielding the satisfac- 
tory result the measured orifice discontinuity is, as 
usually, very nearly equivalent to a series impedance 
(2Z,) as |Z,|>|Z,| and |Z,|> Zp. 

The experimental method adopted here is con- 
sidered to be advantageous from several reasons. 
Firstly, measurements of both the reflected and 
transmitted waves can be evaluated, paying due 
consideration to each set of measurements which 
may improve the accuracy. Secondly, @ priori as- 
sumptions about an inter-dependence between input 
and output measurements do not have to be postu- 
lated, and this may be very useful in extreme cases, 
e.g. when the orifice plate is very thick. Finally, 
Zp and Zep must both be positive. This is a sensitive 
test of the reliability of experimental results. 

To improve the accuracy of the measurements, 
smoothed curves representing measurements at many 
frequencies have to be used. There are two possibili- 
ties, either to smooth the curves for the charac- 
teristic factors or those for the various components 
of the impedance ratios. As the curves in Fig. 6 
appear to be more complex than those in Fig. 2 it 
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was considered preferable to smooth the curves for 
the characteristic factors as long as the accuracy of 


the single measurement cannot be improved further. | 


It may be mentioned here that an accuracy greater 
than that achieved in the present investigation would 
most likely reveal that equation (1) has to be modi- 
fied slightly. 

In the preceding section it was shown that the 
measured values of z1; in Fig. 6a agree well with 
theory. It is interesting to make a similar compari- 
son for the values of z;p, an accurate measurement 
of which is far more difficult. Applying Sivran’s for- 
mula (eq. (4) in[4)]) to the orifice discontinuity meas- 
ured above would yield values of z;;~0.042 VO/x, 
but the actual values of z;z in Fig. 6a are several 
times greater and not directly proportional to /0/z. 
Satisfactory checks have been made concerning non- 
linear effects and the substantial difference cannot 
be explained solely by inaccuracies in the present 
measurements, (e. g. by disturbances caused by the 
probe tube) or by the particular analogy chosen 
(the lattice type circuit). (It can be shown that the 
values of zjz determined from A and a of Fig. 2 
for an analogous circuit like that of Fig. 5 with the 
two arm impedances omitted are still of the same 
order as those shown in Fig. 6a.) The close agree- 
ment between computed and measured curves in 
Fig. 4 is further evidence that the values of A, B 
and a, and thus the derived values of z;p, will not 
be very inaccurate, unless there are hidden syste- 
matic errors. 

There are maxima and minima of zjp in Fig. 6a 
and it seems possible that these may be explained 
by Incarp’s method [5]. Unfortunately the diagrams 
given, in his paper (Fig. 24) cannot be used to find 
zip of the orifice discontinuity under discussion. 

All measurements reported in the present paper 
have been deliberately restricted to a frequency 
range at which only the principal mode or plane 
waves can be propagated along the duct and higher 
order modes are confined to the vicinity of the ori- 
fice plate. This makes it possible to measure plane 
(reflected and transmitted) waves in the duct only 
and to specify the orifice discontinuity by the cor- 
responding characteristic factors. The actual sound 
fields close to the orifice plate are in general rather 
complex because of the superposition of the higher 
order modes [5] and have not been investigated in 
detail by the method adopted here. This is not ne- 
cessary in many practical cases in which the overall 
effect of these higher order modes on the magnitude 
and phase values of the reflected and transmitted 
waves is all that matters. An example of this was 
provided by measurements of the reflection factor 
of the resonator of which the reflecting terminating 
duct was long enough (/)=4 a) to ensure that at all 
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frequencies practically nothing but plane waves were 
reflected at the piston. 

The derivation of the resonator formulae (eqs. (4) 
to (7)) illustrates the usefulness of the proposed 
method, and it is of practical interest to note that 
the complex general theory in such cases can be 
greatly simplified. 


7. Conclusions 


It is shown in a practical example that magnitude 
and phase values of both the reflected and trans- 
mitted plane waves at an orifice discontinuity in a 
duct can be measured accurately and at such low 
intensity level that linearity is valid. 

The reflection and transmission of plane waves at 
certain resonators in a duct—consisting of the ori- 
fice plate terminated by a sufficiently long duct 
portion and a hard piston —can be conveniently re- 
presented by the characteristic factors. 

An analogy for the impedance of an orifice dis- 
continuity in a duct is set up in form of a lattice- 
type circuit, the various impedances of which are 
computed and discussed for the particular example. 


(Received February 13th, 1957.) 
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BEERS, TO THESE DITORS 


SIMPLE CALIBRATION OF SMALL VIBRATION PICK-UPS 


In the frequency range 8 to 800 c/s the calibration 
of vibration pick-ups is rather simple as is well known. 


_ For example there are available electrodynamic vibra- 


tors for which the curient-force constant can be deter- 
mined accurately. When the total mass of vibrator and 
pick-up and moreover the exciting current are known 
the resulting acceleration can be calculated. In practice 
this method is very useful for instance when measure- 
ments must be carried out in the field: before starting 
the measurements the complete apparatus consisting for 
example of pick-up, cathode follower, amplifier, filters 
and decibel-meter may be checked. Some of the many 
types of vibration pick-ups are known to have a constant 
sensitivity over a wide frequency band. A check at one 
or two low frequencies is then sufficient for many ap- 
plications. Nevertheless an accurate check at higher 
frequencies is sometimes necessary, the more so when 


_ measurements must be carried out in a frequency range 


where resonance of the pick-up may occur. In any case 
there is a need for simple reliable calibration methods 
for the range 0.8 to 16kc/s; there are already a few. 
It seems, however, worthwhile to describe a new method 
that has been in use for some time in our laboratory. 
Essentially the pick-up to be calibrated is compared 
with an accurately calibrated condenser microphone. 
The absolute calibration of this microphone does not 
meet with difficulties whatsoever. On top of a vertical 
travelling-wave tube (length 2 m, inner diameter 4.6 cm, 
with oc-termination) a vibrating table (copper disk, 
thickness 1 cm, diameter 6cm) is mounted. An elec- 
trodynamic vibrator excites this table vertically. Between 
the table and the upper tube flange is a small rubber 
ring of, say, four millimeters thickness (same inner dia- 
meter as tube) so as to enable the table to vibrate 
without appreciably exciting the tube wall. It will be 
clear that a travelling wave in air is sent down the tube. 


About half way down the tube the microphone (mem- 
brane mounted flush with tube wall) measures the 
sound pressure of this wave. On the table the pick-up 
is mounted preferably outside the tube. In this way 
the pick-up may be calibrated directly in dB velocity 
(a better designation might be “rapidity”) level rel. 
4.6°10~§ m/s (for: p=ocv, and in air 46 nm/s cor- 
responds to 20 uN/m?). 

Above 4.3 ke/s hydrogen instead of air must be used. 
This is to prevent the higher modes of sound propa- 
gation, that might be excited by an unsymmetrically 
vibrating table, to travel down the tube. If the pick-up 
is sufficiently insensitive to tangential or rotational vi- 
brations the calibration is not affected by not exactly 
vertical vibrations of the table. Up to 16kc/s small 
vibration pick-ups may easily be calibrated in this way, 
the accuracy being estimated as better than +1.5 dB 
Comparison in the common (overlapping) frequency 
range 200 to 800 c/s) of this method and of the “elec- 
trodynamic” calibration shows differences not greater 
than 0.5 dB. The measured sound pressure in the tube 
is always slightly higher than would be expected at 
first sight. This is due to the extra “pumping effect” of 
the small rubber ring between the table and the flange. 
Correction for this effect is very simple as rubber is 
essentially incompressible. It is thought that the method 
lends itself very well for automatic calibration of vibra- 
tion pick-ups. 
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ULTRASONIC ABSORPTION IN BONES 


The absorption of ultrasonic waves in bones has 
been investigated by Hiirer, mainly as to its frequency 
dependence [1]. We have measured further the tem- 
perature change of the absorption coefficient. The ex- 
periment has been performed by using the pulse me- 
thod, the sample (bones of horses) being placed bet- 
ween the transmitting and receiving probes in a water 
bath. The main results are as follows. 


1. Variation of absorption in different specimens 


The density varied markedly with the position of 
the bones, 9 ~1.5...2, even when they are cut out 
from the same specimen, subsequently the absorption 
coefficients showed a similar wide variation. Fig. 1 
shows the transmission loss (dB) vs. thickness of the 
bones at three frequencies 1.43 Mc/s, 2.86 Mc/s and 
4.5 Me/s. In this figure, the three absorption curves 
cut the ordinate at a common point A, and the inter- 
cept 0A (=54dB) represents the reflection loss. This 
value of 5 dB is in accord with the value calculated 


from the value of the density (9 =1.9) and the velo- 
city (V=3700m/s), the latter being obtained from 
the angle of total reflection for the longitudinal wave. 
The average absorption coefficients of 80dB/cm at 
4.5 Mc/s, 40...50dB/cm at 2.86 Mc/s and 22 dB/cm 
at 1.43 Mc/s are obtained from this figure. This result 
indicates that the absorption coefficient of the ultra- 
sonic waves in tissues is proportional to the frequency 
in agreement with the data by Hiirer. Measurements 
were performed also at 0.7 Mc/s and the linear depen- 
dency of absorption on frequency seems to be valid also 
in the frequency range below 1 Mc/s. The accuracy of 
the measurement in this range, however, is not satis- 
factory, because of the disturbance by multiple reflection 
and interference phenomena in samples becoming con- 
spicuous with the decrease of absorption. 


2. The temperature dependence of the absorption coefficient 


Fig. 2 shows the temperature dependence of ultra- 
sonic absorption in bones at 4.5 Mc/s, 2.86 Me/s and 
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1.43 Me/s respectively, the thickness of the bones being 
indicated in the figure. These values include the re- 


flection loss, so the true values of the transmission loss ° 


in the sample are about 5dB lower. These measure- 
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Fig. 1. Ultrasonic attenuation vs. thickness in bones at 
various frequencies; L] 4.5 Mc/s, A 2.86 Mc/s, 
© 1.43 Me/s. 


ments indicate that the temperature dependence of 
ultrasonic absorption is nearly the same in the fre- 
quency range of our experiment, the absorption in- 
creasing about 25% when the temperature changes 
from 0° to 60° C, though the temperature coefficient 
at 4.5 Mc/s is somewhat lower than those at lower fre- 
quencies. 
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Fig. 2. ‘Temperature dependence of ultrasonic attenua- 
tion in bones at various frequencies; ] 4.5 Me/s — 
(thickness 3.3mm), A 2.86 Me/s (thickness 
4.9mm), O 1.43 Me/s (thickness 6.7 mm). 0 


3. Conclusion 


As indicated above, the absorption of ultrasonic ~ 
waves in bones is proportional to the frequency and the ~ 
increase with temperature rise is only small, as com- 
pared with various artificial high polymers. These facts 
permit us to conclude that the mechanism of absorption 
in bones has not to do with the inter- or intramole- 
cular relaxation processes, but may be accounted for | } 


by a hysteresis type loss. 
(Received February 15th, 1958.) 
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Buchbesprechung 


J. Marauscuex, Einfiihrung in die Ultra- 
schalltechnik. VEB Verlag Technik, Ber- 
lin 1957, Din A5, 536 Seiten, 349 Bilder, 22 Ta- 
feln, Ganzleinen DM 36,—. 

Kaum ein anderes Gebiet der Technik ist in seiner 
schnellen Entwicklung von so viel Erwartungen beglei- 
tet gewesen wie die Ultraschalltechnik. Nicht ohne 
Grund sah man in den hochfrequenten Ultraschall- 
schwingungen ein ideales Mittel, die Struktur der Ma- 
terie zu beobachten und zu beeinflussen, d. h. einerseits 
Mef- und Beobachtungsmethoden und andererseits Be- 
*irbeitungsverfahren zu gewinnen. Viele zu hoch ge- 
spannte Hoffnungen wurden enttaéuscht, daneben wurde 
aber eine Fiille von Anwendungen des Ultraschalls in 
den verschiedenen Gebieten der Technik verwirklicht. 
Zu diesen gehdren: Ortungsverfahren in Schiffahrt und 
Fischerei, Materialfehler-Suche in der Fertigungstech- 
nik, Dickenmessung in der Walztechnik, Material- 
bearbeitung harter Werkstoffe, Emulgierung und Mole- 
kiilabbau in der chemischen Industrie, Therapie in der 
Medizin. 

Der Verfasser des neu erschienenen Buches hat sich 
in verdienstvoller Weise die Aufgabe gestellt, eine 
systematische Darstellung des Gesamtgebietes zu geben. 
Ausgehend von den physikalischen Grundlagen: den 
Eigenschaften der Schwingungen und Wellen kleiner 
Amplitude (Kap. 2) und der Wellen endlicher Ampli- 


(s 
}. 
tude (3) behandelt er die technischen Probleme: die — 
mechanische Ultraschallerzeugung (4), die elektro- 
mechanisch-aktiven Stoffe (5) und allgemein die Ultra- 
schallsender (6) und -empfanger (7) sowie die Mef- 
technik und-die Abbildungsverfahren (8). Den letzten — 
Teil des Buches bilden die Anwendungen des Ultra- — 
schalls in Ortungs-, Me8- und Untersuchungsyerfahren — 
(9) und zur Materialverinderung in Metallurgie, Fer- 
tigung, Nahrungsmittelindustrie, Textilindustrie, Che- 
mie, Pharmazie und Medizin. ' 
Auf allen diesen Gebieten werden die wesentlichen 
Dinge dargestellt und praktische Kenntnisse vermittelt; 
es wird aber vermieden, zu weit ins Einzelne zu gehen. 
DafS dadurch an einigen Stellen spezielle Wiinsche offen 
bleiben miissen, ist verstaindlich und zum Teil auch auf — 
die schnelle Entwicklung zuriickzufiihren. Besonders 
klar dargestellt ist die Systematik der Schallwandler. 
Im ganzen gesehen, wurde die gestellte Aufgabe in — 
sehr befriedigender Weise erfiillt. Dazu tragen auch die — 
vielen Diagramme und Tafeln sowie die Berechnungs- 
beispiele bei. Besonders verdienstvoll ist die konse- — 
quente Durchfiihrung des VAms-Mafsystems. Das 
Buch kann allen auf dem Gebiet der Ultraschalltechnik — 
arbeitenden Ingenieuren, Technikern, Physikern und 
den Wissenschaftlern aus den angrenzenden Gebieten, 
aber auch den Studierenden sehr empfohlen werden. 


K. Tamm. 


